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THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
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In what follows we first discuss the basic code structure, including 
the overall program flow and a brief description of all subroutines. This 
is followed by instructions on the preparation of input data, definitions 
of key Fortran variables, sample input and output, and a complete listing 
of the code. 


3. BASIC CODE STRUCTURE 
3.1 Overall Subroutine Flow Chart^ 


Fig. 3.1 shows the overall program flow, which is divided into input 
(SI) and integration (Ml) routines. Input is via cards or from a restart 
file (RS) which is automatically created as the calculation proceeds. If 
an output flowfield file is being created for input to the A.R.A.P. radia- 
tion code (STARAD), LU will either set up the file for starting a run or 
read the file which already contains flowfield information, for a restart 
run. IF, II, and LI are used for inputting and processing inviscid flow- 
field data maps, while IP establishes initial profiles of velocity, temper- 
ature, etc. for program-calculated profiles. IN prints all input data. 

In the integration routine, S3 solves the finite difference equations, 
utilizing output from VI for the turbulent viscosity, CC for the chemical 
reaction rates, and SL to invert the matrix formulated in the implicit 
solution of the species continuity equation. OT is the output routine and 
all the resetting of variables, and step size controls are performed in M2. 

EN calculates the mass entrained at each integration step while DS calculates 
the displacement thickness and position the "effective plume boundary" 
for use in NASA/LRC boattail pressure calculations. TK and LP are inter- 
polation routines while CP is used to save common and the flowfield files 
for use in restarting the program. 

Additional details on these subroutines are given below. 


3.2 Subroutine Description 

BOATCC Uses input on chemical reaction rates from 52 to set up appro- 

priate terms in the matrix, which is to be inverted for solution 
of the species continuity equations. 

BOATCP Saves common and flowfield files (if created) for use in 

restarting calculations. 


^All subroutine names contain BOAT followed by two alpha-numeric quanti- 
fiers. The BOAT precedent is deleted in the descriptive paragraphs below. 
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Figure 3.1 BOAT code subroutine flow chart 
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BOATDS 
BOATEN 
BOAT IF 

BOATII 


BOATIN 

BOATIP 


Calculates "effective" plume shape (due to jet entrainment 
effects) for use in external inviscid subsonic/transonic 
flowfieid calculations. 

Calculates shear layer growth rates from entrainment rules 
yielding values of (PSII) and (PSIE) at the new 
station, x + Ax. 

Generates vector arrays VJET (J, K, L) and VEXT (0, K, L), 
containing the user-specified inviscid exhaust solutions, 
in mapped streamline coordinates. The J index refers to the 
dependent variable (1 = Y,2 = P, 3*T, and 4 = U); the K 
index to the radial grid point (for the jet exhaust, K * 1 
is the axis and K = KMAXJ is the plume interface; for the 
external flow, K = 1 is the plume interface and K = KMAXE 
is an arbitrary upper boundary); the L index refers to the 
axial station XJET(L) for the jet and XEXT(L) for the external 
flow. 

Yields properties from the mapped vector arrays VJET and VEXT 
by interpolati ve procedures. BOATII has the calling sequence 
BOATII (ITYP, XX, PSIX, IMAXV, KMAXV, IV, VECT, XV, PSV, VI, 
V2, V3, V4, PSJ, FID, NREC). ITYP = 1 or 2 indicates that a 
standard jet or external flow interpolation for a local 
p roperty will be made, while ITYP = 3 or 4 indicates that an 
interpolation for axia l gradients of the jet or external flow 
variables will be made. 

XX is the value of X at which the properties are desired and 
PSIX is the value of the streamfunction, 4;. 

IMAXV is the total number of jet (I JET) or external flow 
UEXT) stations and KMAX the number of processed data points 
at these stations (KJET or KEXT). 

IV is the index of the inviscid mapped station such that 
XV(IV-l) < XX <_ XV(IV) where XV is either XJET or XEXT. 

VECT is either VJET or VEXT while PSV is either the jet 
exhaust mass flow, 4 j(PSJET), or the external flow value of 
4;e(PSEXT). 

VI, V2, V3, and V4 are Y, P, T, and U respectively. 

Prints all pertinent initial data, the chemical reaction 
mechanism, etc. 

Calculates initial shear layer or boundary layer profiles 
if the user does not specify initial data profiles. 


- 4 - 


BOATLI Interpolates in mapped vector arrays for property values 

BOATLP Interpolates for thermodynamic properties at the local 

temperature 

BOATLU Used in setting up or reading flowfield files for input 

to the STARAD code (only for NRAD = 1) 

BOATMI Controls the overall integration process via the following 

sequence of calls and operations: 

0 The auxiliary dependent variables CPBAR, HSTAT, and Y 
are determined. 

0 The turbulent viscosity is determined via a call to 
BOATVI. 

0 The allowable step size, Ax, is established. 

0 Edge conditions at the new station are determined 
from the inviscid data map via calls to BOATII. 

0 The integration procedure is initiated via a call 
to B0ATS2. 

0 BOATOT is called to print flowfield profiles at 
the user specified print intervals. 

0 The program is terminated if the axial station 
exceeds XMAX. 

0 The main dependent variable arrays are reset into 
the initial profile locations via a call to B0ATM2. 

0 Run and job times are compared with user specified 
times and the above sequence is repeated if no time 
limits are exceeded. 

B0ATM2 Primary functions are to perform step size checks and to reset 

dependent variables after an integration step has been taken. 
The following specific operations are performed: 

0 A check is per^'ormed to determine whether the temperature 
change along each streamline has exceeded the maximum 
allowable change, TCONT. If so, the step-size. Ax, is 
halved and the integration process is repeated. 

0 A check is performed to determine whether any mole 
fractions become negative in the integration step. If 
so, the integration process is repeated with a halved 
step-size as above. 

0 The mass entrainment for the next integration step is 
determined from the newly calculated profiles via a 
call to BOATEN. 

0 The calculated dependent variables RU, RT, and RALPHA 
are reset into the initial arrays U, T, and ALPHA in 
an interval extended by the mass to be entrained in 
the next integration step. 
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BOATOT 

BOATRS 

BOATSL 

BOATS! 


BOATS2 


BOATS3 


BOATVI 


Prints out properties at al! radial points at print inter- 
vals, PRNT or PRNTXC. 


Reads in cornnon and flowfield file (for STARAD radiation 
code input), via call to BOATLU, for restarting a calcu- 
lation. 


Solves the system of linear equations generated by the 
implicit chemistry calculational procedure using a Gauss- 
Gordon reduction algorithm with diagonal pivot strategy. 

Main initialization routine which establishes profiles of 
the dependent variables U, T, and ALPHA in evenly spaced 
streamf unction coordinates; also reads and 
thermochemical and inviscid flowfield data. The following 
sequence of operations is performed; 

0 Parameters controlling the run type, grid resolution, 
print interval, etc., and array of the dependent 
variables and/or edge conditions at the initial 
station are read. 

0 Inviscid flowfield data are read 

mapped evenly spaced arrays via a call to BOATIK 

0 If initial profiles are not read, a call to BOATIP 
yields either shear layer or boundary layer starting 
profiles. 

0 The dependent variable arrays are recast into evenly 
spaced arrays in streamfunction coordinates. 

0 BOATIN is called, which prints pertinent initial data 
and lists the chemical reactions considered in the 
calculation. 

Together with the subsidiary subroutines B0ATCC,B0ATEF, 
and BOATTK, this routine comprises the chemical integration 
package. B0ATS2 additionally calls B0ATS3 for each grid 
point to integrate the flowfield equations. 

Contains the finite-difference formulation of the axU^^ 
energy; and species diffusion equations. B0ATS3 integrates the 

and%nergy equations by an 'US"/“1aU?na 
species diffusion equation by an implicit procedure, calling 
BOATSL to solve the resulting system of linear equations. 

Calculates the turbulent viscosity. XMU, for allsrid points 
using either: (1) the Prandtl Mixing Length Model with duel 

length scale provisions for velocity maxima or "’inima at 
interior points, (2) the 2-equation, ke2, turbulent kinetic 
energy model or (3) the Donaldson/Gray eddy viscosity formu- 
lation. 
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4. INPUT DATA PREPARATION 


Col umn 

Fortran Name 

CARD 1 

1 

(11) 

ITYPE 

* 0 restart, * 1 new run 

3-8 

(3A2) 

IFNAM 

restart file name 

10-15 

(3A2) 

NAMAS 

flowfield file name - used for 
input to STARAD radiation code 

17-26 

(FIO.O) 

RTMAX 

overall run time (minutes) 

27-36 

(FIO.O) 

RTJOB 

job time (minutes) 

Column 

Fortran Name 

CARD 2 

1-72 

(18A4) 

TITLE(I) 

job identification 

Col umn 

Fortran Name 

CARD 3 

1-5 

(15) 

MPSI 

number of radial data points in user 
specified initial profile (IDELP = 1); 
number of points for run in other 
initialization options (IDELP =0,-1) 
(maximum of 50) 

6-10 

(15) 

NMPSI 

number of points for run if initial 
profile is user specified (maximum of 
50); for other options set NMPSI = MPSI 

11-15 

(15) 

NS 

number of gaseous species (maximum of 
25) 
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Column 

Fortran Name 

CARD 3 (continued) 

16-20 

(15) 

NR 

number of chemical reactions (maximum 
of 25) 

21-25 

(15) 

NT 

number of temperatures at which thermo- 
dynamic data (Card 15) are defined 
(usually, NT « 22, maximum of 30) 

26-30 

(15) 

IDELP 

indicator for specifying initial 
radial profiles , , * j 

= 0; shear layer profile calculated 

internally 



= 1; user specified profile 



»-l; boundary layer profiles calculated 
internally 

31-35 

(15) 

IPRESS 

inviscid structure indicator 
= 0; constant pressure mixing 

= 2; BOAT overlaid on inviscid solution; 
inviscid property data maps of 
plume and external flow must be 
input on Cards 17 and 18 

36-40 

(15) 

IVIS 

turbulence model indicator 
* 0; Prandtl Mixing Length Model 


=1; Donal dson/Gray Model 



«=-l; kc2 two-equation model (initial 
turbulent kinetic energy profile 
calculated internally) 

e-2; kc2 two-equation model (initial 
turbulent kinetic energy profile 
specified on Card 13) 

41-45 

(15) 

IMAXJ 

number of axial stations input for jet 
exhaust inviscid data map (maximum of 
50) 

46-50 

(15) 

KMAXO 

riu:nber of mapped radial stations desired 
in jet map (maximum of 25) 
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Col umn 

Fortran Name 

CARD 3 (continued) 

51-55 

(15) 

IMAXE 

number of axial stations input for 
external flow inviscid data map 
(maximum of 50) 

56-60 

(15) 

KMAXE 

number of mapped radial stations 
desired in external flow map (maximum 
of 25) 

61-65 

(15) 

lOUTl 

chemical production terms (W) output 
indicator 
= 0 no output 

= 1 * terms output for each species 

66-70 

(15) 

I0UT2 

production/depletion (RP/RM) terms 
output indicator 
= 0 no output 

= 1 RP/RM output for each reaction 

71-75 

(15) 

NRAD 

flag to generate radiation output 
tape 

= 0 no output file for radiation code 

= 1 flowfield output file generated 
for input to STARAD radiation cod< 
(file name on Card 1 ) 

Column 

Fortran Name 

CARD 4 

1-10 

(E10.3) 

X 

initial axial station (ft), cannot he 
0 for IDELP » 0, typically X = 0.1 W 

11-20 

(E10.3) 

RJ 

nozzle exit radius (ft) 

21-30 

(E10.3) 

XMAX 

total length of run (ft) 

31-40 

(E10.3) 

PRliT 

print interval (ft) 
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Column 

41-50 

(E10.3) 

Fortran Name 
XCHANG 

CARD 4 (continued) 

change print interval at this axial 
location (ft) 

51-60 

(E10.3) 

PRNTXC 

new print interval (ft) 

61-70 

(E10.3) 

FDL 

multiplies program calculated step size. 
Ax, in order to reduce step sue. Useful 
in initial regions with steep 
(e.Q., initial boundary layers); typical 
value for initial boundary layers; = 

0.2, to supress oscillations. For smooth 
initial profiles set FDL = 1.0 

71-80 

(E10.3) 

DFDL 

at each step. DFDL is added to FDL until 
FDL = 1.0. In problems ^Hh initial 
boundary layers, FDL = .2 and DFDL .05 
should prove adequate 

Col umn 

1-10 

(E10.3) 

Fortran Name 
XLE(l) 

CARD 5 

turbulent Lewis number 

11-20 

(E10.3) 

SIGMA(l) 

turbulent Prandtl number 

21-30 

(E10.3) 

TCONT 

maximum allowable temperature change 
Demitted in an integration step ( K). 
typically. 5” < TCONT < 10”K 

31-40 

(E10.3) 

TKINET 

chemical kinetics cut-off temoerature- 
chemistry assumed frozen below this value. 
(If TKINET » 0, the default value of 
400°K will be used.) 

41-50 

(E10.3) 

CARBON 

use only if writing output ^o file for use 
in radiation calculations with solid carbon 
in Dlume. If CARBON « 1.0. program wi n 
calcul^e normalized 

of inert species for input vO radiation code 
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Column 

Fortran Name 

CARD 5 (continued) 

51-60 

(E10.3) 

CNZINT 

mole fraction of N£ at jet exit (only 
needed if CARBON • 1.0). If using this 
option, set mole fraction of N£ in free 

stream, Xw. ■ .78973 
"2,e 

61-70 

(E10.3) 

CVISC 

multiplies values of turbulent viscosity 
at initial station (default, CVISC • 1.0) 

Column 

Fortran Name 

CARD 6"^ 

1-10 

(E10.3) 

P 

pressure (atm) for constant pressure 
mixing solution 

11-20 

(E10.3) 

U{1) 

jet velocity (ft/sec) 

21-30 

(E10.3) 

U(MPSI) 

external flow velocity (ft/sec) 

31-40 

(E10.3) 

T(l) 

Jet exhaust temperature (®K) 

41-50 

(E10.3) 

T(MPSI} 

external stream temperature (°K) 

Column 

Fortran Name 

CARD 7 (must always input this card even 

if none of the parameters are used) 

1-10 

(E10.3) 

FFF 

ratio of t/6 in Prandtl Mixing Length 
Model in nearfield shear layer region; 
use FFF « .065; must also be input for 
kc2 initialization procedure 

11-20 

(E10.3) 

GGG 

ratio of £/6 in Mixing Length Model in 
fully developed region; use GGG ■ .08 


+In overlaid procedure (I PRESS ■ 2) these values are redundant. 
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Column 


Fortran Name 


CARD 7 (continued) 


21-30 PSID 

(E10.3) 


Input PSID - 1.0 if "effective" plune 
boundary Is to be calculated. Can only 
be used If IPRESS ■ 2; PSID must be set 
« 0 If IPRESS - 0 


31-40 

(E10.3) 


DELJ 


jet side boundary ^ayer displacement 
thickness at nozzle exit plane (ft) 


41-50 

(E10.3) 


DELE 


external boundary layer displacement 
thickness at nozzle exit plane (ft) 


51-60 

(E10.3) 


USTJ 


jet side frictional velocity ratio 
(default value of 1/30 built In) 


61-70 

(E10.3) 


USTE 


external frictional velocity ratio 
(default value of 1/30 built in) 


71-80 

(E10.3) 


RBUOY 


bouyancy indicator; set RBUOY ■ 1.0 to 
include buoyancy term in momentum equation. 


*** Cards 8 and 9 are required only if IDELP < 0 (i.e., they are not 
required for a user specified initial proTile) 


Column 

Fortr.^n Name 

CARO 8 

1-10 

(E10.3) 

ALPHA(l.l) 

irole fraction of first jet exhaust species 

11-20 

(E10.3) 

ALPHA(2,1) 

2nd species, etc. Total of eight species 
per card 

Column 

Fortran Name 

CARO 9 

1-10 

(E10.3) 

ALPHAd.MPSI) 

mole fraction of first external stream 
species 

11-20 

(E10.3) 

ALPHA(2,MPSI) 

2'*d species, etc. 
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Colunffi Fortran Name CARD 10 


1-10 

(tl0.3( 

RIN(l) 

nondimensional radial location (r/RJ) of 
first grid point out of MPSI user specified 
points. This point can be the lower edge 
of a shear layer. Do not input RIN(l) ■ 0 
for axis, start witFRlNTl) « .01 

11-20 

(E10.3) 

RIN{2) 

radial location of 2nd grid point, etc. 



continue with a total of eight values per 
card 

Column 

Fortran Name 

CARD 11 

1-10 

(E10.3) 

T(l) 

temperature of first grid point ("K) 

11-20 

T(2) 

temperature at 2nd grid point (°K), etc. 
eight values per card 

Column 

Fortian Name 

CARD 12 

1-10 

(E10.3) 

U(l) 

velocity at first grid point (ft/sec) 

11-20 
( no . 3) 

U(2) 

velocity at 2nd grid point (ft/sec), etc., 
eight values per card 

*** “ -2 (i-e.. when the TKE option is 

selected with a k.nown initial turbulent kinetic energy profile) 

Column 

Fortran Name 

CARD 13 

1-10 

(E10.3) 

XK(1) 

turbulent kinetic energy at first orid 
point (ftVsecM 

11-20 

(E10.3) 

XK(2) 

turbulent kinetic energy at 2nd grid point 
(ft'/sec*), etc., eight values per card 

Column 

Fortran Name 

CARD 14.1 

1-10 

(E10.3) 

ALPHA(l.l) 

mole fraction of 1^ species at 1^ point 
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Col urm Fortran Name 


CARD 14.1 (continued) 


n-20 ALPHA(2,1) 

(E10.3) 


mole fraction of 2nd species at 1^ point; 
continue to NS species, eight per card 


Column Fortran Name 


-10 ALPHA(1.2) 

.10.3) 


CARD 14.2 


mole fraction of Ut species at 2nd point, 
etc. 

continue to Card 14.MPSI in a similar 
manner 


*** The i-.ext group of cards contain the thermodynamic data. The order of 
the S£ecj|_es must be consistent with that on Cards 9 or 14. FoFiach 
specie: , the first card contains its name, irolecular weiaht and heat 
Of formation; the second and subsequent cards contain the temperature. 

the^JANNAF tables enthalpy, input exactly as presented in 


Column 

Fortran Name 

CARD 15.1.1 

1-4 

(A4) 

AID(l) 

name of 1st species (H 2 O, CO 2 , etc.) 

7-16 

(E10.3) 

WTMOLE(l) 

molecular weight of 1st species 

17-26 

(E10.3) 

HF(1) 

heat of formation - kcal/mole) 

Column 

Fortran Name 

CARD 15.1.2 

1-10 
(Flo. 4) 

TTB(1) 

temperature of 1st species at 1st data 
point (®K) — 

11-20 
(Flo. 4) 

CPTB(l) 

Cp of 1^ species at 1st point, (cal/mol 

21-30 

(F10.4) 

GTB(l) 

Gibbs free energy of 1st species at 1st 
point -(F°-H° 298 )/T. I^/mole-‘’K) — 
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Column 

31-40 

(F10.4) 


41-50 

(F10.4) 


51-60 

(F10.4) 


61-70 

(F10.4) 


Fortran Name 

CARD 15.1.2 (continued) 


HTB(l) 

static enthalpy of 1st species at 
point, (H®-H® 293 ). T^al/mole) 

1st 

TTB(2) 

temperature of 1st species at 2nd 

point 

CPTB(2) 

r 


GTB(2) 

J same as above at 2nd point 



71-80 HTB(2) 

(F10.4) 




Column Fortran Name 


CARD 15.1.3 


1-10 TTB(3) 

(F10.4) 


temperature of 1st species at 3H data 
point, etc. This is continued until the 
data at all ^ points (input on Card 3) 
is specified for the 1st species 


Column Fortran Name 


CARD 15.2.1 


1-4 AID(2) 

(A4) 


name of second species, etc. Repeat the 
sequence of cards 15.1 .1— -15.1 . (NT/2) 
for each of the NS species 


*** The next group of cards contain the chemical reaction mechanism, the 
reaction type indicator and associated rate coefficient data. The 
order here is arbitrary. 


Col umn Fortran Name CARD 16.1 - First Reaction 

1-6 ZID(l) soecies A 

(A6) 

7 + sign 
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Column Fortran Name CARD 16.1 (continued) 


8-13 

(A6) 

ZID(2) 

species B (or M) 

14 

+ sign 

(if needed) 

15-20 

(A6) 

ZID(3) 

(Blank or M) 

21 

= sign 


n-11 

(A6) 

ZID(3 or 4) 

species C 

28 

+ sign 

(if needed) 

29-34 

(A6) 

ZID(4 or 5) 

species D (or M) 

35 

+ sign 

(if needed) 

36-41 

(A6) 

Blank or ZID(5) 

species E (or M) 

42-48 

leave blank 


49-50 

(12) 

IRR(l) 

reaction type, 1 to 10 (see NASA CR-3075) 

51 

(ID 

IRT(l) 

rate coefficient type, 1 to 8 (see NASA 
CR-3075) 

52-59 

(E8.2) 

RC(l.l) 

pre-exponential factor. A, (cm-molecule-s 
units); note that kf = AT"N exp (B/RT) 

60-63 

(F4.1) 

RC(1,2) 

temperature exponent, N 

64-72 

(F9.1) 

RC(1,3) 

activation energy, B (cal/mole) 
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CARD 16.2 - Second Reaction 


repeat procedure of Card 16.1 for each reaction 

**♦ The next group of cards comprise the inviscid flowfield data and 
are required only if IPRESS = 2 


Column 

Fortran Name 

CARD 17.1.1 

1-2 

(12) 

LMAP 

number of data points for 1st jet station 
(maximum of 25) 

6-15 

(E10.3) 

XJET(l) 

axial location of 1st jet data station (ft) 

Col umn 

Fortran Name 

CARD 17.1.2 

1-10 

(E10.3) 

VJET(1,1,1) 

radial location (ft) of 1st point at 1st 
station (must be axis point) 

11-20 

(E10.3) 

VJET(1.2,1) 

radial location (ft) of 2nd point at Ijt 
station, etc. 

Column 

Fortran Name 

CARD 17.1.3 

1-10 

(E10.3) 

VJET(2,1,1) 

pressure (atm) at 1st point 

11-20 

(E10.3) 

V0ET(2.2.1) 

pressure (atm) at 2nd point, etc. 

Col umn 

Fortran Name 

CARD 17.1.4 

1-10 

(E10.3) 

VJET(3,1,1) 

temperature (®K) at 1st point 

11-20 

(E10.3) 

VJET(3,2,1) 

temperature (°K) at 2nd point, etc. 
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Col umn 

Fortran Name 

CARD 17.1.5 

1-10 

(E10.3) 

VJET(4,1,1) 

velocity (ft/sec) at 1st point 

11-20 

(E10.3) 

VJET(4,2,1) 

velocity (ft/sec) at 2nd point, etc. 

Col umn 

Fortran Name 

CARD 17.2.1 

1-2 

(12) 

LMAP 

number of data points for 2nd jet 
station 

6-15 

(E10.3) 

XJET(2) 

axial location of 2nd jet station (ft) 


Continue supplying jet data in this sequence for 


all IMAXJ station? 

1 • 

Col umn 

Fortran Name 

CARD 18.1.1 

1-2 

(12) 

LMAP 

number of data points for 1st external 
flow station (maximum of 25) 

6-15 

(E10.3) 

XEXT(l) 

axial location of 1st external stream 
data station (ft) 

r.ol umn 

Fortran Name 

CARD 18.1.2 

1-10 

(E10.3) 

VEXTd.l.l) 

radial location (ft) of 1st. point at 
1st station (must be at mviscid plume 


interface) 

11-20 

(E10.3) 

VEXT{1,2,1) 

radial location (ft) of 2nd point, etc. 

Col umn 

Fortran Name 

CARD 18.1 .3 

1-10 

(E10.3) 

VEXT(2,1,1) 

pressure (atm) at 1st point, etc. 
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Col umn 


Fortran Name 


CARD 18.1.4 


1-10 

(E10.3) 


Col umn 

1-10 

(E10.3) 


VEXT(3,1,1) temperature (°K) at 1st point, etc. 


Fortran Name CARD 18.1.5 

VEXT(4,1,1) velocity (ft/sec) at 1st point, etc. 


... Continue supplying external flow data in this 
sequence for all IMAXE stations. 
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5. DEFINITIONS OF KEY FORTRAN VARIABLES 
AND PROGRAM WORKING UNITS 


ALPHA(J,I) 


mole fraction of species J (input at 
grid point I); redefined in BOATS! 
as ALPHA(J.I) = ALPHA(J,I)/WTVR where 
WTVR » mixture molecular weight 


CPBAR 

HSTAT 

P 

PSI 


specific heat of mixture, ft /sec /°K 


static enthalpy of mixture, ft^/sec^ 


pressure, Ib^/ft^ (input in atm) 


mass flow (radial) coordinate, 
(slug/sec )% 


RHO 


density, gm/cm^ 


SIGMA 

t" 


turbulent Prandtl number 
temperature, K° 
velocity, ft/sec 


WD0T(J,I) 


chemical reaction rate of species J at 
grid point I, mole/cm^-sec 


WTMIX 


reciprocal of mixture molecular weight 


axial distance, ft. 


XE' 


turbulent dissipation, (ft/sec)Vsec 


XK' 


t. 


turbulent kinetic energy, ftVsec' 


Note that RALPHA, RT, RU, RXE, RXK are values of these dependent 
variables at the end of an integration step. 
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turbulent Lewis number 


turbulent viscosity, slug/ft- 
radial distance from axis, ft 


6 . 


TEST CASE 


Input and output for a sample case are given in this section. The case 
is a nearly perfectly expanded cold air jet, with a total to exit static 
pressure ratio of 2.0, and is one of the cases analyzed in NASA CR-3075. 


- 22 - 
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6.2 Sample Output 


Some terms that appear on the output and are not defined elsewhere, or 
whose definition is not obvious, are defined below. 

PSIBAR Nondimensional mass flow parameter. 

MASS FLOW PARAMETER for JET, Inviscid map, « PSI (JET) 

for EXT, Inviscid map, ■ PSI (EXT). 

ENTH-TKE If kc2 turbulence option used, this column prints the turbulent 
kinetic energy, XK, in ftVsec^ - for all other options it 
prints static enthalpy in cal/gm. 

PR-XE If ke2 turbulenc j option is used, this column prints turbulent 
dissipation, XE. in ( ft/sec )Vsec - for all other options it 
prints static pressure, in atm. 

XBAR Axial distance/jet radius. 

PSID Value of PSI at dividing streamline. 

PSI(X) Value of PSI at edge of shear layer. 

Y(MPSI) Value of Y at edge of shear layer. 

SLOPE dY/dX at edge of shear layer. 

V/u Ratio of normal to axial velocity at edge of shear layer. 

* 

DEL*(X) Local value of 6 

DEL*(X-DX) 6* at previous station. 

RDIV Radius of viscous dividing streaml ine/jet radius. 

REFF Radius of effective body/jet radius. 
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E)(1* INVISCIO HAP 


STATIUN N'JMRfP 


AXIAL LnCAIlaN • 0. 


MASS FLOW PAP.'iMt fFi» • 


PI 

1 

? 

3 

A 

5 

6 


y 

. l?7A9E*00 
.IP910E *00 
,2A06?E *00 
.28750E*OO 
. 33l«5F*OQ 
. 37A6AE*00 


P 

.908PAE *00 
.90A06F *00 
. P9900F *00 
.89677E *00 
, gq<,<, «ir *00 
.HQIlfet *00 


I 

. 79H26E*03 
.797BIE*03 
,20733E^03 
.7970«E*03 
, ?9691E*03 
.7967HE»03 


U 

, 396AHE*03 
.A089AE *03 
, A? 190E»O3 
,A?80«E *03 
.<.3337E*03 
,A3663E*03 


PSIBAR 

0 . 

.^ 0000 E *00 
. AOOOOC 400 
.60000E*00 
,0OOOOE*OO 
,10000E*01 


SIATIUN NUMBER ? 


AXIAL LHCATIUN * ,f>9AA5E-0l 


MASS FLOW PARAMflF* « 


PI 

1 

?. 

3 

A 

5 

6 


y 

.1?^>0OF*OO 
. 106A9E *00 
.P30?7E*OO 
,?8533E*00 
.3?991E*00 
. 37 ^ 90 L♦OO 


P 

.9031 7F *00 
.90105E *00 
. 69P1 ?E *00 
.096O7F *00 
,09A5AE»OO 
.0933AF 400 


T 

.?9773E»03 
. ? 97 S 7 E»Q 3 
,Z9725E*03 
.?9705E*03 
,29f.9lE*03 
, 2960OE*O3 


U 

.A179SE*03 
,A1 796E*03 
.A7A90E *03 
.A798?F*03 
.A33S0E*03 
,A363BE»03 


PS IBAR 

0 . 

,?0000E*00 

,A0000E*00 

.60000F»00 

,0OOOOF*OO 

,10000E»01 


STATIUN NUMBER 


AXIAL LOCATION « ,1A706E*00 


HASS FLOW PARAMETER » 


PT 

1 

I 

3 

A 

5 

b 


y 

. l?2?3E*00 
. 10271E*OO 
. ?3AA6E*00 
. 781 i» 0 L^OO 
. 3?6?3E*00 
.36930E*OO 


P 

, 0979AF *00 
. P9735F *00 
. 89bl OE *00 
. P9A91 f *00 
.P93P8F *00 
, 09799 F *00 


T 

.<?97?3E*03 

.2971HL*03 

.29706E*03 

.^ 969 AE ♦03 

.79b0SE*O3 

.?9b76£*03 


U 

.A?510E*O3 
.A760OE *03 
.A?90BE*O3 
.A3775E*03 
.A3520E*03 
.A3730t*03 


PSIBAR 

0 . 

,20000E*00 
.AOOO0E*OO 
.6C000E *00 
•60000E*00 
. 10000E*01 


STATION NUMBER A 


AXIAL LOCATION • ,23A30E*OO 


MASS FLOW parameter • 


PI 

1 

P. 

3 

A 

b 


y 

. 12025E*00 
. 1 8000 E *00 
. 732SBE *0J 
.?71bAt*00 
. « 0 !) 
. 7A0f ♦ 00 


P 

,B9hS7F*00 
. B9577L *00 
. BO<, 5 BE *00 
.R9371E *00 
.P930?r *00 
. *00 


T 

.29710E+03 
•29702E*03 
.?9691E*03 
.29603E*O3 
,;'9677E*a3 
. 29671F*03 


U 

.A209SE*O3 
.A 309AE *03 
.A33bAE*03 
,A3^6AE*03 
,A 373 I F *03 
.A3H7HF *03 


PSIBAR 

0 . 

.20000E*00 
. A0000E*00 
.60000E*00 
,800001*00 
,10000E*G1 


. 39175E*00 


.391 A2E*00 


, 390b6E*00 


. 39016E*00 


station 

NUf»8ER 5 

AXIAL 

LOCATION » 

.33333E +00 

MASS FLOW PARAMETFP » .19077E+00 

PT 

Y 

P 

T 

U 

PSIBAR 

1 

.11?22E*00 

.89357F^00 

.29682E+03 

.43607E^03 

0. 

2 

.179b6t*00 

. 89333E *00 

. 29679E +03 

3662t *03 

.20000E*00 

3 

.231 33E*00 

. 89289E ♦OO 

. 29675E^03 

.<•37676 ^03 

.<*0000E^00 

<t 

.2 78<f<VE»Oi) 

,892<i^E400 

.29671E^03 

.<.3d71E^03 

.60000E+00 

5 

.3231^E^OO 

. 89203E ♦OO 

. 29667E+03 

.<• 3968E+03 

.0OOOOE ^00 

6 

.3663 7E + 00 

. 89163E +00 

. 29663E ^03 

.<•<•0596 + 03 

. lOOOOE^Ol 


STATION 

NUMBER 6 

AXIAL 

LOCATION • 

.<»<*6<i3E ♦OO 

NASS FLOW PARAMtTtP . .39097E*00 

PT 

Y 

P 

T 

U 

PSIBAR 

1 

.U835E + 00 

. 0921OE *00 

. 29668E^03 

.<V3952E^03 

0. 

2 

.17801F +00 

.R9193F>00 

. ?9666E^03 

.<•39936 *03 

. 20000E^OO 

3 

.230<.5E*03 

. 0916<.E >00 

.2966<.E^03 

.'♦‘•060E^03 

.<»0000E+00 

<• 

.27756E^OO 

.8013«E+00 

. 29661E^03 

.<»<tl21E^03 

.60000E^0O 

5 

.32225E^03 

. P9l 1 3F +00 

,29659E^03 

.<.<.180E^03 

.BOOOOE^OO 

6 

.36550E^00 

.B90P9E ♦OO 

.29656E^03 

.<*<♦2356 ♦OS 

.lOOOOE^Ol 


STATION 

NUMBER 7 

AXIAL 

LOCATION ■ 

.57692E ♦OO 

MASS FLOW parameter « ,391?9E^00 

PT 

Y 

P 

T 

U 

F SIBAR 

1 

.11791E^00 

. R9089F ♦OO 

. 29656E^03 

.'•'•235£^03 

0. 

2 

. 17835F^00 

.OROeOE^OO 

. 29655E^03 

.<i<i257E^03 

.20000E^00 

3 

. 22997E^00 

. 8906<iE ♦OO 

.2965<^E^03 

.<*<3293E^03 

.<30000E^00 

<• 

.277076^00 

.R90<^9E +00 

.29652E^03 

.'♦<i327E^03 

.60000E^00 

5 

.32179E^00 

. 89035E^00 

.29651E^03 

.<t<V360E^03 

.aooooE^oo 

6 

.36506E^OO 

, 8 <30? IF >00 

.29650E^03 

.^<»393E^03 

.lOOOOE^Ol 


STATION 

NUMBER 8 

AXIAL 

LOCATION ■ 

.35000E ♦Ol 

MASS FLOW PARAMETER « .83569E^00 

PT 

Y 

P 

T 

U 

PSIBAR 

1 

. 12515E^00 

. BP7QAF ♦OO 

. 29628E+03 

.<»<i917E^03 

0. 

2 

.SOOllE^OO 

. 0879AE ♦OO 

.29628E *03 

.<•<•9176 ♦OS 

.20000E^00 

3 

.<3<^55<iE*00 

. 887<J<iF ♦OO 

. 29628E^03 

.<•<391 8E ♦OS 

,<.0000E^00 

<• 

.58392£^0n 

.8P793E^00 

.2962HE^03 

.‘♦‘•91HE^03 

.60000E^00 

5 

. 71916E+C0 

.8P79 3E ♦OO 

,2962HE^03 

.<.<3919E^03 

.eooooE^oo 

6 

.P529HF^03 

. pP703r ♦00 

.29628E^03 

.<•<•91 9E^03 

. loooor ♦Ol 


JtT SIDE BOUNDARY LAYER PROFILE FROH I ■ 1 TO 3 


ENTHALPY 

CP 

MOL WT 

GAS CONST 

GAMMA 

SOUND SPEED 

MACH 

.086326*05 

, 108936*05 

.208506 *02 

.310506*05 

.139056*01 

.101696*05 

. 10999F*01 

R IN 

U 

T 

TERM 

UX 

TX 

DUMA 


OUMB 


.lOOOOE^Ol 
. SOOOOfc^OO 
. lOOOOE-02 


. 11186E»0<i 
. 10037E ♦O'l 

0 . 


.?3B2<.E *03 
• *03 

.?9567E*03 


0 . 

. 71<.95E-01 
.57150E+00 


.11 IBfeE^O^. 

. 1 1 106E*O«f 

.11106l+0<i 


,2382'.E^03 

.2382<»E*03 

.2382‘tE^03 


0. 

. l?07<iE-Ol 
,359MF-01 


0 . 

,3'>7<»eE-Ol 

.32092E*00 


PHYSICAL Y DISTRIOUTIUN 


. 1 I993F+00 
. 12252E*00 
.12511E^00 


EXTERNAL BOUNDARY LAYER PROFILE FROM I 


ENTHALPY 

-.08632E»O<i 

RIN 

. 100001-02 
.22222E-01 
. A<*<»4.^E-01 
.66667E-01 
. 88089E-O1 
. lllllE^OO 
.13333E«-00 
. 18556E*00 
.17778E+00 
. 20000E»00 
.22222E+00 
.2<.'f'i<.E*00 
. 26667E^00 
.28089E+OO 
.31111E»00 
. 33333E*00 
. 35556E+00 
.37778T ♦OO 
. iiOOOOE^OO 
.‘.2222L+00 
L* 00 


CP 

. 10927E*05 
U 

O. 

.25525E^03 
.27993E*03 
.29«.51E»03 
. 30502E «-03 
.3133«.E*03 
. 32031E *03 
. 32636E ^03 
. 33177E*03 
. 33670E *03 
. 3<.127E+03 
. 3->555E^03 
. 3<f962E >03 
. 353«.9E»03 
. 35722E *03 
. 36082F *03 
. 36<» 31E+03 
. 36770E >03 
. 37100E »03 
. 37<i23E*03 
. 377386 »03 


f^OL WT 
.208«.OE»O2 
T 

.305'*9Ff03 
. 30251F *03 
.30191E+0 3 
.30153F*03 
,30l2<.E+03 
.30100E+0 3 
. 3O0POE*O3 
. 30062E403 
.300«.6E*03 
. 30031E»03 
.30017E+03 
. 300036*0 3 
. 299906 *03 
.2997BE*03 
,29966E*03 
.2995«.E*03 
.299<.2F *03 
.299316*03 
.29920E+03 
.299096*03 
.298986*03 


TO 

GAS CONST 
. 310'«0E*0‘. 

TERM 

.571506*00 
. 707066*00 
. 306856*00 
.350806*00 
. 31 1 236*00 
.288986*00 
.270886 *00 
.255536*00 
.25181E*00 
.229596*00 
.218166*00 
.207596*00 
.197616*00 
.180106*00 
.178986*00 
.170196*00 
.161676*00 
. 153396*00 
. 155336*00 
. 137556*00 
. 129766*00 


59 

GAMMA 
. 1 3968E*01 
UX 

.529016*03 
.52901E*03 
.529016*03 
.52901E*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.52901 F*03 
.529016*03 
.52901E*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 
.529016*03 


SOUND SPEED 
.113596*05 
TX 

.29707E*03 
.297076*03 
. 297076*03 
.297076*03 
.297076*0? 
.297076*03 
.297076*03 
.297076*03 
.297076*03 
.2970 76*03 
.297076*03 
.297076*03 
.297076*03 
.297076*03 
. 29707F»O3 
.297076*03 
,297076*03 
.297076*03 
.297076*03 
. 297076 *03 
.2970 76*03 


MACH 

. 370026*00 

DUMA 

0 . 

. 10865P-03 
.388056-03 
.805556-03 
.133956-02 
.197926-02 

.27129P-02 
.35310F-02 
.552506-02 
.538896-02 
,650906-02 
. 758516-02 
,86051P-n? 
.9765T6-02 
.10957^-01 
,121756-01 
.135116-01 
.156606-01 
.15915F-01 
.171706-01 
. 185186-01 


DUMB 

0 . 

. 1 50226-01 
,236196-01 
,31 1926-01 
. 38108E-01 
.555306-01 
.505506-01 
.562266-01 
.616006-01 
.666996-01 
.715586-01 
,761616-01 
.005526-01 
.057326-01 
.00 7096-01 
.925916-01 
.960056-01 
.996936-01 
.102726*00 
. 1 057HE *00 
.108666*00 


. <i6667t»00 
,‘*e089E+OO 
.51111E400 
. 53333F+00 

,5777Ht*00 
,60000E»00 
.6?2?2E^00 
, 6<»^<.<.E+00 
.66667E»00 
.68089E»OO 
.71lllE^00 
.73333E^OO 
. 7^>t>56E + 00 
.77778E*00 
.HOOOOEfOO 

.e2222F>00 

E*00 
,B6667E»00 
.a8«89E*00 
.911 1 lE^OO 
.93333EfOO 
.95556E+0Q 
.97778E+00 
.lOOOOE^Ol 

PHYSICAL Y 


. 380<.5E + 03 
. 383<t6E *03 
. 3e6<iOE *03 
. 3H926E *0 3 
. 39206E»03 
. 39^ 7HE *03 
, 397'*3E*03 
. <.0000E *03 
. <*02^9E *03 
.A0<.89E«-03 
.‘.0721E*03 
,'i09^3F*0 3 
. <.1157E *03 
. <tl 360E *03 
. ^S'tE +03 
.«.1 738E»03 
.«il9llE*03 
,^P07^B*03 
.<»2225E+03 
. <i2366E+03 
.«.2<.96E *03 
.<i261<VE^03 
.<^2721E*03 
.<i2016E*O3 
. <*2901E^03 


.29807E+O3 
.29877F*03 
.29866F403 
.?985*>c»03 
.298<»^)E »03 
.2983ftF»03 
.29A27F403 
. 29T 1 7F403 
.2980PF 40 3 
.29799P403 
.29791 E403 
.29767F403 
.2977<*F 403 
,297h7F403 
.29759F 403 

.29752F 493 

.297<»6E403 
.29739F 403 
.2973<»E403 
.29728F 403 
.29723F403 
.29719E403 
.2971<»F403 
. 29711F 403 
.29707F403 


. 12223E400 
.11<.87E400 
. 10768E400 
.1006<.E400 
.93768E-01 
. 87065E-01 
.80537F-01 
.7A190E-01 
.68030E-01 
.62065E-01 
. 5630<iE-01 
.5075‘*E-01 
.<»5‘i23C-01 
.<i0320F-01 
.35<»52E-Ol 
. 30827E-01 
.26'.51E-01 
.22 332E-01 
. 18<»7<»E-0I 
. 1<*882E-01 
.11 562E-01 
.8bl50E-02 
.57<.<«7E-02 
. 32521E-02 
. 10377E-02 


.^2901 E403 
.^2901E403 
.<i 2901E403 
.'V2901E403 
.<.2901E403 
,<.2901 E+03 
.<.2901E403 
.<i2901E403 
.<i2991E403 
.<.2901E403 
.<»2901E403 
.<*2901 E403 
,<.2901E403 
.<*2901E403 
,<.2901E403 
.<.29010403 
.<.2901 E403 
.<*2901E403 
,<.2901 E403 
.<.29010403 
.<i290iE403 
.<.29010403 
.<.2901 0403 
.-429010403 
.<429010403 


.297070403 
.297070403 
.297070403 
.29707F403 
.2970 70403 
.297070403 
. 29707F 403 
.297070403 
,297070403 
.2970 70403 
.297070403 
.297070403 
.297070403 
.297070403 
.297070403 
.297070403 
.29707E403 
.297070403 
.2970 70403 
.2970 70403 
.297070403 
.297070403 
.297070403 
.29707E403 
.297070403 


.196'55C-01 
.2087<4F-01 
,220F)9F-01 
.23236P-01 
,2<437O0-O1 
,25<.65F-01 
.265190-01 
.275250-01 
.28<.«2F-01 
.293860-01 
. 302330-01 
.31022F-0T 
. 31750F-01 
. 32<.160-O1 
.330200-01 
.335600-01 
,3<.036F-01 
. 36<»690-Ol 
.36800F-01 
. 350900-01 
. 353210-ni 
,35-4950-01 
.356150-01 
.356850-01 
. 357070-01 


.125300400 
.129090400 
.13306E400 
.137030400 
,1<4099E400 
. !<. <4960400 
.158930400 
. 1 52900400 
.156870400 
.160850400 
. 16581E400 
.168770400 
. 1 72750400 
.176710400 
. 180680400 
,185650400 
.188620400 
.192590400 
.196560400 
, 200520400 
.205590400 


.11 1380400 
. 1 1 3930400 
.11 6320400 
.11 8560400 
.1 20650400 
.122580400 
.1 25370400 
.126020400 
.127530400 
.128910400 
.130160400 
.131290400 
. 13230F400 
.13319E400 
.1 33980400 
.135660400 
.11525E400 
.135750400 
.136160400 
.136590400 
.136750400 
.136950400 
.137060 400 
. T 37150400 
.137160400 




.208<V6E*00 

.?16‘i0t»00 
. ?^ 037£♦00 
*00 

.22830F»00 
.23227t»00 
.2362<t6+00 
.2<»02l E»00 
.2<»«tl8F*00 
.2‘*815E*00 
.25212E*00 
. 25608E«-00 
.26005E«-00 
.2b<i02E*00 
,26799E>00 
.27196E*00 
.27593E*00 
«27990E*00 
.20386E»OO 
.287H3E400 
. 29180E+00 
.29577E+00 
.29979E*00 
.30371E+00 


AERONAUTICAL RESEARCH ASSOCIATES OF RRINCETON 
AX I SYM1ETR IC MIXING WITH NON-E OU I L I BR I U« CHFM7STRY 

TEST CASE 1 - WITH INVISCIO FLOW MAPS 


initialization of boat RUN» restart file - , BOAT-SPECRA FUR - 78/09/20. 

PPESSUPF ( INIT I AL ) » .8961563E*00 ATMOSPHERES 


NC7ZLE RADIUS* . 1 2 50000F ♦ DO FEET 
LEWIS NUMBER < CONS T ANT I ♦ .lOOOOOOE+Ol 
X INIT lAL IFEFT) * .lOOOOOOE-01 

PRINT INCREMENT* .5000000E*00 

KE? TURBULENCE MODFL 


BUOYANCY FACTOR * 0, 

PRANOTL NUMBER (CONSTANT I f .lOOOOOOE+01 
X FINALIFEETI* ,3000000P*01 
MINIMUM STEP SIZE* .1009000F-09 


TEMPFPATURF (DEG. KELVIN) 
VELOCITY (FFE I/SECOND ) 
MOLE fraction 02 

MOLE FRACTION N2 


JET EDGE 

.2970733E«^03 
.11I8556E+0A .<*290072E+03 
.2100000E+00 .2100000E*00 
.7900000E*00 .7900000b»00 














«. .IOOOOOOE-01 EFU I F •; T C*^t 1 - WIIH IMVISCIO FlUM 

«/» JFII* « FCFI PRESS!. Tn> 

.90')0')OF-OJ .><7<MlF-03 .w'JblObt^OO 


Pt 

Y/P 

vri IJC 1 1 Y 

TF5PFR*ru«! 

0! NS I T Y 

n*C>! NO. 



F EE T/ St C 

K 

ON/CC 


1 

.9696 

.111H56F‘06 

.738735F *03 

. 1377751-0? 

.1099971 *01 

7 

.9717 

. 106877E *06 

.7650 39E *0 3 

• 17H61 71 -02 

.1016911*01 

3 

.9897 

.5338981*03 

.7 708601*03 

. 1 159 3»!1-0Z 

.6971 761*00 

6 

1 .6768 

.771 1 7<iE *03 

.107176F *03 

. 10 3551 ! -02 

.736953E *00 

5 

1 .0B7r 

.7930661 »03 

.3016791 *03 

. 1058981-07 

.7607311*00 

6 

1.1111 

.317716! *03 

.3010361 *03 

. 105991! -02 

.7733031*00 

7 

1 . 1 769 

. 37706OF *03 

. 300769E *0 3 

. 106)561-07 

.7870651*00 

B 

1.7703 

. 3798l6t*03 

.3005171*03 

. 1061071-07 

.7889561*00 

9 

1.9617 

.336377E*03 

. 3003181*03 

.1061 371-07 

,796798! *00 

10 

1.3015 

.367070E '"I 

.3001671*03 

. 1061 66! -0? 

.7998351*00 

1 1 

1 . 3600 

. 3670/7E *03 

.799983E *03 

. 1061861-07 

. 3063691 *00 

1 7 

1.1773 

. 361 7C5E »03 

.7998351 *03 

.1067011-07 

. 3086831*00 

1 3 

1.6135 

. 3659961 *03 

.7996961*03 

. 1067 151-07 

.3173161 *00 

16 

1 .6688 

. 360006E *03 

.7995661*03 

. 1067751-07 

.3159051*00 

1 6 

1 .6837 

,3637P7E*03 

.7906 391*0 3 

,1067361-07 

.3197901*00 

16 

1 .5169 

. 36737 3F *01 

.7993191*03 

, 1067611-07 

. 3775011*00 

17 

1.5698 

.3 70 78 31 ‘01 

.7992061 *0 3 

. 1067661 -07 

.3755611*00 

1 8 

1.6871 

.3760671 *03 

,7990931 *03 

. 1067501-07 

.3786071*00 

19 

1 .6138 

.377180! *03 

, 7 9898' t *0 3 

, 1067571-07 

.3317961*00 

70 

1 .6650 

. 380I86E *03 

.79888 1 1 *0 3 

. 106751 1-07 

.3339901*00 

7 1 

1 .6756 

. 38 3073E *03 

.7987801 *03 

.1067681-07 

. 3365661*00 

7? 

1.7058 

.306866k*O3 

.798682! *03 

.1067661-07 

.3390831 *00 

73 

1 .7355 

• 383533E *03 

.7985871*03 

.1067391-07 

.3616911*00 

76 

1 .7668 

.391116F*03 

. 7986951 *0 3 

. 1067331-0? 

,36381 31 *00 

75 

I. 7937 

.3936071*03 

.7986061*03 

.1067 76 F -0 7 

.3660511*00 

76 

1 .3773 

.396000! *03 

.7983191 *03 

. 10671 81 -07 

. 3687101 *00 

77 

1.8505 

.3983111*03 

.7987351*03 

. 1067091-07 

. 3507901*00 

?M 

1 ,8786 

• 60053»t *03 

.7991561 *03 

. 1 061 991-07 

. 3577951*00 

79 

1 . 9060 

.6076781*03 

.7980751 *0 3 

. 1061 88F-0? 

. 3567761 *00 

30 

1 .9333 

.606 737E *03 

,7979991 *03 

. 1061 771-07 

.3560781*00 

31 

1 .9603 

.606697F *03 

.6979761*03 

.1091661-07 

.3578501*00 

3? 

1 .9871 

.603687! *03 

.7978561*03 

.1061501-07 

. 3595511*00 

33 

7.0137 

.6103901*03 

. 6977H8F *01 

. 1061 15F -07 

. 3611671 *00 

36 

7 .0600 

. 6 1 7 1 7 1 F ♦ 0 3 

.79777 31 *0 3 

.1061 301-07 

.3677651 *00 

35 

7.0661 

.61 37771 *03 

.7976611 *03 

. 1061 311-07 

. 3667611*00 

36 

7.0970 

, 6 1 5 36 II 1 *03 

.797601 ( *03 

.1061 301-07 

. 3656601*00 



fsiH-tKE vrsmsiTv 

C«l/G*! IR/FT/SFr 

0. 0. 

.32<i6/3E*09 .*3399 Sc_oz 

.1^6H1^E♦05 

.->5B5«9E*03 .311939r_n3 

.<.3fc737E*0!. . ’fe<)7R0F-0 1 

.1 753>73E*0‘. .I7U31F-01 

.106?8F>E»0<i .n3»79r-()l 

.776<i5<iE*03 .nilR^F-Ol 
.61 3«.95E*03 .I0n5l'-01 

.60B09<»E*03 

,6'.6330E*03 .86?7nP-0^ 

.3'79B3<.£*03 .Bl*676F-6^ 

.3h600BE»03 .7B3376C-02 

.339637E*03 ,76'.665F-0Z 

.317B>7BE*03 ,730»36F-0? 

.799900E»03 ,70»3O7C-02 

I 6E *03 ,^667R9F-0^ 

.^70771E♦0^ .673'>76C_r,; 

.^SPIO^E^OS .66B367F-92 
,?66B00f»03 .663616F-0? 

.736969E»03 . fc”71 76F-02 

.7?5567E*03 .61B|70F-0? 

.716576E*03 .BOl’HF-O? 
.705736E*03 .6B76%6r_^2 

.196119E*()3 ,6735'.3r-o? 

.186698F*03 .5*!O9»rf-02 

.177396F*03 .66B363F-07 

. 16B720F »0 3 . 9 3 167 E -r ) 2 

•158616E*03 .S1S676C-62 

.l'.8397E*03 .<1986507-67 

.13BB79E*03 .6<»77 6iir-0 2 

.130115E»03 .6668087-02 

.1716961*03 .651706F-02 

.H 3 309t»03 .6 355367-0? 

.lO6688E*03 .61»763r-07 

.»67?95E»02 ,3607677-62 


PAGE 1 70709/70. 


PSI PR-KE 

.707797E*00 0. 

.7O5?05E»OO .175657E*I0 

.707773E600 .670769E*O9 

.71076?f»00 .1?6B56F»07 

.?17750E*00 .670566E*07 

.715738E*00 .107079E607 

.717726E*00 .506876E*06 

.?70?I6E*00 .315757EF06 

.7Z7707F300 .771617E*06 

.775191E«00 .166876E*0b 

.??7679E*00 .13667IE*0b 

.730167E*00 .U6696E*06 

.737655E600 .10707SE«0b 

,735163E»00 .911717E»05 

.737631E*00 .875B78E*05 

.760170E*00 .756767E*05 

,767608E*00 .698092E*O5 

.765096E»00 .669037E605 

.?675B6E*00 .606917t*05 

.75007ZF.600 .566939E705 

.75756lf»00 .5Z8090E709 

.755069E»00 .693676E*05 

.757537E*00 .661036E»05 

.760075E»00 .629987E*05 

.767513E600 .600769E»05 

.765001E*00 .3717071*05 

.?67690E*00 .366266E*05 

.769978E*00 .317908E*05 

,?77666E*00 .791076E*05 

.776956E*00 .763606E*05 

.777667E*00 .738366£*05 

.?79931E*00 .?16760E*05 

.787619E*00 .19537IE*05 

.7B6907E*00 .175766E*05 

.?B7395E*00 .155677£*05 

.2B98B3E«00 .135607E*05 



37 7.1177 

31 7.1I.37 

39 7.16«b 

<.0 7.1931 

(.1 7.711B 

<.7 7.7. J 7 

13 7.7bH5 

ll 7.7931 

11 7.1177 

1b 7.1171 

17 7,lbb1 

11 7.390b 

19 7.1117 


.11b>*31t *03 
,1|t)719» *01 
.1 I J19 3t *03 
. 1701bot *03 
.1770b‘>l *03 
.17 11 bbf »0 1 
.1717091 *03 
.17'il'lbl *03 
.17b091f ♦03 
. I7b91 bf »03 
.1771711 403 
. 1 7 1 3 7 1 1 ♦ o 3 
.17900 71 *03 


.7971111 403 
.7971904 403 
.7971 39E 40 3 
.7973004 40 3 
. 79 7 311 E 40 3 
.79730U 403 
.’97?bir 4C3 
.’97773E 403 
.797117E403 
.797I11E403 
.7971 7bE 40 3 
.79 7091E 40 3 
.7970731 403 


. 10bl7 lE-07 
. 10b 1 771 -0? 
. lObl 71F -07 
. 1061 701 -07 
.1061 111-07 
.1061 one -07 
. lOhlOlF-07 
.1060931-07 
. 1060111 -07 
. 1060 711-07 
.1060631-0? 
. 1060171 -07 
. 10601 71 -07 


. 3670011 400 
. 3617071 400 
. 36919 71 400 
. 3 706191 400 
.371 7711 400 
. 3777791400 
.3736731 400 
. 3711181 400 
. 371 3101 400 
.3761711400 
. 3 7691 11 400 
. 3771171400 
.3710191 400 


.8691901407 
. 7969101407 
.7770031*07 
.6117711*0? 
.171911 140? 
.1 100911*07 
.1131771 *0? 
.3990171 *0? 
.3317101*07 
.711081 1*07 
.7110831*07 
.7061131407 
0. 


.111in'‘~07 
.3b1’11f-07 
. 31111 7C 07 
. 339’79C-07 
.3101 ooc-07 
.7971617-0? 
.1711M7-0? 
. 7113 » 17_07 
.7310317-07 
.71-3177-0? 
. 70701 77-07 

0 . 


.797371E400 
.7918601400 
.7973181 *00 
.7998 361*00 
.3073711*00 
. 3011 1?E*00 
. 3073011 400 
. 3097891400 
.3177771400 
.3117611400 
.3177131400 
.3197111 400 
.3777301400 


.1187171405 
. 1036611 401 
.9037091 401 
. 7666111 401 
.6317111401 
.1308311401 
.1111611 401 
. 3673171401 
.7877661401 
.7717971401 
.1717011401 
.1367831 401 
0 . 


T f T C « ■> E 1 


UIIM flow 


RAOF 


1 


71/09/?0 


»• . i •?mooof -01 f f r T 


,17 


PI I /» 

til! 


1 

. OvO'.7 

7 

,071'.'. 

1 

. OdOl H 

I. 

i .17*! 7 7 

i 

1.0110/ 

t 

1-11110 

7 

1 . 1 7010 

H 

1 . 7707 7 

0 

1 . 7bl 7 1 

10 

1 . 1015*. 

1 1 

1, l«.t 01 

l.> 

1 . 37770 

1 1 

1 , 1 151 

!«• 

1 , <.<V«70 

!•» 

1.1.1177 

Ife 

1.51617 

1 7 

1 . 55017 

1 H 

1 . 5171 1 

10 

1 . 61 11 1 

70 

1 ,655'10 

71 

1.67561 

77 

1 . 70510 

71 

1 . 71557 

7*. 

1 . 7551 1 

7^) 

1.71175 

7h 

1 .177 10 

77 

1 .1501 1 

7B 

1 .1 7150 

70 

1 .00510 

10 

1 .01 170 

11 

1 .06011 

17 

1.01711 

33 

7. Cl 161 

i>> 

7.010O1 

15 

7 . 06600 

lb 

7.01171 

17 

7.11 750 


I : : 1 

. 71 O'J )l ♦O') 

. 7 1 00')F 400 
.?100 )l ♦ jO 
.7 lOOOf ‘00 
, 7 loom *00 

. 7 lOO'Jt KJU 
.710011 *00 

.7iaoot*uo 

.7 100 11 •OIJ 
.710001 ‘or 
,7 lOOOI *00 
.7rj >oi *00 
.71000F * 00 
.7 10011* JO 
.7 100 If ♦ 10 
.710001 ♦ JO 
.7 1000F *00 
.71 0001 *00 
.7 lOJlt *00 
.7 1 OOOF *00 
.710001 *00 
.7 lOOOf *00 
.710001*00 
. 7 lOOOt *00 
.7 10O0F *00 
. 7 lOOOf *00 
.710001 *00 
.7 lOOOF *00 
,7 100 )t *00 
. 710001 *00 
.710001 *00 
.7 moOl *00 
.710001 *00 
. 7 1 0001 *00 
.7100)1 *00 
.7 10001 *00 
.7100)1 *00 


*1 7 

p I 

. 700001 * O'J 
, 700001 * 00 
, 7OO0O r *00 
.700001 *00 
. 700001 *00 
. 7O0001 • 00 
. 700001 * 00 
. 700001 *00 
. 700001 *00 
. 70000' * 01 
. 700001 * 00 
, 70QOO' *00 
. 7Q0001 *00 
. 700001 *00 
. 7O00Q1 *00 
. 700001*00 
. TOOOOf *00 
. 700001 *00 
. 700001 *00 
. 700001 *00 
. 7OC001 *00 

. 700001 ♦ no 

. 700001*00 
, 70000" • JO 
. 700')01*0'J 
. 700001*00 
. 700001 *00 
. 700101 *00 
. 700001*00 
. 7O0001 *00 
. 700001*00 
.700001 *00 
. 700001*00 
. 700001*01 
. 700001 * 00 
, 700001 *01 
. 701001*0 ) 


0 . 
0 . 

0 . 

(. . 

0 . 
0. 
0 . 

0 . 

0 . 
c . 
0. 
0. 
0 . 
c . 
0, 
0 . 
0 . 
0. 
0 . 
0 . 
c . 
0 . 
0 . 
0. 
0 . 

0 . 
0. 
0 . 
0 . 
0 . 
0 . 
c. 
0 . 
0 . 

1 . 

c . 


dlJLl nON'; 


0 . 

1 . 

0 . 

0 . 

0 . 

0 . 

1 . 
0 . 
1 . 
0 . 

1 . 

u. 

0 

0, 
0 , 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
1 . 
0 . 
1 . 
0 . 
0 . 
1 . 
0 . 
0 . 
0 . 
1 . 
0 . 

0. 
0 . 
T • 
0 . 

1 . 


0. 

0. 

') . 

c. 

1 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0, 
0 . 
0. 
u. 
0 . 
0. 
c. 
0. 
0. 
0. 
0. 
0 . 
1 . 
0. 
0. 
0 . 
0. 

0. 
0 . 

1 . 
0. 
0 . 


0 . 

0 . 

0 . 
0 . 
1 . 

0 . 

1 . 
I). 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 

0 , 
0 . 
0 . 
0 . 
0 . 
0 . 

1 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 
0 . 


O. 

0 . 

0 . 
0 . 
o, 
0 . 
0 . 
0 . 
0 . 
0 . 
O. 
0 . 
0 . 
0 . 
0 . 
O. 
0 . 
0 . 
0 . 
o. 
0 . 
O . 

0. 

1. 
0 , 
0 . 
o , 
O . 
1. 
0 . 
O . 
O. 

o. 

O. 

O. 

o . 

O. 


38 

1 . l*.322 

.2 1000E*00 

. 790OOE^OO 

0. 

39 

2c 16858 

.2 lOOOt *00 

. 7900QE *00 

0 . 

<»0 

2. 19378 

.21000E *00 

. 79000F ♦00 

0. 

‘•I 

2.21882 

.21000fc*00 

. 79000F *00 

0. 


2. 2<*373 

.2100QE400 

. 7Q000E *00 

0 . 

<i3 

2.26850 

.2 lOOOE *00 

.79oonr*oo 

0. 


2.2931‘t 

.2 1000E*00 

, 79000F*00 

0. 


2.31767 

. 2 lOOOE +00 

. 79000F *00 

0 . 

*tb 

2. 3‘^208 

.21000E *00 

. 79000E*i)0 

0 . 

«i7 

2.36639 

.2 lOOOE^OO 

. 79000E *00 

0 . 

<.8 

2 . 39060 

.21000F*00 

. 79000E *00 

0 . 

*>9 

2.<»l‘.72 

. 21000E *00 

. 79000E ♦OO 

0 . 


XBAR 

PS ID 

PS MX) 

Y(8PS I) 

SL OPE 

.200992 

.208717 

.325028 

.303353 

.012500 

.«.612<.3 

. 200797 

. 328328 

. 305585 

.065356 

.773669 

.208797 

. 329<i20 

. 305563 

-.01 7112 

1.087<i7l 

.208797 

.329959 

. 30<.916 

-.017985 

1 .<101920 

. 208 79 7 

.330391 

. 306577 

-.006637 

1.717339 

.208797 

.330777 

. 306380 

-.005355 

2.033670 

.208 79 7 

.33113<i 

. 306109 

-.0081 36 

2.350822 

.208797 

. 331«.71 

. 303778 

-.008503 

2.668025 

. 208797 

.33179*. 

. 303630 

-.0086H6 

2.98 790<t 

.200797 

.332106 

. 303206 

-.003912 

3 . 308306 

. 208797 

.332609 

. 303121 

-.006223 

3.630079 

.208797 

.332705 

. 302966 

-.002606 

3.953363 

.208797 

. 332995 

. 302060 

-.002723 


0. 

0. 

0. 


9. 

0. 

0. 

0. 


9, 

0 . 

0. 

0. 


9. 

0. 

0. 

0. 


9. 

0 . 

0. 

0. 


9, 

0. 

0. 

0. 


9. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 


9, 

0. 

0. 

0. 


9, 

0. 

0. 

0 . 


9, 

0. 

0. 

0. 


0. 

0. 

0. 

0. 


9, 

V/IJ 

0EL*(X) 0EL*(X-DX1 

POIV 

PCPF 

-.252026 

.2O80‘’2 

.221033 

1.006016 

1.700872 

-.062037 

,179937 

.102867 

• • 

1.179937 

-.063163 

,155931 

.158930 

1.015665 

1.155931 

-.032026 

.132773 

135059 

1.016627 

1.132773 

-.017017 

.117688 

. 118676 

1.017166 

1 .117688 

-.016760 

.105708 

.106877 

1.018607 

1 .105700 

-.018777 

.093292 

.096601 

1 .020306 

1 .993292 

-.010602 

.080286 

.001579 

1.021660 

1.080206 

-.018 368 

.067326 

.068603 

1.023659 

1.067326 

-.013257 

.058185 

.059109 

1 .0251 95 

1 .058185 

-.013290 

.068929 

..069856 
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LISTING OF BOAT CODE 


A CMiplete Fortran listing of BOAT, as presently used on the NASA/ 
Langley CDC CYBER series Is given In this section. 
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PROGUAM BUAT I INPUT , OUTPU T, TAPE l» T Ap£2 » I APE i>» T APE6»0UTPUT » 
NAIN PROGRAM FOR BOAT 


BUATCM - 
OIMENS ION 

1 CM(25»25) 

2 HCPC2f25» 

3 ISAVEI6I* 
A RCf A0> 3}» 

5 S1GMAC50I 

6 VJETIA>2;» 

7 XFC50)# X 

8 Y0Un50»# 


BOAT CUHMr)N 

All^Olf AID(29}» AL0CI!»0> 
, CPBAH(50), CPTBV(7!>0)» 

* HSIATI50I# HTBV(750>» 
J123A5(5>» PSEXT(50I» P 
RHrj(50»> RMO0UTI50l» RT 

* STARTIl 1* TI50I* TITLE 
»2)» WD(JTJ?5»50I» WMC25I 
EXTI50>» XJET(50I, XK(50 

ZID(3) 


6)f ALPHAI2 
ECC(501» G 
IRR(AOI» IR 
SI(30)» QX( 
(50»» RUI50 
(18)> U<30) 
» WP(23)> W 
It XLECSOI* 


3t30l» CGHV(730t3lt 
C25l» GTBV<750)t 
RR(AOt3)t IRTIAOI* 
25l» RALPHAI25»501# 
It RXE(30lt RXK(30I 
t VEXr«A,25t2lt 
TMiXCSOIt WTK0LEf25 
XMUI50lt Yt50lt 


LOGICAL LHALFtLSWON 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


IJlt Jl?3A5mi 
lALOCdtllt CMfltlll 
lECCt CN(ltl3llt lYOUT 
(CPTBVt CGHVI 
fHCPdtllt wntiii 
ISTARTt TTRIlll 


t CN(ltI3ll 


BUAT 
BOAT 
BUAT 
BOAT 
BOAT 
BOAT 
BOAT 
BOAT 
BUAT 
t BOAT 
BOAT 
ItBUAT 
BOAT 
BOAT 
BOAT 
BOAT 
BOAT 
BUAT 
BOAT 
BOAT 
BOAT 
BUAT 
BOAT 
BOAT 


COHHON 

TTBI30I 

HFC25I 

f 

CPTBV 

9 

GTBV 


HTBV 



BOAT 

Zb 

COriNON 

A 

AID 

P 

ALPHA 

9 

CARB 


CM 

9 

CPBAR 

tBOAT 

Zb 

1 

CRR 

CVISC 

P 

OELPSI 

9 

OFOL 


OPOX 

9 

OX 

tBOAT 

27 

2 

OXNIH 

FOL 


FFF 

9 

G 


GGG 

9 

HSTAT 

tBOAT 

28 

3 

lOELP 

IFCC 

P 

IFINIS 

9 

lOUT 


lOUTl 

9 

I0UT2 

tBOAT 

29 

A 

IPAGE 

IPRESS 

P 

IRR 

9 

IRRR 


IRT 

9 

ISAVE 

tBOAT 

30 

5 

itflg 

IVIS 

9 

JI 

9 

J2 


J3 

9 


tBOAT 

31 

6 

J5 

LHALF 

P 

NMOD 

9 

NPSI 


hxnpt 

9 

nxNPi 

tBOAT 

32 

7 

NPSI 

NR 

9 

NRAO 

9 

NRAS 


NS 

9 

NT 

tBOAT 

33 

8 

P 

PCNT 

9 

PRNT 

9 

PRNTXC 


PS I 

9 

PSIO 

tBOAT 

3A 

9 

QX 

RALPHA 

P 

RBUOY 

9 

RC 


RHO 

9 

RHOOUT 

tBOAT 

35 

1 

RJ 

RT 

9 

RTACU 

9 

RTJAC 


RTJOB 

9 

RTNAX 

tBOAT 

36 

2 

RU 

RXF 

9 

PXK 

9 

SIGE 


SIGK 

9 

SIGMA 

tBOAT 

37 

3 

T 

TCONT 

9 

TF06E 

9 

TEMRM 


tenrp 

9 

TITLE 

tBOAT 

3B 

*» 

TKINET 

V 

9 

UNIT 

9 

MOOT 


WM 

9 

WP 

tBOAT 

39 

5 

WTNIX 

WTMPLE 

9 

X 

9 

XCHANG 


XD 

9 

XE 

tBOAT 

AO 

6 

XINIT 

XK 

P 

XK2 

9 

XLE 


XMAX 

9 

XMU 

tBOAT 

A1 

7 

Y 

710 









BOAT 

A2 

connoN 

FLJTEXllOOOOIt 

NANASnit 

RPRIU2<*00I 

9 

SPACRf 101 


BOAT 

A3 


BOAT 

BOAT 


C 

C 


JET/EXT FLOW FIELD COMMPN SECTION 


n n n 





C 


connaN 

OFLJ 

P 

DFLE 

P 

IJET , 

1 

KHAXJ 

P 

KNAXF 

P 

NRJET , 

2 

PI 

P 

P2 

P 

P3 » 

3 

VJET 

9 

VFXT 

P 

X.IET , 

CUMI10N 

ENOCH 





END OF 

COMMON 

TU 

BE COPIED 

TO RESTART 


CnW^CM FIO(3»5)> IFNA^ni* LSM0MC16 


C 

CALL BOATSl 
CALL BOATHl 
C 

STOP 

END 








BOAT 

A6 

lEXT 

» IMAXJ » 

IMAXE 

>SOAT 

•*1 

NREXT 

» PSJET , 

PSEXT 

»BOAT 


P<» 

» usrj > 

USTE 

»eOAT 

•*9 

XEXT 



BOAT 

50 




BOAT 

51 




BOAT 

52 




BOAT 

53 

FILE 



BOAT 

59 

) 



BOAT 

55 



BOAT 

56 




BOAT 

57 




BOAT 

5B 




BOAT 

59 




BOAT 

60 




BOAT 

61 




BOAT 

62 




BOAT 

63 



r» o o 


C 


C 

C 


C 


c 



subroutine BOATCC (K» I, I, NJ> NEQ» OSIGN» 


i^nr I 


BOATCC > CH CALCULATIONS • GENERALIZED lALSO OX) 

LOGICAL LSIO 
BOATCn - BOAT CllNNPN 

OINENSIQN A(SO)» AI0(?5lt ALnC(S0»6)> ALPHA ( Z S* SO I » CGHVI7SO>3)> 


1 CM(Z5»Z5)* 

2 HCP«2#2Slf 


CPBAR(50)» 
HSTATI 50), 


CPTBVJ750), ECCI50), G(25), GT9VI750), 
HTQV(750)f IRRI<i>0)* IRRRfAOfSI# IRTIAO)# 


BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 


1 

2 

3 

A 

5 

6 
7 

e 

9 


3 ISAVE16), J123<i5l5), PSEXT(50), PSI(5U), OX(25), R ALPHA 1 25, 50 ) , 

BCC 

10 

A RCC<iO,3), RHJ«50), RHOOUTCSO), RTI50), RUC50), 

RXE (50 ), 

RXK(50), 

BCC 

11 

5 SIGNAI50), STARTfl), T ( 50 ) , 

TITLEdS), U(50), 

VEXT(A,25,2), 

BCC 

12 

6 VJETfA 

,25,2), WD0T<25,50), 

WNI25), 

WPI25), WTN1XI50), WTNQLE f 25 ) , BCC 

13 

7 XEI50) 

, XEXTI50), XJFTI50), 

XK 150) 

, 

XLEI50), XNU(50I, YI50), 

BCC 

lA 

8 Y0UTI5O), ZI0I5) 






BCC 

15 









BCC 

16 

LOGICAL 

LHALF,LSWON 






BCC 

17 









9CC 

IB 

EQUIVALENCE Ul, 

J123A5U)) 






BCC 

19 

EQUIVALENCE CALOCIl,!), CHfl 

,1) ) 





BCC 

20 

EQUIVALENCE (ECC 

, CNI1,13) ), 

lYOUT, 

CN(1,15)) 



BCC 

21 

EQUIVALENCE CCPT8V, CGHV) 






BCC 

22 

EQUIVALENCE (HCPIltl), WNIDi 





BCC 

23 

EQUIVALENCE tSTART, TTBll)) 






BCC 

2A 









BCC 

25 

connoN 

TTBI30) 

HFC25) , 

CPTBV 

» 

GTBV , 

HTBV 


BCC 

26 

CONNON 

A 

AID , 

ALPHA 

f 

CARS , 

CN 

CPBAR 

,3CC 

27 

1 

CRR 

CVISC , 

DELPSl 

f 

DFOL , 

OPOX 

OX 

,BCC 

ZB 

2 

OXNIN 

FCL , 

FFF 

f 

G , 

GGG 

HSTAT 

»BCC 

29 

3 

IDELP 

lECC , 

IFINIS 

P 

lOUT , 

lOUTl 

I0UT2 

,BCC 

30 

t* 

IPAGE 

I PR ESS , 

IRR 

f 

IRRR , 

IRT 

ISAVE 

,BCC 

31 

5 

ITFLG 

I VIS , 

J1 

9 

J2 , 

J3 

JA 

,BCC 

32 

6 

J5 

lhalf , 

N«r(iL» 

9 

NPSI , 

P.XNPT 

NXNPl 

,BCC 

33 

7 

NPSI 

NR , . 

nrao 

9 

NRAS , 

NS 

NT 

,BCC 

3A 

8 

P 

PCNT , 

PRNT 

9 

PRNTXC , 

PSI 

PSIO 

,BCC 

35 

9 

QX 

RALPHA , 

PBUOY 

9 

RC , 

RHO 

RHOOUT 

,BCC 

36 

1 

RJ 

RT , 

RTACU 

9 

RTJAC , 

RTJQB 

RTNAX 

,BCC 

37 

2 

RU 

PXE , 

RXK 

9 

SIGE , 

SIGK 

StGNA 

,BCC 

3B 

3 

T 

TONT , 

TEOGE 

9 

TENRH , 

TENRP 

TITLE 

,BCC 

39 

A 

TK INET 

U , 

UNIT 

9 

WOUT , 

WN 

WP 

,BCC 

AO 

5 

WTNIX 

WTNOLF , 

X 

9 

XCHANG , 

xo 

XE 

,BCC 

A1 

6 

XINIT 

XK , 

XK2 

9 

XLE , 

XNAX 

XNU 

,BCC 

A2 

7 

Y 

ZIO 






BCC 

A3 

CONNON 

FLJTEXIIOOOO), NANASC3), 

RPRNI2A00), 

SPACR(IO) 

BCC 

AA 









BCC 

A5 



r> o n 





C 

C 


c 


c 

c 

10 

15 

20 

30 

♦ 0 

50 

60 


Jt:f/EXT FLOW FIELD Cm»«ON SECTION 


connoN OELJ » dflf * ijet » iext » inaxj > 

1 KNAX4 t KrAXf , NRJEf > NRtXT , PS JE T » 

2 * P? » P3 » PA r USrj * 

3 VJET > VFXT » XJET # XEXT 


iNAXE 

PSEXT 

USTE 


CUHNON ENOCN 


END OF COHMON fIJ HF COPIFD TO RESTART FILE 
CONNON FID(3»5)f fFNAMO}, LSW0NI]6I 


SIGN « 1,0 

L510 « IK/5K5 .£Q, K 
OU 100 J > 1» NJ 
IF CJ ,GT. 2) SIGN . -1.0 
IROW « IRRRI I> J| 

IF CJ.E0.2 .AND. L510I IROW • 25 
00 80 lEO > NFO 

GO TO I10> 20> 30» A0> 501* lEO 

IF IL5101 GO TO 15 

IF !ALPHAU1»L».E0. O.OI GO TO 80 

TEHP * SIGNATEPIRP/ALPHAUlvLI 

GO TO 60 

TEHP • SIGNACRR 

GO TO 60 

IF IL510I GO TO 80 
IF (ALPHAU2,LI.E0. 0.0} GO TO 80 
TEHP « SlGN*TEPRP/ALPHAfJ2>L I 
GO TO 60 
CONTINUE 

IF ULPHA(J3*L l.EO.0.0) GO TO 80 
TEMP » -SIGN*TEMRM/ALPHA( J3»LI 
GO T(1 60 
continue 

IF lALPHAUAfL } .EO.O.OI GO TO 80 
TEMP • -SIGN*TEMRM/ALPHAUA»LI 
GO TO 60 
CONTINUE 

IF fALPHAI J5»L I.EO. 0.01 GO TO 80 
TEMP - -SIGM*TFM»M/ALPMAI45»LI 
IJ « J123A5(TLQ} 

CM(IPOW»IJ} » CNiropw^IJI ♦ TEMP 


BCC 

BCC 

*BCC 

»BCC 

>BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

ecc 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

see 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 

BCC 


46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 
H4 

85 

86 

87 

88 

89 

90 


f 


f 


60 

CONTINUE 


BCC 

91 


IF IR.EO. 101 GO TO 

100 

BCC 

92 


IF IK .EQ. 5) GU TO 

65 

BCC 

93 


QXflRnWl * OX(IRUW) 

♦ OSIGN*SIGN»l TENRP ♦ RMFAC*TEHRNI 

BCC 

94 


60 10 90 


BCC 

95 

85 

OXfIROkI > OXflROWl 

♦ OS IGN*S IGN*TENRH 

BCC 

96 

90 

CONT INUE 


BCC 

97 

100 

CONTINUE 


BCC 

98 


RETURN 


BCC 

99 


END 


BCC 

100 


c 

c 

c 

c 


c 

c 


c 


SUBROUTINE BOA TCP 

eOATCP - BOAT CHECK POTNT/RE ST ART SUBROUTINE 


BCP 

BCP 

BCP 

ecp 


1 

2 
3 
A 
9 
6 

7 

8 
9 
1 
2 
3 
A 

5 

6 
7 


1 

2 

3 

A 


BCATCP 

- BOAT CONNON 





BCP 

9 

DIHENSION At90)» 

AI0(29)f ALnC(90f6l 

f ALPHAI29 

f90)f CGHV 

(790f 3lf 

BCP 

(> 

CHt25# 

29lf CPBAK(90lf CPTBV(790lf ECC(90i» G(29lf GTBVI 

790lf 

BCP 

7 

HCPI2» 

29lf HSTATI90lf HTBV(790I» IRR(A01> IRRRIA0f9lf 

I 

RTtAOIf 

BCP 

B 

ISAVE(6I* J123A9I9I* PSEXT(90I» PSI(90I» OX(29l> RALPHA 

(29f 90lf 

BCP 

9 

RCIAO* 

3lf PH0(90lf RHn0UT(90)f RT(90lf RU(90I 

f RXE(90ff 

RXK (901 

f BCP 

10 

SlCnA(90l> STARTIll* T(90)» TITLEIIBI* UI90I* 

VEXr (Af 29 

f 2 If 

BCP 

11 

VJETf A»29*2 if W00T(29t90i. WM(29if 

WP(29lf WTNIX(90if 

WTN(lLE(29lfBCP 

12 

XEI90I 

f XEXT(90I 

f XJFT(90)f XK(90lf 

XLE(90lf 

XNU(901f 

Y 190 If 

BCP 

13 

YCUTISOIf 210(91 






BCP 

lA 








BCP 

19 

LOGICAL 

LHALEf LSWON 





BCP 

16 








BCP 

17 

EQUIVALENCE (Jlf 

J123A9( 11 1 





BCP 

18 

equivalence (ALocdfiif cndfiii 





BCP 

19 

EQUIVALENCE (ECCf 

CN(lfl311f (YOUT* 

CNdfl9ll 




BCP 

20 

EQUIVALENCE (CPTBVf CGHV 1 





BCP 

21 

EQUIVALENCE (HCPdfllt WPd II 





BCP 

2? 

EQUIVALENCE (STARTf TTBCltl 





BCP 

23 








BCP 

2A 

CONNON 

TTBOOIf 

HF(29I f CPT0V f 

GTBV f 

HTBV 



BCP 

29 

CONNON 

A f 

AID f ALPHA f 

CARB f 

CH 

P 

CPBAR 

f BCP 

26 


CRR f 

CVISC f OELPSI f 

OFOL f 

DPOX 

9 

DX 

fBCP 

27 


OXNIN f 

FOL f FFF f 

G f 

GGG 

9 

HSTAT 

fBCP 

26 


lOELP f 

IFCC f IFINIS f 

I OUT f 

lOUTl 

9 

I0UT2 

fBCP 

29 


IPAGE f 

IP»ESS f IRR f 

IRRR f 

IRT 

9 

ISAVE 

fBCP 

30 


ITELG f 

IV1S f J1 f 

22 f 

J3 

9 

JA 

fBCP 

31 


Q9 f 

IHAIF f NNfJO f 

NPSI f 

NXNPT 

9 

NXNPl 

fBCP 

32 


NPSI f 

NR f NRAO f 

NRAS f 

NS 

9 

NT 

fBCP 

33 


P f 

PCNT f PRNT * 

PRNTXC f 

PS I 

9 

PSIO 

fBCP 

3A 


QX f 

RALPHA f RBUrjY f 

RC f 

RHO 

9 

RHOUUT 

fBCP 

39 


RJ f 

RT f RTACU f 

RTJAC f 

RTJOB 

9 

RTHAX 

fBCP 

36 


RU f 

PXE f RXK f 

SIGE f 

SICK 

9 

SIGNA 

fBCP 

37 


T f 

TCONT f TEOGE f 

TENRN f 

TENRP 

9 

TITLE 

fBCP 

30 


TKINET f 

U f UNIT f 

WOOT f 

WN 

9 

HP 

fBCP 

39 


WTHIX f 

WTNOir • X , 

XCHANG f 

XO 

9 

XE 

fBCP 

AO 


XINIT f 

XK f XK2 f 

XLE f 

XNAX 

9 

XHU 

fBCP 

Al 


r f 

210 





BCP 

A2 

CONNON 

FLJTEXdOOOOlf NANASIlif RPRN(2A00)f 

SPACRdOl 

BCP 

AJ 








BCP 

AA 

JET/EXT 

FLOW FIELD CONNON SECTION 





BCP 

A9 


C 

C 


c 


BCf» 





COHHON 

oeu 

* 

OFLE 

» 

IJET 

r 

lEXT 

* INAXJ 

* IBAXE 

*BCP 

67 


1 

KNAXJ 

* 

KNAXE 


NRJET 

p 

NREXT 

f PSiET 

* PSEXT 

*BCP 

60 


2 

PI 

* 

P2 

$ 

P3 

0 

P6 

* USTJ 

* USTC 

»BCP 

69 


3 

VJET 

* 

vrxT 

» 

XJET 

0 

XEXT 



BCP 

50 

c 











BCP 

51 


CONNON 

ENOCH 









BCP 

52 

c 











BCP 

53 

c 

(NO OF 

CONNON 

TO 

BF CnPlFf) 

TP restart 
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Bl 

F 

62 

c 




Bl 

F 

63 


DATA IFL0W/2HJE»2HFX/ 



61 

IF 

6A 

c 

25 RADIAL L0CATIUNS»50 AXIAL LOCATIONS 



Bl 

IF 

65 

c 




Bl 

[F 

66 

c 

STARTING RECORD NO. OF FLOW FIELDS FOR 

JET/EXT (NWPR«1) 


Bl 

IF 

67 

c 




Bl 

F 

66 


NREC«1 



Bl 

IF 

69 

c 




Bl 

F 

70 

c 

READ INVISCIO DATA ARRAYS 



Bl 

IF 

71 

c 




Bl 

F 

72 


DO 100 IP>1*IPRESS 



Bl 

F 

73 


IMAX«IMAXV( IP) 



Bl 

F 

7A 


NRECVUP)-NREC 



Bl 

F 

75 


WRITE (6>150) IFLOWdP) 



Bl 

F 

76 


DO 90 I«1»IMAX 



Bl 

F 

77 


IF (IP.EO.l) READ I5»l<^0) LnAP*XJET(I) 



Bl 

F 

78 


IF (IP.E0.2) READ (5»1<.0) LNAP*XEXTU) 



Bl 

F 

79 


00 10 J>1*^ 



Bl 

F 

80 


READ (5»160) ( VJET( J»K»1 )»K>1»LNAP) 



Bl 

F 

BI 

10 

CONTINUE 



Bl 

F 

62 

c 

VJEni*K>l)«Y;VJET<2>K»l)«P;VJET(3»K»l)«T;VJETIA»K*l)-U 


Bl 

IF 

63 

c 




BI 

F 

BA 

c 

MASS FLOW INTEGRATITN 



Bl 

F 

65 

c 




Bl 

F 

86 


MOLOIl l«0. 



Bl 

F 

87 


IF IIP.E0.2) HULD(1)«PSJET**? 



Bl 

F 

68 


IF IIP.EQ.I) WDUM«S«PR2INS#ALPHAC1»1I» 

WTNOLE) 


Bl 

F 

89 


IF (IP.E0.2) WDUM=SMPP?(NS»ALPHAI1»MPSI IfWTMOLE ) 


Bl 

F 

90 


"W 





DUKl>VJETf 2*1>1 l*WOun*VJET(4>l»U/VJETO>l» n 

BIF 

91 


00 20 K-2.LHAP 

BIF 

92 


0UI12«VJETf2.H»ll*W0UM*VJEn«,K»n/VJET<3»K»ll 

BIF 

93 


TEKM».S*<0U«1*0UM2I 

BIF 

9^ 


TE»»H«TE«M/A2,2a5 

BIF 

95 


HOLDIK )«H0L0(K-1 l♦TrPH*<VJ^Tf l»K,n**2-VJ£Tll»K-l, 1 >♦♦21 

BIF 

96 


0UI11-0UH2 

BIF 

97 

20 

CONTINUE 

BIF 

90 


00 30 tfl*LHAP 

BIF 

99 

30 

HaLO(K)>SORT(HOLO(K I) 

BIF 

100 


IF UP.FO.l.ANO. I.EO.n PSJFT-HOLOCLNAPI 

BIF 

101 


IF fIP.EQ.2i PSEXTm»HniOtlHAPl 

BIF 

102 

C 


BIF 

103 

C 

NONDINENSIONALIZF NASS FLOW 

BIF 

lOA 

C 


BIF 

105 


H0LDS«H0L0<1 1 

BIF 

106 


DO ^0 K>1#LNAP 

BIF 

107 


HOLOCK) •IMOLOIKI-HniOS >/ IHOLOILHAP )-HULOSI 

BIF 

lOB 

^0 

CONTINUE 

BIF 

109 

C 


BIF 

110 

C 

INTERPOLATE ONTO EVEMLT SPACED ARRAY 

BIF 

111 

C 


BIF 

112 


IF (IP.EO.li KHAX»KMAJ(J 

BIF 

113 


IF (IP.E0.2I KNAX-KMAXF 

BIF 

116 


DO 70 LL-l»A 

BIF 

115 


00 50 K«1»LHAP 

BIF 

116 


HOLOlfK )«VJETILL»K> 1 ) 

BIF 

117 

SO 

CONTINUE 

BIF 

IIB 


DO 60 K«1»KHAX 

BIF 

119 


PSOUN-FLOATU-1 » /FLOAT (KNAX-1 ) 

BIF 

120 


CALL BOATLI (PSOUN » VEX T ( LL *K * 1 ) >H0L0»H0L 01 » LNA P > 

BIF 

121 

60 

CONTINUE 

BIF 

122 

70 

CONTINUE 

BIF 

123 

C 


BIF 

126 


PNFLO-PSJET 

BIF 

125 


IF (IP.EQ.2) PMFLO-PSEXT II » 

BIF 

126 


WRITE I6«170l 1*XJEV(I>IPI»PNFL0 

BIF 

127 


DO 80 K*1»KNAX 

BIF 

126 


PSOUn-FLOATIK-1 I /FLOAT f RNAX-1 ) 

BIF 

129 


WRITE I6»160) K, (VEXT (J»K»1 IfPSDUH 

BIF 

130 

60 

CONTINUE 

BIF 

131 

C 


BIF 

132 

C 

WRITE ON DISC 

BIF 

133 

C 


BIF 

136 


CALL PBFDW I F 1 0> liR F C > AAKNA X > VF X T 1 

BIF 

135 




■■■-- ;iiPlff‘wii!i(l^’iii!'i :ir".ii':; lillPrpWip 


C 

90 

100 

C 

C 

C 


110 


C 

C 

c 


120 

130 

C 

1<|0 

150 

160 

170 



BIF 

136 

CONTINUE 

BIF 

137 

CONTINUE 

BIF 

138 


BIF 

139 

READ INITIAL TWO STATIONS OF JET MAP 

BIF 

160 


BIF 

161 

DO 110 I*l>2 

BIF 

162 

CALL PBFOR ( F I D» NR JE T » <**KNAX J » V JE T 11 » 1 f 1 1> 

BIF 

163 

CONTINUE 

BIF 

164 

IJET-2 

BIF 

165 

IF CIPRESS.EO.n GO TO 130 

BIF 

166 


BIF 

167 

READ INITIAL TWO STATIONS OF EXT NAP 

BIF 

168 


BIF 

169 

DO 120 1>1#2 

BIF 

150 

CALL PBFOR ( F 1 D»NRF XT, <i*KM AXE » VE XT 1 1 » 1» I» ) 

BIF 

151 

CONTINUE 

BIF 

152 

IEXT«2 

BIF 

153 

CONTINUE 

BIF 

156 

RETURN 

BIF 

155 


BIF 

156 

FORMAT (I2,3X,7E10.3) 

BIF 

157 

FURNAT IlHl,^0X,A2,l<iHT,INVISCID MAPI 

BIF 

158 

FORMAT feE10.3l 

BIF 

159 

FORMAT (///15H STATION NUMBER , 1 <i , 1 OX , 1 6HAX I AL 

LOCATION •,EI3.5*10XBIF 

160 

1,21HMASS FLOW PARAMETER « » E 1 3 . 5 //^ X . 2HP T , l<t X, 

1HY,12X,1HP,12X,1HT,1BIF 

161 

22X,1HU,7X,6HPSIPaRI 

BIF 

162 

HJRMAT (IH , I5»2X» 8F13.5* 

BIF 

163 

END 

BIF 

166 


180 


c 

c 

c 


c 

c 


10 


?0 

c 


c 

30 


GO Tn 10 


13» V^^PsirFlO^NPECI 

INVISCIO iNTERPOLATrON POUTIUE 

ir^'' — •* w... 

CALL PBFOP <FIO,NPfC.^*KMAXV,VECTCI,i,an 
IV-IV*1 

1^ < IV. LE. inAXVI 

write <6, 30) 

STOP 

CONTINUE 

IC ! ?!^’^*^'**1*0R.ITYP,E0.1) PSVV*PSVM 1 
OSI-l./FlOATCKNAXV-l) 

IF SI8AR.PSIX/PSVV 

KV-SIBAR/DSUi****’*^^'***^^*^ ’ SI9AR-CPSIX-PSJ l/TPSVV-PSJ) 
IP IKV.LT.l) KV«1 

''“51 

VV<N)-OOH1»RATX*<OUP?-Ouin 

CONtInUe’**^'^* WV(N )«(OUHP-OUNl )/(XVflV)-XV(lV-l ) ) 

Vl«r*V2«p,V3*T»VA«U 
Vl«vvu I 

V2*VVC2)*?ii7, 

V3»VV(3) 

VA-VVIA ) 

return 


fOM.t l,3H BO.ril - „ ,r.T. ,„„V) 


V2»VBII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BI T 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
-DBII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 
BII 


1 

2 

3 

A 

5 

6 
7 
B 
9 

10 

11 

12 

13 

lA 

15 

16 
17 
IB 

19 

20 
21 
22 
23 
2A 

25 

26 
27 
20 

29 

30 

31 

32 

33 
3A 

35 

36 

37 
36 
39 
AO 
A1 
A2 




SUBROUTINE BOATtN 
INITIAL DATA OUTPUT 
0OATCH - BOAT CUhPON 

DIMENSION AI50l» AI0I75I# AL0C(50*6I» ALPHA (25*90}* CGHV(750*3)* 

1 CM(25*25I* CPBAR(50»* CPTBV(750»* ECC(50)> GI25}* GTBV(750i* 

2 MCP<2*25}* HSTAT(50I* HTBV(750I. IRR(<^OI* IRPR(50*9I» IRTI50I* 

3 ISAVEIGI* JI23A5I5I* PSEXT(5P»* PSK50I* OX(251» RALPHAf 2S* 501 * 

A KC(A0*3)» RHIJI501* RH0PUT(90)» RT(501* RUC50I* RXE(501* RXK(50I> 

5 SIGMA(501* STARTCll, TC501* TITLEdBl, U(501* VEX T I A* 25* 2 1 » 

6 VJET(A*25*2I* MDUT(25*501* WM(25)» WPI25)* WTMIX(50)* WTMULE(25I 

7 XE(50)» XEXT(50}* XJFT(501* XA(50I* XLE(50>* XMUI50I* Yf50l* 

8 YUUT(501» ZID(51 

LOGICAL LHALF*LSWON 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


COMMON 

COMMON 


Ul* J123A5I1II 

(AL0C(1»1I* CM(1*1II 

(ECC* CM(l*13n* (YOUT* CM(1*15I1 

(CPTBV* CGHVl 

(HCP(1*1I* UMfin 

(START* TTBTin 


BIN 

BIN 

BIN 

BIN 

BIN 

SIN 

BIN 

BIN 

BIN 

BIN 

BIN 

*BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 

BIN 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

lA 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 


COMMON 


TTBT30I 

MF(25) 

CPTBV 

f 

GTBV 

P 

HTBV 



BIN 

25 

A 

AID 

ALPHA 

# 

CARB 

P 

CM 

P 

CPBAR 

*BIN 

26 

CRR 

CVISC 

0ELP51 

P 

OFOL 

P 

DPOX 

P 

OX 

»BIN 

27 

DXMIH 

FOL 

FFF 

P 

G 

P 

GG6 

P 

HSTAT 

*BIN 

28 

lOELP 

IFCC 

ifinis 

P 

lOUT 

P 

lOUTl 

P 

IUUT2 

»BIN 

29 

IPAGE 

IPRESS 

IRR 

P 

IRRR 

P 

IRT 

P 

ISAVE 

*BIN 

30 

ITFLG 

IVIS 

J1 

P 

J2 

P 

J3 

P 

JA 

>BIN 

31 

J5 

LHALF 

MMOO 

P 

MPSI 

P 

MXNPT 

P 

MXNPl 

*B1N 

32 

NPSI 

NR 

NRAO 

P 

NRAS 

P 

NS 

P 

NT 

*BIN 

33 

P 

PCNT 

PRNT 

P 

PRNIXC 

P 

PS I 

P 

PSIO 

*BIN 

3A 

OX 

PALPHA 

R6UQY 

P 

RC 

P 

RHO 

P 

RHOOUT 

*BIN 

35 

RJ 

RT 

RTACU 

P 

RTJAC 

P 

RTJOB 

P 

RTMAX 

f BIN 

36 

RU 

RXF 

RXK 

P 

SIGE 

P 

SICK 

P 

SIGMA 

»BIN 

37 

T 

TCONT 

TFDGE 

P 

TEMRM 

P 

TEMRP 

P 

TITLE 

*BIN 

38 

TKIMET 

U 

UNIT 

P 

WDOT 

P 

WM 

P 

WP 

*BIN 

39 

WTMIX 

WTMOLE 

X 

P 

XCHANG 

P 

XD 

P 

XL 

»BIN 

AO 

XINIT 

XR 

XR2 

P 

XLE 

P 

XMAX 

P 

XMU 

*BIN 

A1 

Y 

710 








BIN 

A2 

FLJTEX (10000)* NaMASOT* 

RPRMI2A00I 

P 

SPACR( 101 


BIN 

A3 


BIN 

BIN 


AA 

A5 


JET/EXT FLOW FIELD COMMnN SECTION 


Tipn' I 


c 





BIN 

66 


cunnuN OEU » dele » ijet » lEXf 

* 

INAXJ 

* ihaxe 

*B1N 

67 


1 KHAXJ , K>«AXE > N9JET » NREXT 

* 

PSJET 

* PSEXT 

*B1N 

6B 


2 PI » P? » P3 » PA 

* 

USTJ 

* USTE 

*6IN 

69 


3 VJET » VFXT » XJET p XEXf 




BIN 

50 

c 





BIN 

51 


COH'IOM ENOCH 




BIN 

52 

c 





BIN 

53 

c 

END OF COHHON TO BF COPIED TO RESTART FILE 




BIN 

56 

c 





BIN 

55 


CONNON FIDf3»5}» IFNAri3)* LSU0Ntl6l 




BIN 

56 

c 





BIN 

57 

c 





BIN 

5B 


CALL DATE (ROATE) 




BIN 

59 


WRITE I6t300) 




BIN 

60 


WRITE t6»310) 




BIN 

61 


WRITE t6*320l ( T I TL E <1 ) > 1 > 1 » 1 R I 




BIN 

62 


WRITE (6»I50I IFriAN»NANAS»ROATE 




BIN 

63 


IF (IPRESS.NE.OI WRITE (6*3301 P 




BIN 

66 


IF (IPRESS.EO.O) WRITE (6*3A0I P 




BIN 

65 


WRITE (6*3701 RJ*RBUnY 




BIN 

66 


WRITE (6*360) XLE( 1 1* SIGHA( 1 ) 




BIN 

67 


WRITE (6»A10) X*XNAX 




BIN 

6B 


WRITE (6*A20) PRNT*OXNIN 




BIN 

69 


IF (IVIS.LT.OI WRITE (6*160) 




BIN 

70 


IF (IVIS.EO.O) WRITE (6*170) 




BIN 

71 


IF (IVIS.EO.l) WRITE (6*160) 




BIN 

72 


WRITE (6*390) 




BIN 

73 


WRITE (6*360) T(1)*T(NPSI) 




BIN 

76 


WRITE (6*360) U(1)*U(NPSI) 




BIN 

75 


WTNIX( D-O.O 




BIN 

76 


WTniX(NPSI )«0.0 




BIN 

77 


DO 10 J-1*NS 




BIN 

7B 


wrntxd )«wTntx(i)f ALPHA (j*i) 




BIN 

79 

10 

WTN1X(NPSI )>WTNIX(PPSI ) ♦ AL PHA ( J* HP SI ) 




BIN 

80 


DO 20 J>l*NS 




BIN 

B1 


RALPHA( J*1)*ALPHA( J*1 )/WTNIX(l) 




BIN 

82 


RALPHAI J*NPSI )«ALPHA(J*MPSI)/WTNIX(nPSl ) 




BIN 

B3 

20 

WRITE (6*600) AI0(J)*RALPHA(J*1)»RAIPHA(J*NPS1) 



BIN 

B6 


IF (NR.LE.O) GO TO 160 




BIN 

B5 


WRITE (6*190) 




BIN 

B6 


DO 130 I>1*NR 




BIN 

87 


L«1RR( I ) 




BIN 

8B 


GO TO (30*60*50*60*70*60*90*100*110*120)* L 




BIN 

B9 

30 

J1-IRRR(1*1 ) 




BIN 

90 


^0 


w 


60 


70 


80 


90 


100 




BIN 

91 

J3-IRRRI I. 31 


BIN 

92 

J6-IRRR (I»4) 


BIN 

93 

kRITE f6>?00l 

I*AID(J1 )>A(0(J2I*AI0(J3)*A10(J6)* (RC( 1* J)* J*l*3) 

BIN 

96 

GO TO 130 


BIN 

99 

J1«IRRP( (* 11 


BIN 

96 

J2-IRRR n»2) 


BIN 

97 

J3>IRRR<I»3I 


BIN 

9B 

WRITE f6»210l 

If AI0(J1 )»AI0( J2)»AI0(J3)*(RC (I* J)*J*1* 3) 

BIN 

99 

CO TO 130 


BIN 

100 

J1*IRRR(I*1 1 


BIN 

101 

J2*IRRR(I>Z} 


BIN 

102 

J3*IRRR(I>3 1 


BIN 

103 

J6-|RRRU»6I 


BIN 

106 

J9«IRRRII»3) 


BIN 

109 

WRITE (6f?20l 

IfAKHJl )»AIO(J?)*AID( J3)* AID (J6)*A10(J9 )* (RC(1*J 1* JOIN 

106 

1*1*3} 


BIN 

107 

GO TO 130 


BIN 

lOB 

J1*IRRRU*11 


BIN 

109 

j;?-IRRRf 1*21 


BIN 

110 

J3«[RRRII*3} 


BIN 

111 

WRITE (6*2301 

I*AI0(J1 )*AI0(J2)*AI0(J3)* (RC(I* J)*J«1*3) 

BIN 

112 

GO TG 130 


BIN 

113 

Jl-IRRRI I* 1) 


BIN 

116 

J?-IRRRII»3) 


BIN 

119 

J3*IRRR(1*6I 


BIN 

116 

WRITE (6*2601 

I*AI0(J1)*AI0(J2)*AI0(J3)*(RC(I*J)*J«1*3) 

BIN 

117 

GO TO 130 


BIN 

IIB 

J1>IRRR(I»1) 


BIN 

119 

J2>1RRR(I*2) 


BIN 

120 

J3*IRRR(1*3) 


BIN 

121 

J6«IRRR(I»6) 


BIN 

12? 

WRITE (6*290) 

I*AI0(J1 }*AI0( J2)* AI0( J3)*AI0(J6 )*(RC( I* J)* J»l*3) 

BIN 

123 

GO TO 130 


BIN 

126 

J1«IRRR( 1*11 


BIN 

129 

J2-IRRR(I*2I 


BIN 

126 

J3*IRRR(I*3) 


BIN 

127 

WRITE (6*260) 

I*AI0(J1 )»AI0( J?l*AI0( J3)* (RC(I» J)*J«1» 3) 

BIN 

1?B 

GO TO 130 


BIN 

129 

Jl-lRRRI 1*1) 


BIN 

130 

J2-IRRR( 1*2) 


BIN 

131 

J3*IRRR( 1*3) 


BIN 

132 

J6-IRRR(1*6) 


BIN 

133 

J9«IRRR( 1*9) 


BIN 

136 

WRITE (6*270) 

l*Ain(Jl)*AID(J2l*AIO(J3)*AlD(J6)*AIO(J9)*(RC(I*J)* 

JBIN 

135 


110 


120 


130 

1^0 

C 

150 

160 

179 

leo 

190 

200 

210 

220 

230 

290 

250 

260 

270 

280 

290 

300 

310 

320 

330 


1>1»3I 
GO TO 130 
Jl«IRRRf l»ll 
J2-IRRRC 1*21 
J3-IRRRtI»3} 
URITE (6>2e01 
GO TO 130 
J1*IRRRI I» II 
J2-IRRRU.2I 
J3>IRRRf 1» 31 
MRITE <6»299) 
CONTINUE 
CONTINUE 
RETURN 


l»A10UU»*ID(J2l»AIDf431*IRCII»J)>J«l*3) 


IfAIDUl l»AIO( J2ltAI0I J3I> IRCII* JI»J«1» 31 


UNITS 1 1 


FORMAT llH0»22X»93MtNITIALlZATinN OF 8UAT RUN, RESTART FILE - 
1,21H, BOAT^SPECRA FILE - *3A?,SX»A101 
FORNAT tlH0»22X#20HKE? TURBULENCE NODE L I 
FORMAT I1H0,22X#27HPRANDTL "IXING LENGTH MODELI 
FORMAT f lHO»22Xf 20H0nNAL0S0N/GRAY MODELI 
FORNAT f 1H0,19X,26HPFACTI0NS BEING CONS lOEREO, 6X, 1 5H 
1/,9HT**N, 7Xt 1HA,BX, 1HN,9X, 1HB»7X,23HIM0LECULE-NL-SEC 
FORMAT II9»9X,A<>,2H« ,A6»RX*2Ha #A6»2H» ,A6#18X#1E10 
IFlO.ll 

FORNAT f I9f9X»A6,2H« , A6» 3H«> M,5X#2H« fA6,3H« M,23X, 
U2X,F10.1I 

FORNAT n9,9X,A6»2H«^ tAGfeRiZH* ,A6»2H« »A6,2H* ,A6» 
1.1,2X»FI0.II 
FORMAT I 19»9X,A6»2H« 

FORMAT f I9>9Xf A6>3H» 

12X,F10, II 

FORMAT I I9»9Xt A6*2H» ,A6»8X#2H- 
10. If 3X, 16HUHE WAT REACTION I 

FORMAT f l9f9Xt A6f 2H« f A6f 3H« Mf5Xf2H> fA6f3H€ Nf 23X# E9.3f 2Xf 
IXf FlO.lf 3Xf 16HUNE HAY REACTtONI 
FORMAT II9f9Xf A6f2H* fA6fBXf2H« fA6f?H9 fA6f2H« fA6flOXf£9.3 
l.lf2XfF10.1f 3Xf IGHfiME WAY » •CTION) 

FORMAT < l9f9Xf Abf 2M« fA6fBX.2H« f A6f 26Xf E9. 3f 2Xf F9 . I f 2Xf FIO • 
16H0NE HAY REACTIIJNI 

FORMAT ( I9f9Xf A6f 3H« M»13Xf2H> fA6f2H* fA6f3H« M f 1 5Xf E 9. 3t 2X 
12XfF10. If 3Xf 16H0NE WAY REACTION} 

FORMAT lIHlf 37Xf95HAEROFAUTICAL RESEARCH ASSOCIATES OF PRINC 
FORMAT l35Xf 50HAXISYMMETRIC MIXING WITH NON-F OUI L I BR 1 UM CHEN 
FORMAT llH0f29Xf 1BA9) 

FORMAT tlH0f?2Xf 19HrRESSURE f INITIAL } •E15.7fl2H ATMOSPHERES} 


f A6fBXf2H« f A6f26XfE9.3f 2XfF9.1f 2Xf FlO.l 
Mfl3Xf2He fA6f2H« fA6f3H« Mf 15Xf E9. 3f 2Xf 


f A6f2H« 


BIN 

136 

BIN 

137 

BIN 

138 

BIN 

139 

BIN 

190 

BIN 

191 

BIN 

192 

BIN 

193 

BIN 

199 

BIN 

195 

BIN 

196 

BIN 

197 

BIN 

198 

BIN 

199 

BIN 

150 

• f3A2BIN 

151 

BIN 
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9 

RXK 

9 

SIGE 

9 

SIGK 

9 

SIGHA 

>BLP 

35 

3 

T 

TcnNT 

9 

TEOGE 

9 

TEHRH 

9 

TEHRP 

9 

TITLE 

f BLP 

35 

h 

TKINE X 

U 

9 

UNIT 

9 

WDOT 

9 

wn 

9 

WP 

»BLP 

36 

5 

UTHIX 

WTHOLE 

9 

X 

9 

XCHANG 

9 

xo 

9 

Xt 

»8LP 

37 

6 

XINI r 

XK 

9 

XK2 

9 

XLE 

9 

XNAX 

9 

XNU 

»BLP 

38 

7 

Y 

ZIO 









BLP 

39 

COHHUM 

FLJTEXI 10000) * 

NAHASI3)* 

RPRH(2A00) 

9 

SPACR (10) 


BLP 

50 












BLP 

51 

JET/EXT 

FLOW FIELD CQHHON SECTION 







BLP 

52 












B^P 

53 

COHHDN 

OELJ j 

DELE 

9 

I JET 

9 

lEXT 

9 

IHAXJ 

9 

IHAXE 

»BLP 

55 

1 

KHAXJ 

KHAXF 

9 

NPJET 

9 

NREXT 

9 

PSJET 

9 

PSEXT 

• BLP 

55 




r» r» o 




-“»‘”i!iiilf'!l«mili:fl!fllp 



PI 

VJET 

P 

P 

P2 » 

VEYT » 

P3 » 

XJEI » 

XEXT 

> USTJ 

CONNOM 

FMOCN 






END OF 

CONNON 

TO 

BE CPPIEO 

TO RESTART 

FILE 


CONNON 

FID13» 

51 

, IFNAr(3l» LSW0NC16I 



c 

C 

OIHENSIOH AT8I2I 

NREC • tl-ll*30 ♦ ITKFY 

CALL SFVt'V (CGHVINREC» IFILE-1 »» ATB» 2l 

AX» AI8m*H0l^ ATB»2)*S0T 

RETURN 

END 


USTE »BLR ^6 
BLP ^7 
BLP 

BLP ^9 
BLP 50 
BLP 51 
BLP 52 
BLP 53 
BLP 5^ 
BLP 55 
BLP 56 
BLP 57 
BLP 5B 
BLP 59 
BLP 60 
BLP 61 






noon 




C 

C 


c 


c 

c 


SUBROUTINE nUATLU (INTTI 
BQATLU - BOAT FILf LOOK UP 
80ATCH - BOAT CUMNON 

DINENSIQN A(50lf AID(?5)* ALnCC50»6}> ALPHA I 2 $0 1 * CGHV(730f3>* 

1 CHCP5*25I» CPBARJ5CI, CPTBV(75a), ECC(30I* G»25l» GTBVI730)* 

2 HCP(2»23I» HSTAT(50»* HTBV(730I* IRR(AO)* IKRRIA0»5)» IRT(<iO»» 

3 1SAVE(6)» J12333(3>> PSE*T(50I» PSK30I» OXt23>* R AL PHA { 23* 30 I * 

A RC(<iO»3}» RHQ(30>* »Hn(’UT ( 53 » * R1(30>* RUC50»* RXk(30)» RXKI30I* 

5 SlGnA(30)> STARHlIt T«50l» TlTLEtlB)* U(50)» VEX T I * 23* 2 ) > 

6 VJET(A»23»2)» WOO T ( 2 3 * 3 0 ) » Wff25l* WP(23)» WTNIX(30I* WTf10LE(23l 

7 XH50I* XEXT(30>* XjETISOI* XKf50l* XLE(30)* XNU(30I* Y(50l* 

8 rOUTISOI* ZID(3) 

LOGICAL LHALF*LSWQN 


EQUIVALENCE 

equivalence 

EQUIVALENCE 

equivalence 

EQUIVALENCE 

EQUIVALENCE 


IJl* J123<i3(lll 
CALOCa*!)* CN(1*]|| 

(ECC* CHn*13Il* (TOUT* CN(1*13)> 
CCPTUV* CGHVI 
(HCPU*1>* WN(in 
(START, TTBUn 


BLU 

BLU 

BLU 

8LU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

*BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 


1 

2 

3 

V 

3 

6 

7 

6 

9 

10 

11 

12 

13 

13 

13 

16 

17 

10 

19 

20 
21 
22 
23 
23 


CONNON 

TT8(30» 

MF(25I , 

CPT8V 

9 

GTBV 

9 

HTBV 



BLU 

25 

CONHON 

A 

AID * 

ALPHA 

9 

CARB 

9 

cn 

9 

CPBAR 

*BLU 

26 

1 

CRR 

CVISC , 

OELPSI 

9 

OFDL 

9 

OPOX 

9 

DX 

*BLU 

27 

2 

oxniN 

FOL , 

FFF 

9 

G 

9 

GGG 

9 

HSTAT 

*BLU 

28 

3 

lOELP 

IFCC * 

IFINIS 

9 

lOUT 

9 

lOUTl 

9 

I0UT2 

*BLU 

29 

3 

IPAGE 

IPRESS , 

IPR 

9 

IRKR 

9 

IRf 

9 

ISAVE 

*BLU 

30 

3 

ITFLG 

IVIS * 

J1 

9 

J2 

9 

J3 

9 

J3 

*BLU 

31 

6 

J3 

LHALF • 

HNOO 

9 

flPSl 

9 

nXNPT 

9 

nxNPi 

*BLU 

32 

7 

NPSI 

NR * 

NRAO 

9 

NRAS 

9 

NS 

9 

NT 

*BLU 

33 

8 

P 

PCMT , 

PPNT 

9 

PRNTXC 

9 

PS I 

9 

PSIO 

*BLU 

33 

9 

OX 

PALPHA * 

RBUOr 

9 

RC 

9 

RHO 

9 

RHOOUT 

*BLU 

35 

1 

RJ 

PT * 

RTACU 

9 

RTJAC 

9 

RTJIIB 

9 

RTNAX 

*BLU 

36 

2 

RU 

Rxr * 

RXK 

9 

SIGE 

9 

SIGK 

9 

SIGNA 

*BLU 

37 

3 

T 

TCONT * 

TFOGE 

9 

TENRN 

9 

TENRP 

9 

title 

*BLU 

39 

3 

TKINET 

U , 

UNIT 

9 

WDOT 

9 

WN 

9 

WP 

*BLU 

39 

5 

WTNIX 

WTHOLC , 

X 

9 

XCHANG 

9 

XO 

9 

XE 

*BLU 

30 

6 

XINIT 

XK t 

XK2 

9 

XLE 

9 

XHAX 

9 

xnu 

»BLU 

31 

7 

y 

710 








BLU 

32 

connoN 

FL JTEX(IOOOO)* NAPAS(3I* 

RPRNI2300) 

9 

SPACRI 101 


BLU 

33 


JET/FXT FLOW FILLO CONMON SECTION 


BLU 

BLU 


33 

35 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

10 

c 

c 

c 

c 

c 

c 


20 

30 


c 

‘tO 

c 

iO 

60 


CONNON 

DELJ 


OELF 

9 

IJET 

9 

lEXT 

# INAXJ # 

INAXE 

1 

KNAX J 

P 

knaxe 

9 

NRJE T 

9 

NREXT 

# PSJET # 

PSEXT 

2 

PI 

f 

P2 

9 

P3 

9 

PA 

# USTJ # 

USTE 

3 

VJET 

9 

VEXT 

9 

XJE T 

9 

XEXT 



CONNON 

ENDCN 









END OF 

CUNNON 

TO 

BE COPIED 

TfJ RESTART 

FILE 




COnnUN F10n*5l» IFNAMI), LSW0NI16I 


PEWINO PESTA9T FILE 
REWIND 1 

IF UNIT. GT. 01 GU TO 10 

ZERO RPRN FILE flN STAPT-UP 

CALL SFVFL I0.0»RPRM»30*A0*? I 
CONTINUE 

LOOK UP 0OAT-SPLCPA OUTPUT FILE 
IF CNRAO.EO.O) GU TO AO 

BOAT>SPECRA FILE# FOPNATTEO SEOUENTIAL# POSITION AT END-OF-DATA 

IF UNIT. EO. 01 END FILE 2 

REWIND 2 

READ 12*901 

IF lEOFIzn 30#20 

BACKSPACE 2 

GU TO AO 

WRITE (6>60l NANAS 
STOP 

CONTINUE 

RETURN 

FORNAT (A0A2I 

FORNAT (25H 3UAT-SPCCRA FRF ERROR I 
END 


BLU 

»BLU 

#BLU 

»BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 

BLU 
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BLU 

BLU 
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BLU 

BLU 

BLU 
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BLU 
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A7 

Ad 

A9 

90 

91 

92 
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99 
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97 

98 

59 

60 
61 
62 
63 
6A 


6B 

69 

70 

71 

72 
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7A 
79 

76 

77 

78 

79 

80 
81 
62 
63 
BA 

89 
86 
67 
88 
B9 
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1 


f 

I 

I 


C 

C 


c 


c 

c 


SUBROUTINE BOATNl 





BNl 


NAIN INTEGRATION 

CONTROL ROUTINE 




BNl 

BNl 

1 

3 






811 

A 

BCATCR - BOAT CUHRON 




BNl 

5 

DINENSIQN ktbOi, 

AID(?5I, ALnCl50,6l 

* ALPHA(^5 

•50}. CGHVI750.3}* 

BNl 

6 

1 CNI?5»?5I» CPRARI30I> CRTBVI750lf ECCI50lt CCZ5I> GTBVI750)* 

BNl 

7 

1 HSTATI50I, HTBVC750I» IRRIAOI* IKRR(A0»$I* 

IRT( AOI* 

BNl 

8 

J ISAVEI6I* J1?3AS(^I» RSEXTISOI* PS1I50)* OXI25lt RALPHAf 

BNl 

9 

% RCtAO*)!* RHOtSO)* RHOOUTCSOI* RTI50I» RUI50) 

f RXEI50) 

» RXKr50l 

> BNl 

10 

5 SIGNAI30)* STARMlIt H 50 1 » TIUEU0I* U(50l» 

VEXT (5.25*2}> 

BNl 

11 

!> VJET(AfZ5»2 >» W0UT(25*50)> WN(?5>* 

WP(25)» WTNIX(50I« 

WTH0LE(?5I*BN1 

12 

7 XC(50»» XEXT(50) 

. XJFKSOI* XKI50)* 

XLE (50I» 

XNU(50lf 

TI50I» 

BNl 

13 

9 rOUTISOI# ZI0*5» 





BNl 

lA 






BNl 

15 

LOGICAL LHALFvLSWON 




BNl 

16 






BNl 

17 

EQUIVALENCE tJl» 

J1?3A5( 11 I 




BNl 

18 

equivalence lALOCUfll* c*ici,in 




BNl 

19 

EQUIVALENCE ltCC» 

CNIl»13n» (YOUTf 

CN(l»i51) 



BNl 

20 

EQUIVALENCE fCRTSV* CGHVI 




BNl 

21 

EQUIVALENCE lHCPfl>ll# WNIIM 




BNl 

22 

equivalence (start* TlHdll 




BNl 

23 






BNl 

2A 

CUNNON TTBOOI* 

HF(Z5I » CPT8V * 

GTBV * 

HTBV 


BNl 

25 

CUNNON A » 

AID * ALPHA * 

CARS # 

CN 

CPBAR 

#BN1 

26 

CRR > 

CVISC * OELPSI » 

DFOL t 

OPDX 

DX 

»BN1 

27 

oxniN » 

FOL > FFF • 

G * 

GGG 

HSTAT 

fBNl 

2B 

lOELP » 

lECC » IFINIS , 

I OUT » 

lOUTl 

I0UT2 

• BNl 

29 

IPAGE * 

IPRESS » IRR » 

IRRR * 

IRT 

ISAVE 

• BNl 

30 

itflg * 

IVIS » Jl , 

Q2 * 

J3 

JA 

• BNl 

31 

JS t 

LHALF , MNOO » 

NPSI > 

NXNPT 

NXNPl 

• BNl 

32 

NPSI » 

NP * NRAO » 

NRAS • 

NS 

NT 

• BNl 

33 

P » 

PCNT » PPNT f 

PRNTXC * 

PS I 

PSIO 

• BNl 

3A 

OX * 

RALPHA * PB'JOT t 

RC > 

RHO 

RHOOUT 

• BNl 

35 

R J f 

RT » RTACU > 

rtjac • 

RUOB 

RTNAX 

• BNl 

36 

PU > 

PXF > RXK • 

SIGE * 

SIGK 

SIGNA 

• BNl 

37 

T * 

ICUNT » TE DGE * 

TENRN * 

TENRP 

TITLE 

• BNl 

30 

trinet , 

U » UNIT * 

UOUT > 

wn 

WP 

• BNl 

39 

WTfllX f 

WTNOLF * X , 

XCHANG * 

xo 

XE 

• BNl 

AO 

XINIT * 

XR > XK2 , 

XLE • 

XNAX 

XHU 

• BNl 

A1 

r » 

no 




BNl 

A2 

CUNNON FLJTEXdOOOOl* NANASOI* RPRNfZAOOl* 

SPACRt 101 

BNl 

A3 






BNl 

AA 

JET/EXT FLOW FIFLO CUNNON SECTION 




BNl 

A5 


c 


c 

c 

c 

c 

c 

c 

iO 


20 


30 


50 

C 

C 

C 




BNl 

Ab 

COItnON DELJ f PELF • tJET » lEXT 

» INAXJ • iNAXE 

*BN1 

A7 

1 KPAXJ » KHAXE » NPJET » NREXf 

» PSJET » PSEXT 

*BN1 

AO 

2 PI » P? . P3 > PA 

» USTJ » USTE 

*B51 

A9 

3 VJCT » VrXT f XJET > XEXT 


BNl 

50 



BHl 

51 

connoN ENOCH 


BNl 

52 



BBl 

53 

END OF COHHUN TU BE COPIED TO RtSTARI FILE 


BNl 

5A 



BNl 

55 

cunnuN Fion«5i» lsuoniig) 


BNl 

56 



BNl 

57 



B.Nl 

50 

CONTINUE 


BNl 

59 

IF (IPPFSS.F3.0I GO TO 20 


BNl 

60 

CALL 60ATII < 1 » X * P S I ( 11 » INA X J. KNAX J. I JE T» VJ E T 

»XJET»PS JET* Yl.Plf 

DUN ONI 

61 

1.01JN.PS JET.F lO.MRJETI 


BNl 

62 

CALL BOAT 1 1 (2»X>PSI (NPSI l»IHAXE>KHAXE»IEXT>VEXT»XbXTf PSEXTfY2»P2»BH 

63 

lOun.oun.PSjE t* fio.npex t i 


BNl 

bA 

CONTINUE 


BNl 

65 

00 50 I«1»NPSI 


BNl 

66 

WTNI X( I 1*0.0 


BNl 

67 

oa 30 J>lfNS 


BNl 

6B 

WTNI X( II«wrflIX( I )»AIPHA(J. I ) 


BNl 

69 

IF ((PRESS. EO. 01 GO TQ AO 


BNl 

70 

IF (PSim.LT.PSJETI CALL BOATlI (1>X»PS1(1)» 

INAXJtXnAXj*! JET*VJETBN1 

71 


1 »XJE r»PSJET»0UM»P»0UN>0Un*P5JET»FI0»NRJET» BHl 72 

IF IPSl (I I .GE.PSJETI CALL POATIt t2» X» PSI 1 1)» IPAXE »RNAXE > IE XTt VEXTBfU 73 


»XEXT»PSExr»OUN» P» DUN* DUN* P S JFT* F 10»NR E X Tl 

BNl 

7A 

CONTINUE 

BNl 

75 

RHO( I l*PTB9517.501/T(I l/WTNIX(I 1 

BNl 

76 

RHOOUTI n>L.^O( II/1.9A 

BNl 

77 

DO 70 I*1*NPSI 

BNl 

7B 


BNl 

79 

FREE STREAM VELOCITY WILL BE SET TO 1.0 FPS IF ZERO IS ENTERED 

BNl 

BO 


BNl 

B1 

U(ll*ANAXl(1.0*U(in 

BNl 

«2 

CPBARd 1-0.0 

BNl 

83 

HSTATI I 1*0.0 

BNl 

BA 

TX-TII 1 

BNl 

85 

CALL BOATTK ( TX * TTB * I TK F Y* SDT*H0T*NT1 

BNl 

86 

IF (ITKEY.EO.Ol GO TO 250 

BNl 

87 

DO 60 J«1*NS 

BNl 

88 

CALL BOATLP ( I T«F Y , J,2 . SOT *HOT* AX » 

BNl 

89 

HCP(2» JI«AX*A5055.31 

BNl 

90 



CALL BOATLP t t TKE V» J» A » $07 *MOT» AX ) 

BNl 

91 


HCPUji JI«AX*63036.31 

Bni 

92 


HSTATM l-HSTAri I l«HCP(l» JlAALOHAU.n 

BNl 

93 

60 

CPBAR 1 1 |>CPBAR ( I I^HCp ( 2» J I • AL PMA f J » I ) 

BNl 

95 

70 

CONTINUE 

SNl 

95 


IF « IPPESS.EO.OI Yni>PSIf 1)/SQRTIRH0I11*U( III 

BNl 

96 


IF UPRiSS.NE.OI YUl.Yl 

BNl 

97 


TfOP-YC HPS I I 

BNl 

9B 


O'J BO I« 7 »nPSt 

BNl 

99 


oup-(RHO( I ipun MRHnn- 1 1*01 1-111/2 

BNl 

100 

80 

Tm-SQKT(Y(I-ll«*2*«PSm I♦♦?-PS1 f l-n*«2) /DUN) 

BNl 

101 


IF (YTOP.EU.O.I 00 71 90 

BNl 

102 


XK2» CYI1PSI I-Y70PI /OY 

BNl 

103 

90 

CONTINUE 

BNl 

105 

C 


BNl 

105 

C 

KE2 INITIALIZATION 

BNl 

106 

C 


BNl 

107 


IF 1 IVIS.GE.Ol GO 1C 150 

BNl 

108 


IF (X.GT.XiNtn GO TO 150 

BNl 

109 


IF IIVIS.EQ.-2I GO TO 110 

BNl 

110 

C 

XU NUT SPECIFIED* IVIS--1 

BNl 

111 


IV1S>0 

BNl 

11? 


CALL BOATVI 

BNl 

113 


XKlll «0. 

BNl 

115 


XK INPSI l»0. 

BNl 

115 


XE 11 )>0. 

BNl 

116 


XEINPSII*0. 

BNl 

117 


OU 100 I-2*NPS1 

BNl 

IIB 


XK II {•XHUI 1 l*IUf t»l l-UM-l ) 1/2. /OELPSIPUI I)*Yn>/.3/PSIfII 

BNl 

119 


XKU l•ASSIXKIIII 

BNl 

120 

100 

CONTINUE 

BNl 

121 


lVIS—1 

BNl 

122 


GO TO 130 

BNl 

123 

110 

CONTINUE 

BNl 

125 

C 

XR SPECIFIED* IVIS--2 

BNl 

125 


XNUm-0. 

BNl 

126 


XNUfNPS I l>0. 

BNl 

127 


00 120 I-2»NPSI 

BNl 

12B 


XNUII l•.^•XK 1 1 MPS t 1 1 l*2.*0ELPSI/U f 1 l/Ytl 1/ fU f 1 * 11-U 1 1 -1 1 1 

BNl 

129 


XNUni>A8SlX<1Uim 

BNl 

130 

120 

CUNT INUE 

BNl 

131 

130 

CUNTIHUF 

BNl 

132 


DO 150 I>2*NPSI 

BNl 

133 


XE II l•.09♦RHOf I MXX (II**2/XNUm 

BNl 

135 


RXKI II «XKI 11 

BNl 

135 


1^0 


no 

c 

c 

c 


160 


170 


1«0 

C 


190 

ZOO 

?10 

ZZO 

Z30 


RlEtlMX€IIl 
COHf INUE 
RXRI 1 |«0. 

RXEIll-0. 

RXK|t 1 PSn- 0 . 

PXEfMPS I)’'0. 

CONTINUE 
CALL aOATVI 

CHECK OIFFUSIUN STEP 5 I 7E 

XO-triNPSl l-Y( m/FLnATtNPSIl 
00 160 t>Z»NPSI 
0UNnr-?.*|TII»l l-Y(l-ll) 

X0*AN|N1 IXOvOUNHYl 
rOL«FOL*OFOL 

IF CFOL.GT.l.. OR.FOL. lt. .Oil F0L«1. 
X0>X0*F0l 


iNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

SNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

B.11 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 

BNl 


on 170 I*Z«NPSI 
0UNNT*A( mi»A< t-1 MA( IMAft I 
0UNNY>0UHNY7Z« 

OUNNY*PSnmOELPSMOELPSI*STr.NAIl)/XLEII)/OURNY/l.S 

XO*ANINlf XOvDUNNYI 
0X>ANINl<0«fX0l 
IF f IPRFSS.EQ.OI GO TO 1 BO 
XOUN-X«OX 

CALL BOAT I I ( 1 » XDUN»PS I ( 1 I » I N AX J»KNAX J» IJE T » V JE T » XJE T » PS JE T> DUN* P3BN1 
1*RTI1I*RU( 1I*PSJET*FI0.NRJFT| BNl 

CALL BOAT 1 1 fZ>XOUN*PSI <NPSn»tNAXE>KHAXE*IEXl*VEXT*XEXT*PSEXT»DUNBNl 
l*P6*PTf NPSII»KUtHPSf I»PSJET*FI0»NPEX1 } r.i 

CONTINUE 

BNl 
BNl 
BNl 
BNl 
BNl 
BNi 
BNl 
BNl 
BNl 
BNl 
BNl 
BNl 
BNl 


INOtZlO* 190 
Z 00 *? 20.??0 


CALL BOATS? 

lOUT-IOUT^l 
IF IIFINISI 
IF |X«XNAX1 
IF (PRNT-PCNTI ?30*?30*?10 
CONTINUE 
GO TO ZAO 
IFfNIS-Z 
XN^^>Z. OAXNAX 

Call buatut 

IF IPSIO.NE.O.I write f6»Z60l 


136 

137 
13B 
139 
lAO 
161 
lAZ 
163 
166 
166 
166 

167 

168 
169 
160 
161 
15Z 
163 
166 
155 
166 
157 
16B 

159 

160 
161 
16 Z 
163 
166 

165 

166 
167 
16B 

169 

170 

171 
17 Z 
173 
176 

175 

176 

177 
17B 

179 

180 


o o r> 






2<*0 


250 

C 

260 


270 


PCNT«0.0 
CONTINUE 
" IN-IFINIS 

(LSWONKil) WRITE (6*270) OX*X 
CAU. B0ATM2 

CHECK FOR SWO UR MAX TIME EXCEEDED HERE 
ITFlG-0 

IF (LSWQNUn ITFLG-1 
TIME-ETIME (TSTRT ) 

IF (TIM6.LT.RTJAC) TIME «T TME+24 .♦60. 

TJAC-TIME-RTJAC 
IF (TJAC.GE.RTJOB) ITFLG«2 
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GO TO 10 
STOP 
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END 
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CONTINUE 
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IF (PSIO.NE.O.) CALL BOATOS fPSIEI 
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DO 60 I-2*NPSI 

PSIT-PS II*FLOAT( I-l )*0ELTSI 

IF <PSI T.LE.PSI (K n GO TO 30 
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GO TO 20 
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IF (I.EQ.2) RATSI-IPSIT-PSIll n/OELPSI 
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IF (I,E 0 , 2 ) XK(NPSn>PXK (NIPSI l^ftATSIP^a 
IF (I.FU.2) Xf (NPSI MPXF(NPSI )*RATSIP*2 
CONTINUE 

oa 50 J.1»NS 

c ON nii J ^ ^ ^ ‘ ‘ ‘ j 
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U( I ) «RU( I } 
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GO TO 80 


CONTINUE 
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GO TO l<iO 
CONTINUE 

WRITE ( 6 »I 50 ) IfTnitRTin 

GO TO 120 

CONTINUE 

IF (LSWONI^n WRITE ( 6,1601 
IF (OX.LT.OxniN) GO TO 130 
0 X«DX/ 2.0 
GO TO 160 
WRITE ( 6,1701 
CALL EXIT 
RETURN 
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SUBROUTINE BOATOT 
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6 J5 » LHAI.F 

7 NPSl # NP 

8 P # PCNT 

9 OX # RALPHA 
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OIHENSION RPRMV(10I» DATA(7) 
data ZCn/^HCO / 

DATA ZC02/^HCJ2 / 

DATA ZHZ0Z<»MH2D / 

DATA ZN2/^HN2 / 

CALL DATE IROATE > 

OQ 10 1»1»NS 

IF (Atom .FQ.ZCU) icn«i 
IF fAIDIII.E0.ZCU2) ICO?»I 
IF (AID! I ) .bO.ZHPO ) TH20*! 

IF I AIDU ) .E0.ZN2) IN2»I 
CONT INUE 
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OJ 20 I«1»MPSI 

YUUr II )>Y( n/RJ 
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WRITE (6»260) X # < T I TL E ( I ) » I • 1 » 18 ) » I PAGE * RDA TE 
WRITE f6«220) 
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POUT-P/2117.0 

WRITE 16*230) xnPj.DX.PCUT 
WRITF (6«2<»0) 

WRITE (6*370) 

OfJ 30 I«l»MPSt 
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HOUT-HSTAT (I )/^ 5055. 3 1 
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CONTINUE 
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BOT 

164 

160 

IF ILL) ie0»170,180 

BOT 

165 

170 

N«N*I 

BOT 

166 


LL-1 

BOT 

167 


GO TO 120 

BOT 

168 

180 

CONTINUE 

BOT 

169 

190 

CONTINUE 

BOT 

170 


RETURN 

BOT 

171 

C 


BOT 

172 

200 

FORMAT (E10.3» II0*60X 1 

BOT 

173 

210 

FORMAT (7E10.3>10X) 

BOT 

174 

220 

FORMAT llHO»0Xf 3HX/R»ex»8H0ELTA X > 4HF E E T> 4 X» 1 OHPR ES S 1 ATM 1 ) 

BOT 

175 

230 

FORMAT I4X»6E14.6) 

BOT 

176 

240 

FORMAT I4H0 PT»5X, 3HV/R, 6X,8HVEL0C I T Y, 4X , 1 IHT EMPER A TORE t 5X» 5H0EMSI BOT 

177 


1#2HTY.7X,8HMACH ND. »5X »0HENTH-TKEf 5X,9HVISC0S ITY .9X , 3HPS I » 10X» 5HPRB3T 

17B 


2-XE 1 

BOT 

179 

250 

FORMAT (3HOPT»3X»5H f/R » 7 1 3X > A4 » 6X 1 » 1 X > 3H PT 1 

BJT 

IPO 


?60 

FORMAT 

270 

FORMAT 

280 

FORMAT 

290 

FORMAT 

300 

FORMAT 

310 

FORMAT 

320 

FORMAT 

330 

FORMAT 

360 

FORMAT 

350 

FORMAT 

360 

FORMAT 

370 

FORMAT 


lX#9HLB/FT/StC I 
END 


{ IHl, //// ^/3H X»E15,7»5M F E E T »6 X » 1 HA<t » 6X » <tHP AGE 1 4 > 7X » AlO ) 
(I3»F9.5» 7E13.5 » 

(IH »//‘tOX#2PMPEACTinN RATES C NaiE 7HL- SE C I / / » 

(I <. » F 1 0 . <» » 9 E 1 , 6 I 
( 1H0» <»<^X» KiHhOlF FRACTIONS! 

( 1H0» 35X» 36HNE T RATE OF PROOUCTICN ( W-OOT /RHO*U !l 
l3HOPT*8X*7l3X»A<i»6>XI 1 
( l 3 »qx, 7 E 13 . 5 »I 3 ) 

UH0»7HPT><fX»3HY/PreX»5(BHRFACT10N# I3» IIX l*7HPT I 
(9X»5( 10X*7HPP>9Xf 7HRN!) 

CI3»1X*11E11.A,IA) 

I10X»8HFEET/SFC»^X»11H K »6X,5HG«/CC*2^X,6HCAL/GH»6B3T 

BOT 
BOT 


BOT 

BOT 

BJT 

BOT 

BOT 

BOT 

BOT 

BOT 

BUT 

BOT 

BOT 


161 

18? 

183 

18 A 

189 

186 

107 

188 

189 

190 

191 
19? 
193 
196 


1 


c 

c 

c 

c 


c 

c 


c 

c 


SUBROUTINE BOATRS ( XTNAX* XTJOB I 
BOATRS - RESTART FNTRY POINT 


BRS 

ORS 

BRS 

BRS 

BRS 

BRS 

BRS 

BRS 

BRS 

BRS 

BRS 


LOGICAL LHALF 


^ON 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


fALOCIlfll* CNIUin 

(Ecc» cp(i>i3n» frouTt cHiuisn 

ICPTDV» CGHVI 

CHCpii*n* wNiiii 
ISlARTf TTB(ll) 


1 

2 

3 

5 

6 
7 

a 

9 

1 

2 

3 

i 

6 

7 


BUATCM - BOAT CUNNPN 

OINENSION AISOl* AIOI25)* ALnCfS0>6l» ALPHA I 2 3 » SO ) » CGHVf7SO#3>» 
CPI2S»25I» CP0AR(50»» CPTBVI750I» ECCI50). GI25I» GTBVI750)» 
HCPI2»25I# H'-TATISni, HTBVI 750I. IRRf<iO)» IRRRI40»S>* IRTC^O)> 
ISAVEI6I* J123<>S(S), PSEXTISO). FSIfSO}* OXI25)» RALRMA 12 S» SO 1 » 
RCI<iO#3>» RHOfSOI» RHtirUTI50»» PTCSOI# RU I SO I » RXEISOI# RXKI50»» 
SIGHAI50»» START(1)» T ( SO » , TlTLEflSI* UCSO)» VE X T ( 4 » 2S» 2 ) » 
VJET(«»*2S»2»» WOOTf25»50l» WNI25I» WP(2S)» NTNlXISOIr WTNOL E 1 2S )» BR S 
XE(SOI» XFXTISOIf XJFTISOI* XKfSOU XLEfSOI* XHU(SOI» Y(50l» BRS 

YCUTISOI# ZI0»S» BRS 

BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
»BRS 
»3RS 
*BRS 
>BRS 
*BRS 
»BRS 
»BRS 
»BRS 
#BRS 
»BRS 
»BRS 
»BRS 
»BRS 
»BRS 
»BRS 
• BRS 
BRS 
BRS 
BRS 
BRS 


CUNNQN 

COMMON 


TTB(30) 

HFI25J 

9 

CPTBV 

9 

GTBV 

9 

HTBV 



A 

A ID 

9 

ALPHA 

9 

CARB 

9 

CM 

9 

CPBAR 

CRR 

CVISC 

9 

OELPSI 

9 

OFOL 

9 

OPOX 

9 

OX 

OXMIN 

FDL 

9 

FFF 

9 

G 

9 

GGG 

9 

HSTAT 

lOELP 

IfCC 

9 

IFINIS 

9 

lOUT 

9 

lOUTl 

9 

I0UT2 

IPAGE 

IPRESS 

9 

IRR 

9 

IRRR 

9 

IRT 

9 

ISAVE 

ITFLG 

IVIS 

9 

J1 

9 

J2 

9 

J3 

9 

J4 

JS 

LHALF 

9 

MMOO 

9 

MPSI 

9 

MXNPT 

9 

MXNPl 

MPSI 

NR 

9 

NRAO 

9 

NRAS 

9 

NS 

9 

NT 

P 

PCNT 

9 

PRNT 

9 

PRNTXC 

9 

PS I 

9 

PSIO 

QX 

RALPHA 

9 

RBUOY 

9 

RC 

9 

RHO 

9 

RHUOUT 

RJ 

RT 

9 

RTACU 

9 

RTJAC 

9 

Rf JOB 

9 

RTMAX 

RU 

RXF 

9 

RXK 

9 

SIGE 

9 

SIGK 

9 

SIGMA 

T 

TC'JNT 

9 

TEOGE 

9 

TEMKM 

9 

TEMRP 

9 

TITLE 

TKINET 

LI 

9 

UNIT 

9 

MOOT 

9 

WM 

9 

WP 

WTMIX 

WTMOLE 

9 

X 

9 

XCHANG 

9 

XD 

9 

XE 

XINI T 

XK 

9 

XK2 

9 

XLE 

9 

XMAX 

9 

XMU 

Y 

710 









FLJTEXf lOOOOIf HAMaSI3>» 

RPRMI2400I 

9 

SPACRIIOI 



COMMON 


JET/EXT FLOW FIELD COMMON SECTION 


1 

2 

3 

*1 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


ooo o o o o oo n o o 


fin. if 


c 

10 

20 


COMMON 

OELJ 

$ 

DFLE 

• 

IJE I 

P 

lEXT 

* IMAXJ 

> IMAXE 

1 

KMAXJ 

P 

HMAXF 

* 

MPJET 

P 

NREXT 

* PSJET 

* PSEXT 

2 

PI 

P 

P2 

f 

P3 

P 

P6 

* USTJ 

* USIE 

3 

VJE T 

P 

VFXT 

f 

X JF 1 

P 

XEXT 



COMMON 

ENDCM 









END OF 

COM MUM 

TO 

BF COPIED 

TO RESTART 

FILE 




COMHIld FI0(3.5I> IFHtnfll, LSWUNfl6> 


Tn*X>XT*1AX 
TJU»*XT JOB 
CALL DATE fROATE ) 

LOOK UP FILES AND READ IM COMHOM 

CALL BOATLU 111 
NWCN>NFWABfSrART*ENOCm«l 
BEAD (1) ( STARTI n» I«1 tNWCPI 

WRITE 16*101 T1UE*R0ATE *IFNAP*NAnAS 

GET TII1E PARAMETERS 

if I THAX.GT.O.OI RTMAX>RTMAX*Tf1AX 
IF (TJOR.GT.O.OI PTJOB-Tjnn 
IF I TJOB.EO.O.O) RT JIB.RTMAK-RTACU 
RTJAC»0.0 

WRITE ( 6*201 RTnAX*RTACU*RTjnn*RTJAC 
RE TURN 

FORMAT I1H1*23H RESTART OF BOAT RUN - * IBAA * 5X » A 10/ /16H RESTART 
ILE - »3A2*21H* BUAT-SPFCRA FILE - .3A2/I 
FORMAT I/15H MAX HUN TIME • * F8 . 2 > 6X *22H ACCUMUL ATFO RUN TIME •*F8 
1/15H MAX JOB TIME « * FB . 5X * 22HACC UMUL A TE D JOB TIME >*FM.2//I 
END 


BRS 
*BRS 
*BRS 
*BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
BRS 
FIBRS 
BRS 
• 2BRS 
BRS 
BRS 


A6 

67 

6B 

69 

60 

51 

52 

53 
5*. 

55 

56 

57 

58 

59 

60 
61 
62 
63 
66 

65 

66 

67 

68 

69 

70 

71 

72 

73 
76 

75 

76 

77 

78 

79 

80 
81 
82 
83 
86 



SUBROUTINE BUATSL (X>A>N) 

BSL 

1 

c 

THIS PRUGRAM FINOS TMF SOLUTIONS TO A SET OF 

N simultaneous linearbsl 

2 

c 

EQUATIONS BY USING THE G AUSS-r.'jRDAN REUUCTION 

ALGORITHN with the BSL 

3 

c 

DIAGONAL PIVOT STPATfGY 

BSL 

4 


DINENSION A(25»?5I» X(?5I 

BSL 

5 


O’) 40 K»l»N 

BSL 

6 


IF IABSU(K»K)I.GI.I.E-10) GO TO 10 

BSL 

7 


WRITE (6^601 

BSL 

B 


G'J TO 50 

BSL 

9 

10 

KP1*K*1 

BSL 

10 


00 20 J-KPl.N 

BSL 

11 

10 

AtK» J)«AIK»JI/A(K»K I 

BSL 

12 


XIK)*X(K) /A(KrK ) 

BSL 

13 


A(K»K )«1.0 

BSL 

14 


DO 40 I»1.N 

BSL 

15 


IF 1 I .EO.K.OR.AI I>X ).E0.0. I GO TO 40 

BSL 

16 


DO 30 J>KP1»I4 

BSL 

17 

30 

A(1>J)-A(I»J)-ACI>P )«AfK».l) 

BSL 

IB 


XlI)«Xm-AII*K)«X(«() 

BSL 

19 


AC I»R)>0. 

BSL 

20 

40 

CONTINUE 

BSL 

21 

50 

CONTINUE 

BSL 

22 


RE TURN 

BSL 

23 

C 


BSL 

24 

60 

FORNAT I22M ERROR SPALL PIVOT 1 

BSL 

25 


END 

BSL 

26 


r 


c 

c 

c 


c 
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SUBROUTINE BOATS I 

BOATSl - BUAT INPUT ROUTINE 

LOCtCAL LRSI 
BUATCN - BUAT CUNHCN 

DIHENSION A(*»OI» AIOCP5I* AL0C<S0»6I» AL PHA ( 0 5# 50 I* CGHVf750t3lf 
I CHI?5,?5I, CPBAR(50I* CPTRVt750», ECCtSOl* G!75lt GTBV(750I» 

^ HCP|?»25I» HSTAT(50)» HTBV(750It 1RR(^0I> IRPRU0>5>* IRTtAOI# 

3 ISAVEI6)» J173^5(5I> PSE)fTr50l» PSII50)» QX(751» R ALPHA ( ?5» 50) » 

A RCtAOfll* PH0I50)» PHn0UTI50)» RTI50I* RUI50)* RXEf 50 )» RXK150I* 

5 SIGNAf50l» STAKnil. U 50 I . TIUEIIB)* UI50I» VEX T I A » Z5> 2 ) » 

6 VJETI A»25»?)t WDUr (P5»50) » WN(25)» UPI25)* WTNIXI50)* WTNOLE (251 

7 XE(50)» XExn50)> XJEn5Q)» XK(50)» XLE(50)» xnU(50l» YISOI# 

8 YCUT(50)» ZI0I5) 


;CICAL LHALF»LSWUN 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

equivalence 

EQUIVALENCE 

EQUIVALENCE 


CH 


(Jl» JI?1A5(1 
(ALUCdtl )» 
IECC» CPIlf)3 
(CPTOVf CGHV) 
IHCPtlfll* WN 
(START* TTBd 


connoN 


1 

2 

3 

A 

5 

6 

7 

8 
9 
1 
2 
3 
A 

5 

6 
7 


) 

d*ii) 

I* (TOUT* 

111 

) 


CHd>15d 


TTB(3Cn 

HF(25I 

f 

CPTBV 

9 

GTBV 

9 

MTBV 


A 

AID 

9 

ALPHA 

# 

CARS 

9 

CN 

9 

CRR 

CVISC 

9 

OFLPSl 

9 

OFOL 

9 

OPOX 

9 

DxniN 

FDL 

9 

FFF 

9 

G 

9 

GGG 

9 

lOELP 

lECC 

9 

ifinis 

9 

I OUT 

9 

lOUTl 

9 

IPAGE 

IPRESS 

9 

IRR 

9 

IRRR 

9 

IRT 

9 

ITFLC 

IVd 

9 

J1 

9 

J2 

9 

J3 

9 

J5 

LHAL F 

9 

P*«00 

9 

NPSI 

9 

NXNPT 

9 

NPSI 

NP 

9 

NRAO 

9 

NRAS 

9 

NS 

9 

P 

PCNT 

9 

PRNT 

9 

PRNTXC 

9 

PSI 

9 

OX 

P ALPHA 

9 

RBUOY 

9 

PC 

9 

RHO 

9 

RJ 

PT 

9 

RTACU 

9 

rtjac 

9 

RTJOB 

9 

PU 

p xc 

9 

RXK 

9 

sige 

9 

SIGN 

9 

T 

TCONT 

9 

TEtIGE 

9 

TENRN 

9 

TENRP 

9 

TK i*«et 

U 

9 

UNI r 

9 

MOOT 

9 

WN 

9 

WTNIX 

WTNPir 

9 

X 

9 

XCHANG 

9 

xo 

9 

XINIT 

XK 

9 

XK2 

9 

XLE 

9 

XNAX 

9 

y 

zro 








FLJTEXdOOOO)* NANAS(3)» 

RPRN(2A00I 

9 

SPACRdO) 


CPBAR 

OX 

HSTAT 

I0UT2 

ISAVE 

JA 

NXNPl 

NT 

PSIO 

RHOOUT 

R1NAX 

SIGNA 

TITLE 

WP 

XE 

XHU 


BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

BSl 

• BSl 
BSl 
BSl 
BSl 
BSl 
BSl 
BSl 
BSl 
BSl 
BSl 
BSl 

asi 

BSl 

BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 

• BSl 
>BS1 

• BSl 

• BSl 

• BSl 
BSl 
BSl 
BSl 


1 

2 

3 

A 

5 

6 
7 

a 

9 

10 

11 

12 

13 

lA 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3A 

35 

36 

37 
3B 
39 
AO 
A1 
A2 
A3 
AA 
A5 


iP« 






"f'JI!:"’*''' 


c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

0 - 

V 

c 


c 

c 

c 

10 


c 


JET/EKT FLOW FIELD COPHOM SECTION 


COMMON 

OELJ 

» 

DELE 

9 

I JET 

9 

lEXT 

» IHAXJ » 

IMAXE 


AMAXJ 

9 

KMAXF 

9 

NRJET 

9 

NREXT 

» PSiET • 

PSEXT 


PI 

9 

R? 

9 

P3 

9 

PA 

> USTJ > 

USTE 

1 

VJt T 

9 

VFXT 

9 

XJET 

9 

XEXT 



COMMON 

ENOCH 









EMD OF 

COMMON 

ro 

BF COPIED 

TO RESTART 

FILE 



COMMON 

FIOI 3. 

b% 

. TFNAMI31 

» LSU0N(16> 




OINFNSIIH CPT«Un. HTPdl* GTBIll 
otnENSiuN ixnoi 

EOUIVALLNCE I *L OC 1 1 » I I »C PTB 1 1 » I * f *LOC C 1 » *HTB I 1 1 I 
EJUIVfcLENCE UL0CU»31»GTSm I 
EQUIVALENCE ( I T JOB* PT J RB I 

OINENSION RINI50I* HniOI50»* H0LD2I501* MOLD3I50I 
DAT* LRSl/.FALSE./ 

CALL JPARAHS UK» 

ITJUi-IXIll} 

RTjaB>RTJ0B/60.0 
RTJUB«AN*X1 IRI JOB-9. 

READ n»AOOI IITPE> IFNAr,NANAS»RTNAX»RlJOB.LSWON 
IF (ITYRE.NE.O! GU TO 10 

SAVE TINE PARAMETERS FOR AFTER RESTART 

RESTART - COPY COMMON 

CALL 80ATRS (RTHAX»RTJOP I 
GO TO 330 

INITIALIIF Ti’^t KEFPING VARIABLES 
CONTINUE 

IF IRTMAX. rO.O.OI *TMAX«10.0 
IF CRTJ08.E0.0.01 RTJOB-RTMAX 
RTACU-0.0 
RTJAC>0.0 

IFINIS-0 
LHALF>. FALSE. 


asi 

AB 

BSl 

A7 

»BS1 

A8 

»BS1 

AQ 

»BS1 

50 

BSl 

51 

BSl 

52 

BSl 

53 

BSl 

5A 

BSl 

55 

BSl 

56 

BSl 

57 

BSl 

5B 

BSl 

59 

BSl 

60 

BSl 

61 

BSl 

62 

BSl 

63 

BSl 

6A 

BSl 

65 

BSl 

66 

BSl 

67 

BSl 

6B 

BSl 

69 

BSl 

70 

BSl 

71 

BSl 

72 

BSl 

73 

BSl 

7A 

BSl 

75 

BSl 

76 

BSl 

77 

BSl 

7B 

BSl 

79 

BSl 

SO 

BSl 

01 

BSl 

02 

BSl 

B3 

BSl 

BA 

BSl 

85 

BSl 

06 

BSl 

87 

BSl 

BO 

BSl 

89 

BSl 

90 




READ 15*370} t TI HE ( tl . I-l* 18 } 

BSl 

91 



READ 15*340) nPSl*NnRSI*NS»NR,Nr»tDElP*IRRESS*lVlS*lHAXJ»KNARJ* 

IHABSI 

92 



lXE>KnAXE*lGim* I0UT2.NRA0 

BSl 

93 



OELPSI-IOELP 

BSl 

94 



nxNPT-npsi 

BSl 

95 



CALL SFVFL ( 0 .0* AL PHA» ?5*HXNPT ) 

BSl 

96 



MXMPl»«XNPT-l 

BSl 

97 



NPSI-HPSI-1 

BSl 

9B 



READ <5*390) X *R J* XHAX*PRNT » XCHANG* PRNTXC*F DL *DFOL 

BSl 

99 



XINIT-X 

BSl 

100 



READ <5*390) XLE < L ) *SI G)1A< 1 ) * TCONT *TKINET*C ARBUN*CNZ INT* C VlSC 

BSl 

101 



READ <5»390) P* U < 1 ) *U ( PPS I ) » T < 1 1 * T 1 MPS I ) 

BSl 

102 



OXniN«l.E-10 

BSl 

103 



DX«.1«RJ 

BSl 

104 



READ <5*390) FF F * GGG* PS 1 0* D<E L J * OE L E *UST J* US TE * RBUOY 

BSl 

105 


c 


BSl 

106 


c 

LOOK UP FILES 

BSl 

107 


c 


BSl 

108 



CALL BOATLU <0) 

BSl 

109 



CARB*CARBUN/<.78973-CNZINT ) 

BSl 

110 



IF <TKINET.E0.0.0) TKINET*400.0 

BSl 

111 


c 


BSl 

112 


c 

THE VALUE OF 30 SECONDS IS TO ALLOW FUR COMPILE TIME 

BSl 

113 


c 


BSl 

114 



UNIT-um 

BSl 

115 

r\. 


IF IDELPSn 20*20*30 

BSl 

116 

20 

READ <5*390) < AL PHA < J * 1 ) * J «1 * NS ) 

BSl 

117 



READ <5»390) < ALPHA < J* MPS I > » J • 1 * NS ) 

BSl 

118 

t 


MM00»«PSI-2 

BSl 

119 

r-' 


LRSI-.TRUE. 

BSl 

120 


GO TO 60 

BSl 

121 


30 

CONTINUE 

BSl 

122 



LRSI-.TRUE. 

BSl 

123 



READ <5*390) <RI N< I ) * I >1 *MPS I ) 

BSl 

124 



00 40 I-1*MPSI 

BSl 

125 


40 

RIN<I )-RIN<I )*KJ 

BSl 

126 



READ <5»390) < T < I ) * I«1 *MPS I) 

BSl 

127 



PEAD <5*390) <U<I)*I«1»MPSI) 

BSl 

128 


C 


BSl 

129 


c 

READ IN INITIAL XK PROFILE 

BSl 

130 


c 


BSl 

131 



IF UVIS.E0.-2) READ <5*390) < XK< I )* t>l*MPS I) 

BSl 

132 


c 


BSl 

133 



00 50 I«1*HPSI 

BSl 

134 


50 

READ <5*390) <ALPHA U* I )*J«1*NS ) 

BSl 

135 


NEW THERMO DATA DATA INPUT IN JANNAF TABLE FORM 


C 

C 

C 

60 


70 


00 

90 


100 

no 

120 

130 

C 

C 

C 

C 


C 

C 

C 

C 


DO 90 I-l»NS 

READ I5*360» AID 1 1 1 *WTM01E tl ) »HF 1 1 ) 
DO 70 IT«1»NT»2 


l#GTBnT*lI»HTBnT»n 
GTBf IT I «-GTBn T )♦! TPI I Tl ♦HF( 1 1*1000. 

GTBIlT*ll«-GTBtn*l ) *T TB ( I T* 1 1 *HF 1 11*1000. 
HTBUTI-CHTBIITMUFCI) 1*1000. 

HTBUT*! »■ IHTBUT + I l*HF ( 1 1 1*1000. 

CONTINUE 

NRECS«< 1-1 1*30*1 

CALL SFVMV ICPTB*CPTBV(NRECSIfNT) 

CALL SFVMV (GTB»GTBVIMRECSI »NTI 
CALL SFVMV IHTB * HTB VI NHF CS I »NT I 
IF (WTMOLE n 1-1 .01 P0>90»00 
lECC-I 
CONTINUE 

IF INR.LE.Ol GO TO 130 

REA0^(5>3801 (ZI0(Jl»J»l»9)*lRRm»IRTn)»IRCmK)»K*n3l 
DO 110 J>1>9 
IRRR n>ji«o 
00 110 L«1»NS 

IF IZI0< Jl-AIOU 1 1 110*100#110 

IRRRIl* Jl-L 

CONTINUE 

CONTINUE 

CONTINUE 


READ IN INVISCID DATA MAP 
IF IIPRESS.NE.OI CALL BOATIF 


IF CDELPSI.LE.O. 1 CALL BOATIP IRINI 
IF IIPRESS.NE.OI P«P^?117. 

R TO PSI MOO* - MAIN SECTION - 7/2/76 ARAP BT RAB 


IF I.NOT.LRSn GO TO 2B0 


CONVERT R TO PSI 

M0L021 1 1*P*SMPR2 (HS»AIPHA( 1* 1 1 # WTMOLE I / A2 *2 85 / T 1 1 1 
TEMZ-H0L02Ill*Um 


BSl 

136 

BSl 

137 

BSl 

13B 

BSl 

139 

BSl 

140 

BSl 

141 

BSl 

142 

BSl 

143 

BSl 

144 

BSl 

145 

BSl 

146 

BSl 

147 

BSl 

14B 

BSl 

149 

BSl 

150 

BSl 

151 

BSl 

152 

BSl 

153 

BSl 

154 

BSl 

155 

BSl 

156 

BSl 

157 

BSl 

158 

BSl 

159 

BSl 

160 

BSl 

161 

BSl 

162 

BSl 

163 

BSl 

164 

BSl 

165 

BSl 

166 

BSl 

167 

BSl 

16B 

BSl 

169 

BSl 

170 

BSl 

171 

BSl 

172 

BSl 

173 

BSl 

174 

BSl 

175 

BSl 

176 

BSl 

177 

BSl 

178 

BSl 

179 

BSl 

180 




iJiiliHiiiiiiriliiiiiMiiiif*' 


Tlliillililllfllllll'iillllimMlIr 




psillipi'" 









HaLO<l>«TEMZ*RINU 

BSl 

181 


IF 1 IFRESS.EO.OI GO TO 170 

BSI 

182 


lDUn«KnAXJ 

BSl 

183 

1^0 

FSA-FLOATC IOUM-1 ) / FI HA T 1 KH AX J-l »*PSJET 

BSl 

184 


CALL BOAT 1 1 ( 1 » X » P $ At 1 HAX J *KH AX J » I JE Tf V JE T » X JE T» PS JE T» YA» OUn»DUH» 

DBSl 

1R5 


lUMtPSJETtFIOtMRJfcT) 

BSl 

1B6 


IF C lOUM.EO.KiAXJ ) GO TO 150 

BSI 

187 


IF TRINI 1) .GE. YA) GO TO 160 

BSI 

188 

150 

10Un»I0UP-I 

BSl 

189 


PSB«PSA 

BSI 

190 


YB-YA 

BSl 

191 


GO TO lAO 

BSl 

192 

160 

HQLDm»PSA*IRINm-YA»/TYB-YAI*fPSB-PSAI 

BSl 

193 


HULOIl l«H0L[)in**2 

BSl 

194 


TEHZ»H0LDll)/(RINm**2T 

BSI 

195 

170 

CONTINUE 

BSl 

196 


IPITER-O 

BSl 

197 

180 

CONTINUE 

BSI 

198 


IF UPITER.EO.OI GO TO 100 

BSl 

199 


OUNPSI-HOLOl 1 1 

BSl 

200 


TENZ-IMOLOm/RIMl H**2 

BSl 

201 

190 

CONTINUE 

BSl 

202 


DO 200 I*2tNPSI 

BSl 

203 


TENN*TEI1Z 

BSl 

204 

C 

HUL02 CONTAINS THE DENSITY IN APPROPRIATE UNITS FOR PSI 

BSl 

205 

C 

ALPHA USED HERE AS HOLE FRACTIONS 

BSl 

206 


IF IIPITER.GT.OI P*P1*(H0L0( I»-0UNPSII/ CHOLOIMPSn-DUNPSI)*IP2-Pl»BSl 

207 


H0LD2I n-P*SHPP2(NSfALPHAIl»nf WTN0LE>/A2.2 85/T( 1) 

BSI 

208 


TEHZ»HDLD2(1 )*U(U 

BSl 

209 


DUN*>.5*<TEMM*TEMZ) 

BSI 

210 


IF f n.E0.2>.AN0«fIPITEP.GT.0n HDLDI 1) «H0L0l 1 ) ««2 

BSl 

211 


MUL0!II*H0LD( I-l)*DUN*IRINm**2-RINII-ll**2l 

BSI 

212 

200 

CONTINUE 

BSl 

213 


DO 210 I>1»NPSI 

BSl 

214 

210 

HOLOm>SQRT(HOLDII n 

BSI 

215 


IF (IPRESS.EO.O) 60 TO 2?0 

BSl 

216 


IF (IPITER.GT.O) GO TO 220 

BSl 

217 


CALL BOAT 1 1 ( 1 1 X tHOLO f 1) » I NAX Jt KNAX J> IJETt V JE T > X JE T» PS JE Tt DUPItPlt 

DBSl 

218 


lUNtOUNtPSJETtFIDtNR JET ) 

BSl 

219 


CALL BOAT 1 1 ( 2 1 X » HOLOf NP S M > INAXE* KNAX Ef 1 EX Tt VE X T> XE XT » PSEXT*OUN> 

PBSl 

220 


12>0UN»DUn»PSjET*FI0»NREXT) 

BSl 

221 


Pl-Pl/2117, 

BSl 

222 


P2-P2/2117. 

BSl 

223 


IPITER-1 

BSl 

224 


GO TO 180 

BSl 

225 





220 

CONTINUE 

BSl 

276 


IF IIPRESS.EO.OI PSJET«RJ*SORTCH0102I1 )*Un n 

BSl 

22 7 


RSIO-PSJET 

BSl 

22B 


YINTR«RJ 

BSl 

229 


IF UPRESS.NE.O) CALL BOATII ( 1 > X» PS JET» IHAX J > KNAX I JET* VJET» 

XJETBSl 

230 


1*PSJET* YINTR*0UM»0UN*0UM*PSJ6T*FI0*NREXT» 

BSl 

231 


IF UPRESS.NE.O) CALL ROAUI ( Y INTR* PSI D* R I N* HOLD* NPS I ) 

BSl 

232 

C 


BSl 

233 

C 

NOTE REDEFINITION OF DELPSI 

BSl 

236 


OELPSI»IHOLO(MPSI )-HOlO(I ) )/FlOAT(NPSI) 

BSl 

235 

C 


BSl 

236 

C 

INTERPOLATE INPUT VARIABLES FOR * PSI SPACING 

BSl 

237 

c 


BSl 

23B 


CALL SFVnV (U*HOLD?»NPSI ) 

BSl 

239 


CALL SFVNV 1 T»HOL03*NPS I) 

BSl 

260 

c 


BSl 

261 


IF IIVIS.EQ.-2) CALL SFVNV ( XK* H0LD1*NPS I ) 

BSl 

262 

c 


BSl 

263 


NSAV-HPSI 

BSl 

266 


IF IIDELP.NE.I) GO TO 230 

BSl 

265 


IF (NNPSI.EO.NPSI) GO TO 230 

BSl 

266 


HPSI-NHPSI 

BSl 

267 


NPSI-NPSI-1 

BSl 

268 


0ELPSI»0ELPSI*FL0AT(NSAV-1 ) / FLOAT « NPS I -1 ) 

BSl 

269 

230 

CONTINUE 

BSl 

250 

C 


BSl 

251 


00 2A0 I«1»HPSI 

BSl 

252 


PSKI )»HOLO<II ♦OELPSI*FLOATn-l) 

BSl 

253 


CALL BOATLI IPSI f I I *U(I ) *HOL0*HOLO2»HS AV) 

BSl 

256 


CALL BOATLI ( PS I < I ) * T (I ) *H0LD*HOL03*NS AV ) 

BSl 

255 

C 


BSl 

256 


IF UVIS.E0.>2) CALL BOATLI ( RS I f I >* XK 1 I )*HOLD* HQLOl * NSA V > 

BSl 

257 

C 


BSl 

258 

2^0 

CONTINUE 

BSl 

259 

C 


BSl 

260 


DO 270 J«1*NS 

BSl 

261 


DO 250 I«1*NSAV 

BSl 

262 


H0L02U l-ALFHAU*! ) 

BSl 

263 

250 

CONTINUE 

BSl 

266 

C 


BSl 

265 


OU 260 I>1*NPSI 

BSl 

266 


CALL BOATLI I PSI ( I ) » ALPHAU* I )*H0LD*H0LD2*NSAVI 

BSl 

267 

260 

CONTINUE 

BSl 

268 

C 


BSl 

269 

270 

CONTINUE 

BSl 

270 


S'iff . ,||- 


c 


c 

END OF nAiN SECTION » TO PSI HQO. 

BSl 

271 

280 

CONTINUE 

BSl 

272 

C 

initialize delta star calculation 

BSl 

273 


OELJ-OELE/RJ 

BS 1 

274 

C 


BSl 

275 


00 310 I>1«NPSI 

BSl 

276 


WTVR-0.0 

BSl 

277 


00 290 J»1*NS 

BSl 

2 78 

290 

WTVR.WTVR*ALPHA( J, I I•WTN0LFCJ1 

BSl 

2 79 


00 300 J«1»NS 

BSl 

280 


ALPHA! 1 1 « ALPHA Uf I IZWTVP 

BSl 

281 


RALPHAt I l>ALPHA( I » 

BSl 

282 

300 

CONTINUE 

BSl 

2B3 

310 

CONTINUE 

BSl 

2B4 


DU 320 I-1»NPSI 

BSl 

285 


Rum-ufii 

BSl 

266 


RTII »-T(I| 

BSl 

287 


XLEf I)«XLE(1) 

BSl 

288 

320 

SIGNA(Il«SlGNAm 

BSl 

289 


CALL BOATIN 

BSl 

290 


P-2117. 0*P 

BSl 

291 


OPOX-0.0 

BSl 

292 

330 

RETURN 

BSl 

293 

C 


BSl 

294 

340 

format *1515) 

BSl 

295 

350 

FORMAT (8F10.4) 

BSl 

296 

360 

FORMAT (A4>2X»7E10.3) 

BSl 

297 

370 

format (1BA4) 

BSl 

298 

380 

390 

format 

BSl 

BSl 

299 

300 

400 

FORMAT ni»lX»3A2,lX,3A2»lX,2F10.0»16< 1X,L1 >) 

END 

BSl 

BSl 

301 

302 



BSl 

303 




c 

c 

c 

c 


c 

c 


c 


c 

c 




SUBROUTINE BO*TS2 
BOATS? - reaction CALCULATION 


BOATCN - BOAT C INK ON 
OINENSION AfSO)# AlDf?9»» 


ALOCI SO* 


CN(?S*?5I* CPBARI50), CFTBVf750I# 
HCFI?*?5»# HSTAUSOI* HT0VC75OI* 
lSAVE(6)f Jl?3<i5(5)* PSEXT(50I* P 
RCf40»3>* RHOfSOIf RHOnunSOI* RT 
SIGHAISOI. STARTfll. TC50), TITLE 


6)* ALPHA(?S*SO)* CGHVI750»3)> 
ECC(50I» Gf?S»* GTBVfTSOI* 
IRR(40I* IRRRI40*SI* IRTUO)* 
SK50)* OX(?S)* RALPHAf?S»SO)* 
(SOI* RUISOI* RXEfSOI* RXK(SO)* 
(181* 0(501* VEXT(4*25*?I* 

* WP(?5I* WTNIX(SO)* WTN0L£(?5) 


VJET(«^>?S»?I* W0m(?5*50l» W«(25) 

XE(SOI* XEXT(SO}> XJET(SOI* XK(50I* XLE(50I* XNU(SOI» Y(SO|* 
YOUT(SOI* ZIDI5) 

LOGICAL LHALF*LSWOM 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

equivalence 

EQUIVALENCE 


(Jl* Jl?345(in 
(AL0C(1*1I* CN(t*lll 
(ECC* OMl*13ll* (TOUT* Cn(l*lS)l 
(CPTBV* CGHVI 

(HCPdfii* wnmi 

(START* TTB(ll) 


BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

*BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 

BS? 


CONNON 

CUNNON 


CONNON 


TTBOO 

A 

CRR 

OXNIN 

lOELP 

IPAGE 

ITFLG 

JS 

NPSI 

P 

OX 

RJ 

RU 

T 

TKINET 

WTNIX 

XINIT 

Y 

FLJTEX 


1 

? 

3 

« 

5 

6 
7 

e 

9 

10 

11 

1 ? 

13 

1^ 

15 

16 

17 

18 

19 

20 
21 
22 
?3 
?A 


$ 

HF(?5I 

0 

CPIBV 

0 

gtbv 

0 

HTBV 



BS? 

25 

p 

AID 

0 

ALPHA 

0 

CARB 

0 

CN 

0 

CPBAR 

#BS? 

?6 

0 

CVISC 

0 

OELPSI 

0 

OFOL 

0 

DPOX 

0 

OX 

*BS? 

?7 

0 

FOL 

✓ 

FFF 

0 

G 

0 

GGG 

0 

HSTAT 

*BS? 

?8 

0 

IfCC 

0 

IFINIS 

0 

lOUT 

0 

lUUTl 

0 

(OUT? 

*BS? 

?9 

0 

IPRfSS 

0 

IRR 

0 

IRRR 

0 

IRT 

0 

I SAVE 

*BS? 

30 

0 

IVIS 

0 

Jl 

0 

J? 

0 

J3 

0 

JA 

*BS? 

31 

0 

LHAIF 

0 

NHOD 

0 

NPSI 

0 

MXNPT 

0 

NXNPl 

»BS? 

3? 

0 

NR 

0 

NRAO 

0 

NRAS 

0 

NS 

0 

NT 

*BS? 

33 

0 

PCMT 

0 

PRNT 

0 

PRNTXC 

0 

PSI 

0 

PSID 

*BS? 

3A 

0 

RALPHA 

0 

R6U0Y 

0 

RC 

0 

RHO 

0 

RHOOUT 

*BS? 

35 

0 

PT 

0 

RTACU 

0 

RTJAC 

0 

RTJOB 

0 

RTNAX 

»BS? 

36 

0 

RXE 

0 

RXK 

0 

SIGE 

0 

SIGK 

0 

SIGNA 

*BS? 

37 

0 

TC'INT 

0 

TEOGE 

0 

lENRN 

0 

TENRP 

0 

TITLE 

»BS? 

38 

0 

U 

p • 

UNIT 

0 

WDOT 

0 

UN 

0 

UP 

*BS? 

39 

0 

WTNPLE 

0 

X 

0 

XCHANG 

0 

XO 

0 

XE 

»BS? 

AO 

0 

XK 

0 

XK? 

0 

XLE 

0 

XNAX 

0 

XNU 

*8S? 

A1 

0 

ZIO 









BS? 

A? 

100001* 

NANAS(3I* 

RPRN(?AOOI» 

SPACRdOl 


BS? 

A3 


JET/EXT FLOM FIELD COMMON SECTION 


BS? 

03 ? 


AS 



I 




r> n M o r> o r» n rt 






!l|l||Uriilll|i'[]NHW^ 




c 

c 

c 

c 

c 


10 


0 


COMHOM DELJ * DELE 
1 KHAXJ f KHAXE 

^ PI p P? 

3 VJET » VEXT 

CQHMOM ENOCH 


IJET 

NRJET 

P3 

XJET 


lEXT 

NRtXT 

XEXT 


IHAXJ 

PSJET 

USTJ 


* IHAXE 
» PSEXT 

# USTE 


END OF COHHUN TO BE COPIED TO RESTART FILE 
COHHON FI 0 ( 3 » 5 )» IFNAHm, LSW 0 NC 16 I 


R« 8?.06 
AV- 6 . 025 E 23 
DO 330 L- 1 #NPSI 
RRT» 1 ,<> 86 *T(L I 
ROOTT-SQRTITIL I) 

TX-T ILI 

CALL ROATTK I TX, TTB » I TKF Y, SOT,HDT, NT I 
IF fITKEV.EO.OI CALL FXIT 
CALL SFVFL ( 0 . 0 »WP»?B) 

CALL SFYFL ( 0 . 0 »WH,P 5 ) 

CALL SFVFL ( 0 . 0 >CH» 6 ? 5 I 
CALL SFVFL I 0 . 0 » 0 X» 25 } 

00 10 1 « 1 >NS 

CALL BOATLP f I TKE Y » I » 3 »S DT#MDT» AX ) 

G( 1 )*AX 
C INTI NO E 

reaction calculation 

REACTION KINETICS CONTINUE DOWN TO AOO DEGREES K 
UNLESS TKINET IS SET TO A VALUE OTHER THAN AOO K 
REACTION KINETICS FOR ALL REACTIONS CONTINUE DOWN TO TKINET 

IF tNR.LE.Ol 60 TO ?80 

IF (T(LI-TKINET) 280 *P 80»20 

CONTINUE 

DO 270 I* 1 *NR 

TENRP» 0.0 

TENRN- 0.0 

KK>IRTf I) 

REACTION CONSTANT TYPE 


BS 2 

» 8 S 2 

»BS 2 

«BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

8 S 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 

BS 2 


A 6 

A 7 

A 8 

A 9 

BO 

51 

52 

53 
5 A 

55 

56 

57 

58 

59 

60 
61 
62 
63 
6 A 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
BO 
B 1 
82 

83 

84 

85 

86 
87 
B 8 

89 

90 


- 





GO TO C30t40>50#00» 70# 80*90*100)# KK 

BS2 

91 

30 

RATE*RCtIf 1)*AV 

BS2 

92 


GO T(l 110 

BS2 

93 

40 

RATE«RC<I*1 l/f IL )*AV 

BS2 

94 


G) TO 110 

BS2 

95 

50 

RATE >RC U* 1) /TIL )/T IL )«AV 

BS2 

96 


Gil TO 110 

BS2 

97 

60 

RATE>RCII#1)/R00TTAAV 

BS2 

98 

70 

GO TO 110 

BS2 

99 

RATE-RCI It 1)«EXP(RC <I*1)/RRT)«AV 

BS2 

100 

80 

GO TO 110 

BS2 

101 

RATE*RCII»1I*EXP(RCII# 3) /RR T 1 /T I L ) «A V 

BS2 

102 


GO TO 110 

as 2 

103 

90 

RA1E*RC(I*1 )/TfL)/POOTT*AV 

BS2 

104 

100 

GO TO 110 

BS2 

105 

RATE«RC( I* DAEXPIRCII* 3)/RRT)/(TIUA«RC(lt2))AAV 

BS2 

106 

110 

CONTINUE 

BS2 

107 


X«1RRU ) 

BS2 

108 

C 


BS2 

109 

C 

TYPE OF REACTION 

BS2 

110 

r 

« 


BS2 

ill 


00 120 N-l*5 

BS2 

112 

120 

J12345IN)*IRRR(ltN> 

BS2 

113 


GO TO I 150tl60*170»130*140t200t210t220tie0»190)» K 

BS2 

114 

130 

CONTINUE 

BS2 

115 


E«eOATEF(K*3*RRT) 

BS2 

116 


CRR«RATE*RHOOUTIL) 

BS2 

117 


TENRP-CRRARHOUUTIL )«ALPHAUltL)«ALPHAf J2tL ) 

BS2 

118 


TENRN-C RR* ALPHA (J3 * U /F/R/ TIL 1 

BS2 

119 


CALL BOATCC IK t I * L t 3t 3 1 1 .Ot 0 .0 ) 

BS2 

120 

140 

CONTINUE 

BS2 

121 


42«25 

BS2 

122 


E-eOATEFlKt4tRRTI 

BS2 

123 


CRR*RATE*RHOOUTIL lARHOOUTIL )*WTNIX ID 

BS2 

124 


TENRP-CRRAALPHAI JltL) 

BS2 

125 


TENRN«CRR«R«TIL)ARHnOUTIL)*ALPHAf J3tL)«ALPHAI J4tL)/E 

BS2 

126 


CALL BOATCC IKt I * L * 4t 4* -1 .Ot 0.0 ) 

BS2 

127 


GO TO 240 

8S2 

128 

150 

CONTINUE 

BS2 

129 


E«B0ATEFlKt4tRRT) 

BS2 

130 


CRR>RATE*RHOr)UT IL lARMOOUTIL) 

BS2 

131 


TCNRP«CRR«ALPHAUltL)*ALPHA|J2tL) 

BS2 

132 


TEHRH>CRRAALPHAI J3tU* ALPHA IJ4tL)/E 

BS2 

133 


CALL BOATCC 1 Kt 1 1 L t 4» 4 » 1 .0 1 >1 .0 I 

BS2 

134 


GO TO 240 

B3 2 

135 







c 


BS2 

136 

160 

CONTINUE 

BS2 

137 


l-BOATEFfK* 3#RRr 1 

BS2 

138 


CRR*RATCARHOOUTlLMPHnOUTm*WTHIXILI*AV 

BS2 

139 


TEMRP-CRRRRHUOUT III ♦ALPHA U1»L »*ALPHAf JZ#L 1 

BS2 

lAO 


TEHRH»CRR*ALPHA» J3»l I/E/R/TCLI 

BS2 

lAI 


CALL OOATCC (K» 1»L>3*3»1 .OrO.OI 

BS2 

1A2 


GO TO 250 

BS2 

1A3 

C 


BS2 

lAA 

170 

CONTINUE 

BS2 

1A5 


E*BUATEF(K»5»KRT} 

BS2 

1A6 


CRR-RATE*RHOOUnU*PHOOUT(Ll 

BS2 

1A7 


TEnKP>CRR*ALPHAUl»L )♦ALPHAfJ?»L ) 

BS2 

1A8 


TEHRN«CRR*ALPHAU3»LIAALPSA( JA, LI* ALPHA CJ5*L»*RH00UT ILI*R*TIL I/E 

BS2 

1A9 


CALL BOATCC CK» I,L» 5*5»1 .0»-Z.0 1 

BS2 

150 


GO TU 230 

BS2 

151 

C 


BS2 

I5Z 

180 

CONTINUE 

BSZ 

153 


CRR*RATE*RHQOUTILI 

BS2 

15A 


TEPRP«CRR*RHUOUTIL I♦ALPHAU1*LI•ALPHAC J2»LI 

BSZ 

155 


TENRN-O.O 

BSZ 

156 


CALL BOATCC IK» I »L » 3»2> 1 .0*0.0 1 

BSZ 

157 


GO TO 250 

BSZ 

15B 

C 


BSZ 

159 

190 

CONTINUE 

BSZ 

160 


J2-25 

BSZ 

161 


CRR>RATE*RHOOUTILI*RHOnUTf L lAUTNlX ILI 

BSZ 

162 


TENRP«CRR*ALPHAt JlfL 1 

BSZ 

163 


TENRN-0.0 

BSZ 

16A 


CALL BOATCC (K* I *L * A* 1 * 1 .0*0.0 1 

BSZ 

165 


GO TO ZAO 

BSZ 

166 

C 


BSZ 

167 

ZOO 

CONTINUE 

BSZ 

168 


CRR>RATE*RHOOUTILI*RHOOUTfL 1 

BSZ 

169 


TEnRP«CRR*ALPHAIil*U*ALPHAUZ*LI 

BSZ 

170 


TENRN>0.0 

BSZ 

171 


CALL BOATCC IK* I *L • A*2 * 1 .0*0.0 1 

BSZ 

172 


GO TO ZAO 

BSZ 

173 

C 


BSZ 

17A 

210 

CONTINUE 

BSZ 

175 


CRR«RATE*RHOOUTIL l♦»HOrJUT(L l♦WTnIX^LI♦AV 

BSZ 

176 


TEHRP«CRR*RH00UT(L1*ALPHA( J1*LI*ALPHAI J2*LI 

BSZ 

177 


TENRN-0.0 

BSZ 

178 


CALL BOATCC I K* I *L * 1*2* 1 .0*0. 0 1 

BSZ 

179 


GO TO 2 50 

BSZ 

160 




c 


0S2 

181 

Z20 

CONTINUE 

BS2 

182 


C8R>RATE*RHa0UTtLI«8H0nUTfU 

BS2 

183 


TENRP«CRR*ALPIIAU1»LI*AIPMAI J?f L) 

BS2 

I8A 


TENRP-0.0 

BS2 

183 


CAIL BUATCC (K» t t L » 9*2f 1 .0»0.0 1 

BS2 

186 

C 


eS2 

187 

C 

CALCULATt UDUf 

8SZ 

188 

C 


BS2 

189 

230 

WPIJ5I«WP( J3I»TEMRP 

BS2 

190 


J>>MTENRn 

BS2 

191 

2%0 

MPUA)«MP( JAI^TENRP 

BS2 

192 


UNIJA l•WPI JAl ♦TEHRN 

BS2 

193 

2S0 

MPtJ3l>WPCJ3)»TEMRP 

BS2 

19A 


UNf JJl-Wni J3I»TENRP 

BS2 

195 


WPIJ2l«ypl J2I«TENPH 

BS2 

196 


MNU21>4N( J2MTENRP 

BS2 

197 


WPfJl)«WPUll»TENRH 

BS2 

198 


WNIJlt-tfNf JIMTENRP 

BS2 

199 


IE flSWONmi GO TO 260 

BS2 

200 


IF II0UT2.LE.0I GO TO 270 

BS2 

201 


IF f IFINIS.EO.Ot GO TO 260 

BS2 

202 


IF (X.GE.XNAX) GO TO 260 

BS2 

203 


IF tPRNT.GT.PCNTI GO TO 270 

BS2 

206 

260 

NREC»ll-ll*25*I 

BS2 

205 


CALL SFVNV f TENRPf P FRPf NREC I »2 1 

BS2 

206 

270 

CONTINUE 

BS2 

207 

280 

CONTINUE 

BS2 

208 


!r f .N0T.L$NUN(3I 1 GO TO 300 

BS2 

209 


WRITE l6f3AOI T(LI •RHOOUTtLI 

8S2 

210 


DU 200 I-1>NS 

BS2 

211 


NRCC«U-ll*25»l 

BS2 

212 


CALL SFVrV IRPRP.tNRECI»TEHPP»?| 

BS2 

213 


WRITE I6»3A0I ALPHA t I» L 1 » TINRP, TEHRH 

BS2 

21% 

200 

CONTINUE 

8S2 

215 

300 

CONTINUE 

BS2 

216 


DO 310 J*1»NS 

BS2 

217 

310 

WDQTf J»L)«fWP( J)-WNU1 l/PHOOUTfLI/UILI 

BS2 

218 


DO 320 J«1»NS 

BS2 

219 


CALL BOATLP UTKE Y» J*?>SOTf HOT* AX } 

BS2 

220 


HCP12* JI>AX«A9069.31 

BS2 

221 


CALL BOATLP (I TREY* J» A* SOT»HOT» AX 1 

BS2 

222 


HCPtl*JI«AX«A5055.31 

BS2 

223 

320 

CONTINUE 

BS2 

22% 

C 


BS2 

225 



SUBiOUrtNE BOATSl tU 


C 

C 

C 

C 


C 

C 


B0AIS3 • INTEGRAflON ANO SLOB 
BOATCn - BOAT cohhoh 

OIHENSiatl AI90I» AIDt?51* AL0C(»G*6)» ALPHA 5» 90 1 » CGHVI790»3lf 

1 CIM29t?9lf C<»BAR(901» CPTBV(790I* ECCI50I» G(29>» GTBVf790lf 

2 HCPI?»29»> HSrATfSOIf HTAV(750t* IRRfAOl* IRRRIA0»9I» IRTfAOU 

3 ISAVEf6)» J123A9I91* RSEXTI90I# PSII90}* 3X1291* RALPHAt29*90}* 


BS3 
BS3 
BS3 
BS3 
9S3 
8S3 
BS3 
B9 3 
B93 


1 

2 

3 

A 

9 

6 

7 

B 

9 


C 

C 


A RC(A0»3I* RHOI90I* RHIinur(90l» RTI90I* RUI90I* RXEI901* 

RXR(90I* 

BS3 

10 

9 S1GHAI90)* STARTm* fl90l 

• TITLEIIBI* U(90l* 

VEXTIA* 

29* 

21* 

BS3 

11 

6 VJETtA 

»29»2I> M0OT(29*90l* 

HN(29)* 

NPI29I* HTN1XI90I* 

UTNQLE(29I»BS3 

12 

7 XEI90I 

• XEXTI901* XJFTI901 

• XK(90I 

9 

XLEI901* 

RNU(901> 

7(901* 

BS3 

13 

a 70071901* 210(91 








BS3 

lA 











BS3 

19 

LOGICAL 

LHALF»LSHQN 








BS3 

16 











BS3 

17 

EQUIVALENCE Uli 

J123A9I1II 







BS3 

IB 

EQUIVALENCE (ALUC(1*1I* 

CN(l*tll 






BS3 

19 

EQUIVALENCE (ECC* CNU*13lt 

» (TOUT* 

CH(1*19II 




BS3 

20 

EQUIVALENCE ICPTBV* CGMVI 







BS3 

21 

EQUIVALENCE (HCP(l*lt» VNUt) 






BS3 

22 

EQUIVALENCE (START* TTBUI) 







BS3 

23 











9S3 

2A 

CONNON 

TTB(301 

HFI29I 

» 

CRTBV 


GTBV » 

HTBV 



BS3 

29 

CQNNON 

A 

AIO 

9 

ALPHA 


CARS * 

CN 

» 

CPBAR 

»BS3 

26 

1 

CRR 

CVISC 

• 

OELPSI 


DFDL > 

OPDX 

* 

DX 

*BS3 

27 

2 

oxniN 

FOL 

9 

FFF 


G * 

GGG 

* 

HSIAT 

*BS3 

28 

3 

lOELP 

IFCC 

» 

IFINIS 


(OUT * 

lOUTl 

* 

I0UT2 

*BS3 

29 

A 

IPAGE 

IPRESS 

9 

IRR 


IRRR > 

IRT 

» 

(SAVE 

»BS3 

30 

9 

ITFLG 

I»*IS 

9 

J1 


J2 * 

J3 

* 

JA 

• BS3 

31 

6 

J9 

IHALF 

9 

NNOD 


NPSI * 

HXNPT 

* 

NXNPl 

*BS3 

32 

7 

NPSI 

NR 

9 

NRAO 


NRAS » 

NS 

* 

NT 

*BS3 

33 

a 

P 

PCHT 

9 

PRNT 


PRNTXC * 

PSI 

» 

PSIO 

*BS3 

3A 

9 

QX 

RALPHA 

9 

RRU07 


RC * 

RHO 

* 

RHOOUT 

*BS3 

39 

1 

RJ 

RT 

9 

RTACU 


RTJAC » 

RTJOB 

> 

RTNAX 

*BS3 

36 

2 

RU 

R7E 

9 

RXK 


SIGE * 

SIGK 

» 

SIGNA 

*BS3 

37 

3 

T 

TCONT 

9 

TEOGE 


TEHRN * 

TENRP 

• 

TITLE 

*BS3 

38 

A 

TRINE T 

U 

9 

UNIT 


HOOT * 

UN 

> 

HP 

>BS3 

39 

9 

WTNIX 

iitnole 

9 

X 


XCHANG * 

xo 

» 

XE 

*BS3 

AO 

6 

XINIT 

XR 

9 

XR2 


XLE * 

XNAX 

* 

XNU 

*BS3 

At 

7 

7 

ZID 








BS3 

A2 

CONNON 

FL »*:x(iooooi* 

NANAS(3I* 

RPRNIZAOOI* 

SPACRIIOI 


BS3 

A3 











BS3 

AA 

JET/EXT 

FLOW FIELD CONNON SECTION 






BS3 

A9 














c 






BS3 

66 


COKHON OEU » DELE • tJET , lEXf 

• 

IHAXJ 

» 

INAXE 

*BS3 

67 


1 KNAXJ 9 KHAXE • NffJET » fUKXf 

* 

PSJET 

» 

PSEXT 

*BS3 

66 


?. El 9 PI f M 9 P*» 

* 

UST3 

* 

USTE 

*BS3 

69 


3 VJET * VEXT » XJET * XEXT 





BS3 

50 

c 






BS3 

91 


CONNOH ENOCH 





BS3 

52 

c 






BS3 

53 

c 

END OF CONNON FO RE COEIEO TO RESTART FUE 





8S3 

56 

c 






BS3 

55 


CiJHNON FI0t3»5lf IFNANfll* LSWCNI16I 





SS3 

56 

c 






BS3 

57 

c 






BS3 

9B 

c 

DEBUG EX1» EXb» EX7f EX9 





BS3 

59 

c 






BS3 

60 


C1-SEACRC6I 





BS3 

61 


C2*SEACRI71 





BS3 

62 

c 






BS3 

63 


I>1 





BS3 

66 


DEDX>0. 





as 3 

65 


IF IIERESS.EO.OI GO TO 10 





SS3 

66 


IF (RSlfll.lT.PSJETt CALL RnATIl ( 3»X»ESI II )» |NAXJ»KHAXJ» 

1JET*VJETBS3 

67 


lfXJeT»ESJET>TAX»OPOX»TAX»UAX>ESJET»FID»NRJETI 





BS3 

6S 


IF fPSIlII.GE.PSJCTI CALL BOATII 1 X»PSI 1 1 l» INAXE»RMAXE* 

1EXT*VEXTBS3 

69 


l»XEXT»PSEXTf T0X*0P0Xf TBX»UBX»PSJET»FI0»NREX1> 





BS3 

70 

10 

CONTINUE 





BS3 

71 

c 






BS3 

72 

c 

HCP IS WORKING VECTOR OF |H» CPI UP TO NS PAIRS 




BS3 

73 

c 






BS3 

76 


IF (I.EO.ll GO TU BO 





as 3 

75 


IF 1 .NQT.LSWUNfZII GO TO 20 





BS3 

76 


WRITE 16*160) PSllll»AfI)>UII)*RHOf I)*XLEIII* 

SIGNAID* 

CPBARIDfRHOBSl 

77 


lOUTIII 





BS3 

7B 


WRITE 16*1601 iWOnifJ*! )*J«)*NSI 





BS3 

79 


WRITE 16*160) IHCF I2*J)*J«I*NS) 





BS3 

80 


WRITE 16*160) lALPHAU*! )* J-1»NS) 





BS3 

61 

20 

CONTINUE 





BS3 

B2 


EX1«PSIIII*0ELPS1*«2/0X 





BS3 

81 


IX11«.9*|AII)»A| l♦ll) 





a>3 

86 


EX12*.9PlAf 1 I^AI 1-1)1 





BS3 

89 

c 






BS3 

86 

c 

INTEGRATE NUNENTUN EOUATIHN 





BS3 

87 

c 






BS3 

88 


PUfI )*tEXll*f llI*l)-llf I ))»E«l?»fUII-l)-UII) ))/EXl*UIII 



BS3 

89 


RUII l*PUf I)-OXAOPDX/RHOI ll/UI I)*RBUOTPOX*tRIIOINPSI l-RHOIS 

))*32. 

17AB33 

90 






i/RHom/um 

TERHl»fEXll*CUfl*ll-linn*E)a2*IU(I-l)-UIIIII/EXl 8S3 

TERH?««»BUOY*OX*IRHO<i1PSI I-RHOII » »♦ 32 . 1 74/RHOI II/lMl » BS3 

IF USMONC13)) WRITE I6»170l X » U T ERNl » TERM? BS 3 

EX3-0.0 ®53 

EX-i-0.0 9S3 

DO 30 ®53 

IX3»IX3*HCPU*JI*WDPT(J»I) BS 3 

EX<^«EX^*HCP(2» J I* t ALPHA( J» 1*1 »- ALPHA* 4 1 1-l » I BS3 

EX2«EX1#CPBARI n BS 3 

IX5-XLim*A(I WSIGKAfl I BS 3 

EX6«.5A*EX5*XLE(I»l)*ACI*n/SIGMA< I*l) ) BS3 

EX7«.5*CEX5*XLEn-l >*AI I-n /SIGH AC I-II » BS3 

IXB-CPBARC II^ACIl/SIGMACn BS3 

EX9-.5* CEXB^CPBARC I*n*An*l»/SIGKAII*l)) BS3 

IX10«.5*CEXB*CPBAP CI-l )*A* I-U/SIGMAC I-IM 8S3 

CX14-EX4«EX3/4.0 BS3 

BS3 

INTEGRATE ENERGY EQUATION 8S3 

BS3 

RTin«IUC I*1I-UC I- 1 ))** 2 *ACII/EX 2 /^. 0 *DX*DP 0 X/RHUC n/CPBARI II»TCIIBS3 
l*CIEXV*EXl^l*TCIUI*(EX10-EXl<H*TCI-ll-CEX9+EXl0»*Tn II/EX2-EX3A0XBS3 
2/CPBARtU BS3 

RHCUIX-OX/IRHUOUTC I >*UCI n BS3 

BS3 

INTEGRATE SPECIES TOUATIOMS BS3 

BS3 

IF fLSU0Nl2n WRITE (6,160) t OX ( J ) , 4«l ,NS ) 8S3 

DO AO J«1,NS BS3 

OX CJ l■IEX6♦C ALPHA* 4, I*1)-ALPHA( 4, I ) I *EX ?♦ C ALPHA* 4, I-l I -ALPHA* 4, 1 1 IBS 3 
I |/EX1*ALPHA*4,1U0X *4)*RH0UIX BS3 

DO 50 M-1,NS BS3 

00 50 N«1,NS BS3 

CH*H#**l»CH*H,N)*RHnUIX BS3 

IF CN.EO.NI Cn*N,N)»CM*M,M)+1.0 BS3 

CONTINUE BS3 

CALL BOATSL *aX,CH,NS) BS3 

BS3 

BS3 

DO 60 4«1,NS BS3 

RALPHA*4,1)>QX*4) BS3 

BS3 

integrate TKE EOUATIOMS BS3 

BS3 

ic f;;i Tr 70 BS3 


integrate TKE EOUATIOMS 
IF CIVIS.GE.O) GU TO 70 



. j 




70 

C 


80 


C 

C 

C 


90 


100 


C 

C 


C 

C 

C 

110 

C 


C 

C 


RXK(n«XKl n»(EXll*(XK tl^l (•XKini ♦i:X12«(XK(I-l >-XKf II 1 1 /E XI / SI6K*B$ 3 
Iiun«ll-U(l-1> )*«2*A( n/EXlM.-XU n/UlIMOX BS3 

RXEf I l>XE (I !♦ tEXll« IXE (141 |-.xb I n I^EXieAlXE n-t l-KE (1 1 1 1 /EXI/S 1GE4BS 3 
iCl*IUIl4n-un-l))««?*XF(I)*A(l l/EXl/XKI II/A.-C2 4XE ( 1 1 /U ( 1 1 / XK (BS3 


2II*DX 
CONTINUE 


60 TO 150 
CONTINUE 

IE (PSI (1) .GT.O. I GO Tn 130 
EX3«^,0AXnU( n*DX/DFLPSI/OELPS I 
RH0UIX«0X/ (RHOOUTI J 1 40(1 II 


COHPUTE U AT CENTER LINE 


BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 


RU(1I«EX3*(U(2I>UI 1 I14U(1 I-0X*DPDX/RH0(1 l/U(l I 4R BUQY«>OX* (RHO (NPS 1 1 B$ 3 
1-RHOCl I l*32,17A/RH0(l )/U(l I BS3 

TERH3«EX3*(U(2I-U( 111 BS3 

IERNA-RBU0V*0X*|RHn(*1P5I l-RHO(l ll♦32.17A/RH0(l l/Ulll BS 3 

IF (LSW0N(13II WRITE <6>170l Xt I»TERN3> TERNA BS3 


EXA-0.0 
00 90 J«lfNS 

EXA-EXA4HCPI1. JI*W00T(J»1I 
RALPHAf J#NPSI )-ALPHA(J»NPSII 


8S3 

BS3 

BS3 

BS3 


0X(JI-EX3AXLEm*(ALPHAU,2l-ALPHAU,lI I/SIGHA(II4ALPHAC J,ll40X(JIBS3 


l*RMaUIX 

00 100 n-i»Ns 

00 100 N«1»NS 
CN(N»NI«CH(N»NI«RHOUIX 
IF (N.EO.NI Cn(NfNI«CN(N>N|4l.o 
CONTINUE 

CALL BOATSL (OXfCN»NSI 


00 110 J-I»NS 


COHPUTE SPECIES AT CENTER LINE 


RALPHA( Jf 1 |>OX(Jt 


CALCULATE TEMP. AT CENTER LINE 


BS3 

BS3 

BS3 

BS3 

B$3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 

BS3 


RTI1I-EX3P (T(2I-T(1 ll/SIGMA(l |4T(1 |40X*DP0X /RHO( 1 1 /CP9AR (1 |-EXA«0XBS3 
l/CPBARCll BS3 

BS3 


136 

137 
I3B 
139 
lAO 
lAl 
1A2 
1A3 
lAA 
1A5 
1 A6 
IA7 
lAB 
1A9 

150 

151 

152 

153 
15A 

155 

156 

157 

158 

159 

160 
161 
162 
163 
16A 

165 

166 
167 
16B 

169 

170 

171 

172 

173 

174 

175 

176 

177 
17B 

179 

180 







C 


t * ‘ ' n 


c 

IHTEGPATE TKE EQUATIONS AT CENTERLINE 

BS3 

181 

c 


DS3 

182 


IF (IVIS.GE.OI GO TO 120 

BS3 

163 


RXK( 1)«XKI1MEX3«I XK(?)>XKf 11 ) /S IGK-XE ( 1 ) /U( 1 1 *0X 

BS3 

186 


IF (XKI1I.lt. l.E-20) XK(U*l.E-20 

BS3 

185 


RXEm«XE(l»^EX3A(XE(2)-XE (1M/SIGE-C2*XE(1 )*A2/UC1 )/XK( 1 »*DX 

BS3 

186 


DUNK-RXKtl )/RXK(2> 

BS3 

187 


IF I0UNK,LT.,1» RXK (1 »«PXK(2) 

BS3 

168 


DUnE-RXE(ll/RXE(2l 

BS3 

189 


IF (DUNE.LT..1) RXE (1 )«RXE (2) 

BS3 

190 

120 

CONTINUE 

BS3 

191 

C 


BS3 

192 


GO TO 150 

8$ 3 

193 

130 

CONTINUE 

BS3 

196 


on 160 J«1»NS 

BS3 

195 


RALPHA( J»1 )«ALPHAU*1I 

BS3 

196 


RALPHAIJfNPSn •ALPHA IJ.HPSn 

BS3 

197 

1<^0 

CONTINUE 

BS3 

190 


IF (IVIS.GE.OI GU TP 150 

BS3 

199 


RXK(ll-0. 

BS3 

^00 


RXE(1I*0. 

BS3 

201 


PXKdIPSI l>0. 

BS3 

202 


RXE(HPSn«0. 

6S3 

203 

150 

CONTINUE 

BS3 

206 


RETURN 

BS3 

205 

C 


BS3 

206 

160 

FORNAT (IH »10E13.5) 

BS3 

207 

170 

FORNAT (IH »E13.5*I6>2E13.5) 

BS3 

208 


END 

BS3 

209 





■■■ , "■ 


SUBROUTINE BOATTK I T» T TB # I TKE Y» SDT #HOT »NT I 
DIMENSION TTBOOl 
NT1«NT-1 
DO 10 IT»l,MTl 
DT-TTBfIT*l»-TTB(m 
SDT« fT-TTB(IT) »/DT 
HDT-ITTBI IT+ll-TI/OT 
IF riSDT*H0T».6E.0.0) GO TO 
10 CONTINUE 

WRITE (6»30) T»IT 

ITKEY-O 

RETURN 

20 ITKEY-IT 

RETURN 
C 

30 FORMAT CIH ,28H TEMPERATURE OUT OF RANGE »EU.5»I5) 

END 


BTK 1 

BTK 2 

BTK 3 

BTK A 

BTK b 

BTK 6 

BTK 7 

BTK B 

BTK 9 

BTK 10 

BTK 11 

BTK 12 

BTK 13 

BTK 19 

BTK 19 

BTK 16 

BTK 17 








'ItliaimJlirt.iniii. 




Ilppp^pi 


WIWIIH'I : (!ipii|((!igifl|jt|fji|ll';iiHp^^ 



c 

c 

c 

c 


c 

c 


SUBROUTINE BOATVl 
TURBULENT VISCOSITY R0UTIf4E 
ROATCM - BOAT COMNON 

DINENSION A(!»0)» AID(Z9)f ALCCIS0»6I* ALPHA (25* 50 )» CGHV(750*3)* 

1 Cni23pZj}» CPBAR(50)* CPTBV(750)* ECCI50I* G(25>* GTBV(750I* 

2 HCP(2*25)* HSIAT(50I» HTBV(750)* IRR(<*0I* IRRR(<iO»5}* IRT(^OI* 

3 ISAVE(61. J123<i5(5l. PSEXT(50»* PSI(50»* 0X(25I# RALPHA ( 25* 50) * 

A RC(<kO*3)» RHfJ(50)* RH0'JUTI50I* RT(50)* RUfSOI* RXtCSOI* RXKI50)* 
5 SIGPA(50)» STARTin* T(50)* TITLE(IB). U(50)* VEXT ( A* 25*2 ) * 

b VJET(R»25*2 ). WDnT(25*50l» WNC25)* WP(25)* WTNIX(50I* 

7 XE(50)f XEXT(50I* XJET(50). XK(50»* XLE(50)» XNUC50)* 

8 Y0UT(50)» ZI0(5» 

LOGICAL LHALF*LSW0N 


CQUIVALLNCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


(Jl* J123A5(in 
(ALUC(l*n* CN(l,in 
(ECC* CP(1*13)I* (YUUT* CM(1*15)) 
(CPTBV* CGHV) 

(HCP(1*D* WN(in 
(START* TTB(l)) 


BVI 1 

BVI 2 

BVI 3 

BVI A 

BVI 5 

BVI 6 

BVI 7 

BVI a 

BVI 9 

BVI 10 

BVI 11 

WTN0LE(25I*BVI 12 

YI50I* BVI 13 

BVI lA 

BVI 15 

BVI 16 

BVI 17 

BVI IB 

BVI 19 

BVI 20 

BVI 21 

BVI 22 

BVI 23 

BVI 2A 


C 

c 


COMMON 

TTB(30I 

HF (25) 

CPTBV 

* 

GTBV 

* 

HT8V 



BVI 

25 

COMMON 

A 

AID 

ALPHA 

* 

CARB 

* 

CM 

> 

CPBAR 

*BVI 

26 

1 

CRR 

CVISC 

OELPSI 

* 

DFDL 

» 

DPDX 

* 

OX 

*8VI 

27 

2 

OXMIN 

FDL 

FFF 

* 

G 

* 

GGG 

* 

MSTAT 

*BVI 

26 

3 

lOFlP 

lECC 

IFINIS 

* 

lUUT 

* 

lOUTl 

* 

I0UT2 

*BV1 

29 

A 

IPAGE 

i PRESS 

IRR 

* 

IRRR 

* 

IRT 

» 

ISAVE 

*BVI 

30 

5 

ITFLG 

IVIS 

Jl 

* 

J2 
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» PRNT 

» 

PRNTXC 


PS I 

* 

PSIO 

*BEF 

9 

QX 

RAIFHA 

• RNUOY 

» 

RC 


fho 

> 

RHOOUT 

»BEF 

1 

FJ 

FT 

> PfACU 

* 

RTJAC 


RTiOB 

» 

RTNAX 

*BEF 

2 

RU 

RXE 

* RXR 

» 

SIGE 


SIGN 

* 

SIGN A 

• BEF 

3 

1 

TCTNT 

> TFOGE 

» 

fSNRN 


TENRP 

» 

HUE 

• BEF 

4 

tkinet 

U 

* UNlf 

» 

WDOT 


UN 

# 

UP 

• BEF 

5 

UfNIX 

UTNOLE 

* X 

* 

XCHANG 


XO 

* 

XE 

• BEF 

b 

XINIT 

XK 

* XR2 

» 

XLE 


XNAX 

» 

XNU 

• BEF 

7 

7 

210 








BEF 

CONNON 

FlJTEXflOOQOl* NARASI31* 

RPRN 124001 

* 

SPACRf 101 


BEF 











BEF 

JET/EXf 

FLOW FIELD CONNCH SFCTION 







BEF 


1 

2 

3 

5 

5 

6 
7 
B 
0 

10 

11 

12 

13 

lA 

15 

lb 

17 

IB 

19 

20 
21 
22 
23 
2% 
25 
2b 
27 
2B 

29 

30 

31 

32 

33 
3A 
35 
3b 
37 
3B 

39 

40 

41 

42 

43 
4% 
45 


c 


fFI 





C 

C 

c 

c 

c 


to 

c 


20 

30 

%0 


CONnUH 


1 

2 
1 


OELJ 

KHAIJ 

P\ 

VJU 


out 

irnAltF 

PP 

VEKT 


IJtf 

NI»JI:T 

P^ 

XJE 1 


» lEXf 
f nrext 

f XEXf 


COfMION ENOCH 

END OF CONHON TO BF COPIED TO RESTART FILE 
CONNON FIDI 3 »SI» |FNAN« 3 I» LSM 0 NI 16 I 


• INAXJ 
, FSJET 
t USTJ 


F* 0.0 
SIGN-t.O 
03 10 N*lfNJ 

IF IN.FO. 2 .AH 0 .R.E 0 . 9 ) CO TO 10 

IF IN.GT. 2 I SICH--I .0 

J«J 123 ASt*ll 

E-SICN* 6 t J 1 »E 

CONTINUE 

E-E/RRT 

IF lABSIEI .IT.nO.I CO TO 30 

IF lE.lT.O.OI GO TP 20 

IF lE.GT.O.OI E«EXP 180 . 0 I 

GO TO AO 

E>EXR(oB0.01 

GO TO AO 

E«EXF(EI 

BOATEFi^E 

RETURN 

END 




BEF AB 
• IHAXE *BEF AT 

» PSEXI »BEF AS 

, USTE »BEF A 9 

BEF 50 
BEF 51 
BEF 52 
BEF 53 
BEF 5 A 
BEF 55 
BEF 96 
BEF 57 
BEF 50 
BEF 59 
BEF 60 
BEF 61 
BEF 62 
BEF 63 
BEF 6 A 
BEF 65 
BEF 66 
BEF 67 
BEF 6 B 
BEF 69 
BEF 70 
BEF 71 
BEF 72 
BEF 73 
BEF 7 A 
BEF 75 
BEF 76 
BEF 77 
BEF 7 B 





I 

I 

i 


! 

j FUNCTION ETINF (STI*'E» 

I C ELAPSED TlHt FUNCTION 

I ^ ETIME ■ SECOND ISTINEI/60. 

I RETURN 

END 


ij 




ETI I 

ETI 2 

ETI 3 

tn A 

ETI 5 

ETI b 

ETI 7 



FUNCTION NFMAB lA* Bl 
C 

C NFUAB - NUMBER OF STANDARD FORTRAN WORDS FROM A TO B 
C 

NFWAB • LOCFCBI - LOCF(A) 

RETURN 

END 


NFW 1 

NFW 2 

NFW 3 

NFW A 

NFW 5 

NFW 6 

NFW 7 



f 





C 

C 

c 

c 

c 


c 

c 


c 


c 


SUBRaUTINE PBFOR (DUHHY* nro> NW0> VEC) 

DUHHY DISK READ ROUTINE - NAY BE REPLACED BY 
REAL DISK READ ROUTINE. 


1 

2 

3 

A 

5 

6 

7 

8 


BCATCH - BOAT CCIMHON 

DIHENSION A(50»> AID«25»» ALnC*50»6>» ALPHA (25. 50l» CGHVI750.3I. 
CM(25.25». CPBAR(50». CPTBV(750>» ECCtSOJ. G(?51» GTBV(750I» 
HCP(2»25>. H3TAT(50I. HTBV(750I» IRR(AO». IRRRCA0.5I. IRTIAOI. 
ISAVEI6). J123A5I5J. PSEXT(50I. PSK50I. QX(25». R AL PHA (25. 50 » » 
RC(A0»3). RH0(50I> RH0nUT(50J» RT(50». RU(50)» RX6(50I» RXK(50)» 
SlGnA(50). STARTdI. T(50>. TITLEdBl. U(50l» VEX T ( A . 25 » 2 I . 
VJET(A.25*2). W00T(?5.50I* WN(25I. WP(25». WTHIXI50I# HTMOLE(25> 
XE(50I. XEXTI50). XJFT(50>. XK(501. XLE(50I» XMU(50I» Y(50). 
Y0UT(501. Z13(5) 


LOGICAL LHALF.LSWON 


EQUIVALENCE 

equivalence 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


CONNON 

CONNON 


(Jl. J123A5(in 
(ALOCd.d. CNd.ld 
(ECC. CHd.nn. (TOUT. 
(CPTBV. CGHVI 
(HCP(l.l). WNdd 
(START. TTBdd 


CIK1.15M 


CUNHON 


TTB(30). 

HF(25) 

p 

CPTBV 

P 

GTBV 

HTBV 


A . 

AIO 

p 

ALPHA 

P 

CARB 

CN 

. 

CRR . 

CVISC 

p 

OELPSI 

P 

DFOL 

DPOX 

. 

OXNIH . 

FOL 

p 

FFF 

P 

G 

GGG 

. 

lOELP . 

lECC 

p 

IFINIS 

P 

lOUT 

lOUTl 

. 

IPAGE . 

IPRESS 

p 

IRR 

P 

IRRR 

IRT 

. 

ITFLG » 

IVIS 

p 

Jl 

P 

J2 

J3 

. 

J5 . 

lhalf 

p 

NNOD 

P 

NPSI 

nXNPT 

. 

NPSI . 

NR 

p 

NRAO 

P 

NRAS 

NS 

. 

P . 

PCNT 

p 

PRNT 

P 

PRNTXC 

PS I 

. 

OX . 

PALPHA 

p 

PBUOY 

P 

RC 

RHO 

. 

RJ . 

PT 

p 

RTACU 

P 

RTJAC 

PTJOB 

. 

RU . 

RXE 

p 

RXK 

P 

SIGE 

SIGK 

. 

T . 

TCONT 

p 

TEDGE 

P 

TENRN 

TENRP 

. 

TKINET . 

U 

p 

UNIT 

P 

WDOT 

WN 

. 

WTNIX . 

WT^HLE 

p 

X 

P 

XCHANG 

XD 

. 

XINIT . 

XK 

p 

XK2 

P 

XLE 

XNAX 

. 

Y . 

FL JTEXd 

zn 

0000). 

NAHASO). 

RPRN(2A00) . 

SPACR( 10 


CPBAR 

OX 

HSTAT 

I0UT2 

ISAVE 

JA 

NXNPl 

NT 

PSIO 

RHOOUT 

RTNAX 

S IGNA 

TITLE 

WP 

XE 

xnu 


POR 

PJR 

POR 

POR 

PDR 

POR 

POR 

PDR 

POR 

PDR 

POR 

POR 

.PDR 

PDR 

POR 

POR 

POR 

POR 

PDR 

PDR 

POR 

POR 

PDR 

POR 

PDR 

POR 

.POR 

.POR 

.POR 

.POR 

.POR 

»POR 

.POR 

.PDR 

.POR 

»POR 

.PDR 

»POR 

.POR 

.PDR 

*POR 

.POR 

POR 

POR 

POR 


1 

2 

3 

A 

5 

6 

7 

8 
V 

10 

11 

12 
13 
lA 

15 

16 

17 

18 
IQ 
20 
21 
22 
23 
2A 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3A 

35 

36 

37 

38 

39 
AO 
A1 
A2 
A3 
AA 
A5 
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T’lpT 


JET/tKT FLOW FIELD COMMON SECTION 

COMMON OELJ # DELE » IJET » lEKT 

1 KNAXJ $ KMAXE * N»JET > NREXT 

2 PI » P? » P3 » PA 

3 VJET , VFXT , XJET , XEXT 

COMMON ENOCM 

EMO OF COMMON TO BE COPIED TO RESTART FILE 
COMMON FID(3»S)» TFNAMI3|» LSWLN(16I 


CALL SFVMV IFLJTEXtNPO)* VEC* NWOI 

NRO « NRO-»NWO 

RETURN 

END 




" ' """ '' ■'H,-' :h,li 




POR 

A6 



POR 

A7 

IMAXJ » 

INAXE 

»POR 

A8 

PSJET p 

PSEXT 

#POR 

A9 

USTJ » 

USTE 

»PDR 

SO 



POR 

51 



POR 

52 



POR 

53 



POR 

5A 



POR 

55 



POR 

56 



POR 

57 



POR 

58 



POR 

59 



POR 

60 



POR 

61 



POR 

62 



POR 

63 







'^1 





IIH'MIIB 




SUBRUUTINE PBFOW IDUNMY* MROt NWOt VEC I 


]i!«'"”pwiiiif ;:r 


C 

c 

c 

c 

c 


c 

c 


DUMY DISK WRITE KPUTINF 
REAL DISK WRITE ROUTINE. 


NAY BE REPLACED BY 


BOATCN - BOAT CONNON 

DIHENSION Af50l» AIO(25l> ALnC(50#6># ALPHA ( Z5* 50 > » CGHVt750»3)» 

1 CN(25»25)» CPBARI50I» CPTBVI750I# ECCC50)# GI25I» GT8VC750»» 

2 HCP<2»25)f HSTATI50lf HT9VI750I# IRPI<iOJ» IRRRf^0>9l» IRTCAO)# 

3 ISAVEf6)» J12345(5)» PSEXTf50)» PSlf50l> QXf25l# RALPHA 125* 50 1> 
RCIA0»3)» RHUt50l» RHn0UTI50l> RTI50)* RUI50I* RXEI50I* RXKl50t> 

5 SIGNAI50I* STARTin* T«50)f TITLE(1B)> Ut50)» VEXT t5»25>2 

6 VJETf5»25*2l» WU0T(25»50I» WH(25)» WPI25I* WTHIXI50U WTHOLEC25) 

7 XE(50I» XEXTI50I» XJET(50)» XKI50I» XLEC50W XNUC50I> Yf50»» 

8 Y0UTI50)* ZI0I5I 

LOGICAL LHAlFtLSWOr 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


CONNQN 

CUNNON 


fJl* J12355<in 
IALOC(l»l»f CNn*U) 
IECC» CNIltlBUf CYQUTf 
ICPTBV. CGHVI 
CHCPii*n> wNun 
(START* TTBtin 


CN(l*15n 


POW 

POW 

POW 

ROW 

PDW 

POW 

POW 

PDW 

POW 

PDW 

POW 

POW 

*POW 

POW 

POW 

POW 

POW 

PDW 

POW 

POW 

POW 

POW 

POW 

POW 

POW 


1 

2 

3 

5 

6 
7 
6 
9 

10 

11 

12 

13 

1« 

15 

16 
17 
IB 

19 

20 
21 
22 
23 
2^ 
25 


cnnnoN 


TTB( 30> 

HFI25I 

P 

CPTBV 

9 

GTBV 

» 

HTBV 



POW 

26 

A 

AID 

P 

ALPHA 

9 

CARS 

f 

CN 

* 

CPBAR 

*PDW 

27 

CRR 

CVISC 

0 

DELPSI 

9 

OFOL 

» 

DP OX 

* 

OX 

*POW 

2B 

OXNIN 

FDL 

P 

FFF 

P 

G 

t 

GGG 

p 

HSTAT 

*PDW 

29 

lOELP 

IFCC 

P 

IFINIS 

9 

lOUT 

* 

lOUTl 

p 

I0UT2 

»POW 

30 

IPAGE 

IPRESS 

P 

IRK 

9 

IRRR 

f 

IRT 

p 

ISAVE 

»PDW 

31 

ITFLG 

IVIS 

P 

J1 

9 

J2 

t 

J3 

p 


*POW 

32 

J5 

lhalf 

P 

NMOO 

9 

rtPSI 

* 

NXMPT 

p 

NXNPl 

*POW 

33 

MPSI 

NP 

9 

NRAO 

9 

NRAS 

t 

NS 

p 

NT 

»POW 

35 

P 

PCMT 

9 

PRNT 

9 

PRNTXC 

* 

PSI 

p 

PSIO 

»POW 

35 

OX 

R ALPHA 

P 

RBUOY 

9 

RC 

f 

RHO 

p 

RHOOUT 

*POW 

36 

RJ 

RT 

P 

rtacu 

9 

RTJAC 

* 

RTJ09 

p 

RTNAX 

*POW 

37 

RU 

PXE 

P 

RXK 

9 

SIGE 

* 

SIGK 

p 

SIGHA 

f POW 

38 

T 

TCONT 

9 

TEDGE 

9 

TENRN 

* 

TENRP 

p 

TITLE 

fPOW 

39 

TKINET 

U 

9 

UNIT 

9 

WOOT 

* 

WN 

p 

WP 

»POW 

50 

WTNIX 

WTNniF 

9 

X 

9 

XCHANG 

* 

XO 

p 

XE 

#POW 

51 

XINIT 

XK 

9 

XK? 

9 

XLE 

p 

XNAX 

p 

XMU 

*POW 

52 

Y 

710 









POW 

53 

FLJTEXClOOOOIf 

NANAS(3)» 

RPRNI2500I 

p 

SPACR( 101 


POW 

55 


PDW 


•*5 


I 
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JET/FXT Flow field common SFCTION 

COMMON DELJ * DELE » * NREXT 

1 KMAXJ , «»MAXF , NRJET » NREXT 

\ PI » F? » , 

3 VJET » VFXT t XJET » XEXT 

COMMON ENOCH 

end OF COMMON TO BE COPIED TO RESTART FILE 
COMMON F10f3»5l, IFNANI3I. LSWONlie) 


CALL SFVMV IVEC» F L JTE X CNRO I , NWOI 
NRO • NRO^NNO 

return 

END 



r> o o 


1 


SUBROUTINE SEVFL IVAL» VEC* NW) 

SFVFl - STANDARD FORTRAN VECTOR FILL 

OINENSION VECm 
DO 10 I • 1* MW 
VECIII • VAL 
10 CONTINUE 
RETURN 
END 


FVF 1 

FVF 2 

FVF 3 

FVF A 

FVF 5 

FVF 6 

FVF 7 

FVF 8 

FVF 9 

FVF 10 



'I 


r> o n 


• 


SUBROUTINE SFVHV JVEClf VEC?, IIWI 

SFVNV - STANDARD FORTRAN VECTOR HOVE 

DIRENSION VEClIll* VEC?fl) 

DO 10 I • I, NU 
VECZIll • VEClin 
10 CONTINUE 
RETURN 
END 


FVN 1 

FVN 2 

FVN 3 

FVN A 

FVN 5 

FVN 6 

FVN 7 

FVN 8 

FVN 9 

FVN 10 


